Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “distributed simulation”

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 613 records · Page 34

Simulation–Extrapolation for Bias Correction with Exposure Uncertainty in Radiation Risk Analysis Utilizing Grouped Data

Summary In observational epidemiological studies, the exposure that is received by an individual often cannot be precisely observed, resulting in measurement error, and a common approach to dealing with measurement error is regression calibration (RC). Use of RC, which requires assumptions about the distribution of unknown error-free (true) variables, leads to concern about the possibility of bias due to misspecification of that distribution. The simulation–extrapolation (SIMEX) method, in contrast, does not require a distributional assumption. However, analyses of large cohorts may be performed by using grouped or person-year data, and application of SIMEX to grouped data is not straightforward, particularly when there is a mixture of classical and Berkson measurement errors. We compared RC and SIMEX with grouped data analyses to assess robustness of the RC method to misspecification of the true dose distribution. We also applied SIMEX assuming mixtures of classical and Berkson errors and compared the results with those obtained by using RC for classical error only. SIMEX had less bias than RC and performed well regardless of the true dose distribution, whereas RC based on a misspecified true dose distribution showed greater bias than when based on the correctly specified true dose distribution.

Misumi, Munechika↗

SMA OBSERVATIONS OF CLASS 0 PROTOSTARS: A HIGH ANGULAR RESOLUTION SURVEY OF PROTOSTELLAR BINARY SYSTEMS

We present high angular resolution 1.3 mm and 850 {mu}m dust continuum data obtained with the Submillimeter Array toward 33 Class 0 protostars in nearby clouds (distance < 500 pc), which represents so far the largest survey toward protostellar binary/multiple systems. The median angular resolution in the survey is 2.''5, while the median linear resolution is approximately 600 AU. Compact dust continuum emission is observed from all sources in the sample. Twenty-one sources in the sample show signatures of binarity/multiplicity, with separations ranging from 50 AU to 5000 AU. The numbers of singles, binaries, triples, and quadruples in the sample are 12, 14, 5, and 2, respectively. The derived multiplicity frequency (MF) and companion star fraction (CSF) for Class 0 protostars are 0.64 {+-} 0.08 and 0.91 {+-} 0.05, respectively, with no correction for completeness. The derived MF and CSF in this survey are approximately two times higher than the values found in the binary surveys toward Class I young stellar objects, and approximately three (for MF) and four (for CSF) times larger than the values found among main-sequence stars, with a similar range of separations. Furthermore, the observed fraction of high-order multiple systems to binary systems in Class 0 protostars (0.50 {+-} 0.09) is also larger than the fractions found in Class I young stellar objects (0.31 {+-} 0.07) and main-sequence stars ({<=}0.2). These results suggest that binary properties evolve as protostars evolve, as predicted by numerical simulations. The distribution of separations for Class 0 protostellar binary/multiple systems shows a general trend in which CSF increases with decreasing companion separation. We find that 67% {+-} 8% of the protobinary systems have circumstellar mass ratios below 0.5, implying that unequal-mass systems are preferred in the process of binary star formation. We suggest an empirical sequential fragmentation picture for binary star formation, based on this work and existing lower resolution single-dish observations.

79 ASTRONOMY AND ASTROPHYSICS↗

Self-consistent field model for strong electrostatic correlations and inhomogeneous dielectric media

Electrostatic correlations and variable permittivity of electrolytes are essential for exploring many chemical and physical properties of interfaces in aqueous solutions. We propose a continuum electrostatic model for the treatment of these effects in the framework of the self-consistent field theory. The model incorporates a space- or field-dependent dielectric permittivity and an excluded ion-size effect for the correlation energy. This results in a self-energy modified Poisson-Nernst-Planck or Poisson-Boltzmann equation together with state equations for the self energy and the dielectric function. We show that the ionic size is of significant importance in predicting a finite self energy for an ion in an inhomogeneous medium. Asymptotic approximation is proposed for the solution of a generalized Debye-Hückel equation, which has been shown to capture the ionic correlation and dielectric self energy. Through simulating ionic distribution surrounding a macroion, the modified self-consistent field model is shown to agree with particle-based Monte Carlo simulations. Numerical results for symmetric and asymmetric electrolytes demonstrate that the model is able to predict the charge inversion at high correlation regime in the presence of multivalent interfacial ions which is beyond the mean-field theory and also show strong effect to double layer structure due to the space- or field-dependent dielectric permittivity.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Assessment of grid-friendly collective optimization framework for distributed energy resources

Distributed energy resources have the potential to provide services to facilities and buildings at lower cost and environmental impact in comparison to traditional electric-gridonly services. The reduced cost could result from a combination of higher system efficiency and exploitation of electricity tariff structures. Traditionally, electricity tariffs are designed to encourage the use of ‘off peak’ power and discourage the use of ‘onpeak’ power, although recent developments in renewable energy resources and distributed generation systems (such as their increasing levels of penetration and their increased controllability) are resulting in pressures to adopt tariffs of increasing complexity. Independently of the tariff structure, more or less sophisticated methods exist that allow distributed energy resources to take advantage of such tariffs, ranging from simple pre-planned schedules to Software-as-a-Service schedule optimization tools. However, as the penetration of distributed energy resources increases, there is an increasing chance of a ‘tragedy of the commons’ mechanism taking place, where taking advantage of tariffs for local benefit can ultimately result in degradation of service and higher energy costs for all. In this work, we use a scheduling optimization tool, in combination with a power distribution system simulator, to investigate techniques that could mitigate the deleterious effect of ‘selfish’ optimization, so that the high-penetration use of distributed energy resources to reduce operating costs remains advantageous while the quality of service and overall energy cost to the community is not affected.

Pensini, Alessandro↗

A temperature-dependent coarse-grained model for the thermoresponsive polymer poly(N-isopropylacrylamide)

A coarse-grained (CG) model is developed for the thermoresponsive polymer poly(N-isopropylacrylamide) (PNIPAM), using a hybrid top-down and bottom-up approach. Nonbonded parameters are fit to experimental thermodynamic data following the procedures of the SDK (Shinoda, DeVane, and Klein) CG force field, with minor adjustments to provide better agreement with radial distribution functions from atomistic simulations. Bonded parameters are fit to probability distributions from atomistic simulations using multi-centered Gaussian-based potentials. The temperature-dependent potentials derived for the PNIPAM CG model in this work properly capture the coil–globule transition of PNIPAM single chains and yield a chain-length dependence consistent with atomistic simulations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Effects of Non-Sphericity on Mineral Dust Modeling

The dependency of non-sphericity on gravitational settling of mineral dust particles is parameterized for prolate ellipsoids and Reynolds number lower than 2. The settling speed is numerically solved from the momentum equation as a function of particle diameter and aspect ratio. The reduction of settling speed due to non-sphericity is included in the GO-CART model to simulate dust size distribution for April 2001. Two numerical schemes for solving sedimentation are compared. For particles of diameter greater than 5 micron the simulated size distribution is sensitive to the numerical sedimentation scheme. Changing the particle shape from spherical to non-spherical with lambda=2, makes little difference to the simulated surface concentration and size distribution except at the periphery of the dust sources. However, when very elongated particles (lambda=5) are simulated the differences between non-spherical and spherical particles are significant. With limited in-situ measurements reporting most frequent lambda around 1.5, the overall effects on global modeling is rather negligible and the essential benefit is to relax the CFL condition of Eulerian settling schemes.

Ginoux, Paul↗

Development and Evaluation of Chemistry‐Aerosol‐Climate Model CAM5‐Chem‐MAM7‐MOSAIC: Global Atmospheric Distribution and Radiative Effects of Nitrate Aerosol

Abstract An advanced aerosol treatment, with a focus on semivolatile nitrate formation, is introduced into the Community Atmosphere Model version 5 with interactive chemistry (CAM5‐chem) by coupling the Model for Simulating Aerosol Interactions and Chemistry (MOSAIC) with the 7‐mode Modal Aerosol Module (MAM7). An important feature of MOSAIC is dynamic partitioning of all condensable gases to the different fine and coarse mode aerosols, as governed by mode‐resolved thermodynamics and heterogeneous chemical reactions. Applied in the free‐running mode from 1995 to 2005 with prescribed historical climatological conditions, the model simulates global distributions of sulfate, nitrate, and ammonium in good agreement with observations and previous studies. Inclusion of nitrate resulted in ∼10% higher global average accumulation mode number concentrations, indicating enhanced growth of Aitken mode aerosols from nitrate formation. While the simulated accumulation mode nitrate burdens are high over the anthropogenic source regions, the sea‐salt and dust modes respectively constitute about 74% and 17% of the annual global average nitrate burden. Regional clear‐sky shortwave radiative cooling of up to −5 W m −2 due to nitrate is seen, with a much smaller global average cooling of −0.05 W m −2 . Significant enhancements in regional cloud condensation nuclei (at 0.1% supersaturation) and cloud droplet number concentrations are also attributed to nitrate, causing an additional global average shortwave cooling of −0.8 W m −2 . Taking into consideration of changes in both longwave and shortwave radiation under all‐sky conditions, the net change in the top of the atmosphere radiative fluxes induced by including nitrate aerosol is −0.7 W m −2 .

54 ENVIRONMENTAL SCIENCES↗

TH-C-12A-03: Development of Expanded Field Irradiation Technique with Gimbaled X-Ray Head

Purpose: The Vero4DRT has a maximum field size of 150×150 mm{sup 2}. The purposes of this study were to develop an expanded field irradiation technique using a unique gimbaled x-ray head of Vero4DRT and to evaluate its dosimetric characteristic. Methods: The expanded field irradiation consisted of four separate fields with 2.39 degree gimbal rotation around orthogonal two axes. The central beam axis for each field shifted 40 mm from the isocenter for longitudinal and lateral directions, and thus, the field size was expanded up to 230×230 mm{sup 2}. Adjacent region were created at the isocenter (center-adjacent expanded-field) and 20 mm from isocenter (offadjacent expanded-field). To create flat dose distribution in the combined piecewise-fields, the overlapping and gaps regions on the isocenter plane were adjusted with the gimbal rotating and the MLC. To evaluate dosimetric characteristic of the expanded-field, films inserted in water-equivalent phantoms at 50, 100 and 150 mm depth were irradiated and the field size, penumbra, flatness and symmetry were analyzed.In addition, the expandedfield irradiation technique was applied to IMRT. A head and neck IMRT field, which was planned for the conventional linac (Varian Clinac iX), was reproduced with the expanded-field of the Vero4DRT. The simulated dose distribution for the expanded IMRT field was compared to the measured dose distribution. Results: The field size, penumbra, flatness and symmetry of center- and off- adjacent expanded-fields were 230.2–232.1 mm, 7.8–10.7 mm, 2.3–6.5% and –0.5–0.4% at 100 mm depth. The 82.1% area of the expanded IMRT dose distribution was within 5% difference between measurement and simulation, which was analyzed upper 50% dose area, and the 3%/3 mm gamma pass rate was 98.4%. Conclusions: The expandedfield technique was developed using the gimbaled x-ray head. To extend applied targets, such as whole breast irradiations or head and neck IMRT, the expanded-field technique would be effective.

60 APPLIED LIFE SCIENCES↗

The Formation and Evolution of Star Clusters in Interacting Galaxies

Observations of globular clusters show that they have universal lognormal mass functions with a characteristic peak at ∼2×10{sup 5} M{sub ⊙}, but the origin of this peaked distribution is highly debated. Here we investigate the formation and evolution of star clusters (SCs) in interacting galaxies using high-resolution hydrodynamical simulations performed with two different codes in order to mitigate numerical artifacts. We find that massive SCs in the range of ∼10{sup 5.5}--10{sup 7.5} M{sub ⊙} form preferentially in the highly shocked regions produced by galaxy interactions. The nascent cluster-forming clouds have high gas pressures in the range of P/k∼10{sup 8}--10{sup 12} K cm{sup −3}, which is ∼10{sup 4}--10{sup 8} times higher than the typical pressure of the interstellar medium but consistent with recent observations of a pre-super-SC cloud in the Antennae Galaxies. Furthermore, these massive SCs have quasi-lognormal initial mass functions with a peak around ∼10{sup 6} M{sub ⊙}. The number of clusters declines with time due to destructive processes, but the shape and the peak of the mass functions do not change significantly during the course of galaxy collisions. Our results suggest that gas-rich galaxy mergers may provide a favorable environment for the formation of massive SCs such as globular clusters, and that the lognormal mass functions and the unique peak may originate from the extreme high-pressure conditions of the birth clouds and may survive the dynamical evolution.

79 ASTRONOMY AND ASTROPHYSICS↗

Simulation of Vacuum Chamber Pressure Distribution with Surrogate Modeling and Uncertainty Quantification

A major challenge in understanding differences in electric propulsion performance in ground tests and in space operations concerns the pressure distribution within the test vacuum chamber. The chamber backpressure is much higher than experienced in space, modifying thruster performance and plume dynamics. Numerical simulation is a key element to determining the background conditions in non-ideal vacuum chamber environments. An important parameter for the accurate simulation of chamber backpressure is the probability that an atom will stick to a cryogenic panel used to pump away the plume gases. This quantity can be used to model vacuum pumps in particle-based kinetic numerical methods. In this work, a three-dimensional direct simulation Monte Carlo code is used to model neutral xenon atoms flowing from the anode of the H9 Hall Effect Thruster within the University of Michigan’s Large Vacuum Test Facility. Simulated pressures are compared with ion gauge pressure measurements to infer the effective sticking coefficient of the chamber’s vacuum pumps. A pressure predicting surrogate model is developed for inference of pump sticking coefficients and for uncertainty quantification. This information enables accurate and useful kinetic simulations of electric propulsion thruster plasma plumes in vacuum chambers.

DSMC↗

Simulation and analysis of an intermediate frequency (IF) distribution system with applications for Space Station

Simulation and analysis results are described for a wideband fiber optic intermediate frequency distribution channel for a frequency division multiple access (FDMA) system where antenna equipment is remotely located from the signal processing equipment. The fiber optic distribution channel accommodates multiple signals received from a single antenna with differing power levels. Performance parameters addressed are intermodulation degradations, laser noise, and adjacent channel interference, as they impact the overall system design. Simulation results showed that the laser diode modulation level can be allowed to reach 100 percent without considerable degradation. The laser noise must be controlled as to provide a noise floor of less than -90 dBW/Hz. The fiber optic link increases the degradation due to power imbalance yet diminishes the effects of the transmit amplifier nonlinearity. Overall, optimal operation conditions can be found to yield a degradation level of about .1 dB caused by the fiber optic link.

Costello, Thomas A.↗

Statistical analysis of power-size-redshift distributions of extragalactic jets

This paper investigates whether a hot, sparse, yet cosmologically significant intergalactic medium is consistent with data collected from extragalactic radio sources. This is done by use of Monte Carlo simulations which employ previously run pseudohydrodynamical simulations to cover an observational parameter space. These observational parameters include the scale height, central density, and temperature of a (isothermal) galactic halo, and the power of the central engine which drives the jet. The Monte Carlo simulations generate distribution of sizes in bins of (received) power and redshift, which have been compared with observational data using Kolmogorov-Smirnov tests. Results of this analysis are consistent with the existence of an IGM with temperature and density mentioned above. In addition, this analysis suggests that the active lifetime of powerful extragalactic radio sources decreases with increasing power.

Rosen, Alexander↗

Automating the Coupling of ORIGEN with GADRAS via the Fallout Analysis Tool for National Technical Nuclear Forensics

Nuclear forensic teams will be deployed to collect and evaluate fallout samples on the ground in the scenario of a low-yield nuclear detonation in a heavily populated area. Quick non-destructive methods of predicting the quality of the sample before it is analyzed in detail are essential for efficient post-event collections. In this work, the process of exporting Defense Land Fallout Interpretive Code (DELFIC) results into Gamma Detector Response and Analysis Software (GADRAS) has been automated within the Fallout Analysis Tool. This coupling allows for the simulation of detector responses to fallout samples with varying degrees of fractionation. The degree to which the samples are fractionated depends on the location of the samples in the fallout field. In the following study, this phenomenon is examined, as its understanding is important to the investigation of debris distribution. The simulated detector spectra from GADRAS can be used to compare peak ratios of volatile-refractory isotope pairs in order to determine the degree of fractionation. Simulated fractionated fallout samples from DELFIC for a 10 kt, pure 235U fission surface burst were modeled for distances ranging to 256 km out from ground zero, and for times up to 1 week from detonation. The fractionation ratios, also known as r values, from isotope concentrations, photon lines and peak areas of four volatile-refractory pairs were calculated and compared. Fractionation prediction via the peak areas method was evaluated for each pair by comparing the results with the simulated radionuclide inventory.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Evaluating Domestic Hot Water Distribution System Options with Validated Analysis Models

A developing body of work is forming that collects data on domestic hot water consumption, water use behaviors, and energy efficiency of various distribution systems. Transient System Simulation Tool (TRNSYS) is a full distribution system developed that has been validated using field monitoring data and then exercised in a number of climates to understand climate impact on performance. In this study, the Building America team built upon previous analysis modeling work to evaluate differing distribution systems and the sensitivities of water heating energy and water use efficiency to variations of climate, load, distribution type, insulation and compact plumbing practices. Overall, 124 different TRNSYS models were simulated. The results of this work are useful in informing future development of water heating best practices guides as well as more accurate (and simulation time efficient) distribution models for annual whole house simulation programs.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Characterizing the Directionality of Gravitational Wave Emission from Matter Motions within Core-collapse Supernovae

Abstract We analyze the directional dependence of the gravitational wave (GW) emission from 15 3D neutrino radiation hydrodynamic simulations of core-collapse supernovae (CCSNe). Using spin weighted spherical harmonics, we develop a new analytic technique to quantify the evolution of the distribution of GW emission over all angles. We construct a physics-informed toy model that can be used to approximate GW distributions for general ellipsoid-like systems, and use it to provide closed form expressions for the distribution of GWs for different CCSN phases. Using these toy models, we approximate the protoneutron star (PNS) dynamics during multiple CCSN stages and obtain similar GW distributions to simulation outputs. When considering all viewing angles, we apply this new technique to quantify the evolution of preferred directions of GW emission. For nonrotating cases, this dominant viewing angle drifts isotropically throughout the supernova, set by the dynamical timescale of the PNS. For rotating cases, during core bounce and the following tens of milliseconds, the strongest GW signal is observed along the equator. During the accretion phase, comparable—if not stronger—GW amplitudes are generated along the axis of rotation, which can be enhanced by the low T /∣ W ∣ instability. We show two dominant factors influencing the directionality of GW emission are the degree of initial rotation and explosion morphology. Lastly, looking forward, we note the sensitive interplay between GW detector site and supernova orientation, along with its effect on detecting individual polarization modes.

Astronomy & Astrophysics↗

Vertical Distribution of Saharan Dust from the Goddard Three Dimensional Transport Model: Comparison with LITE Data

Spatial distributions of Saharan dust are simulated by the Goddard 3D chemistry transport model driven by meteorological fields from the Goddard Earth Observing System Data Assimilation System, GEOS-2 DAS. An aerosol module accounts for gravitational settling and the aerosol source function is based on surface moisture content and wind speed. The simulated aerosol fields are compared with Total Ozone Mapping Spectrometer (TOMS) by P. Ginoux in this session. The main focus is how ventilation of the assimilation winds will impact the simulated aerosol vertical distribution. Excessive ventilation would transport too much aerosol in the upper troposphere and lead to errors in the long range transport as well as errors in any radiative forcing calculations derived from the simulated aerosol fields. Simulated aerosol fields using GEOS-1 DAS wind fields for September 1994 will be compared with lidar observations during the Lidar In-space Technology Experiment (LITE). The lidar return signal clearly delineates the vertical extent and horizontal variability of the aerosol layer allowing rigorous comparison with the simulated fields. A second focus is the simulation of aerosol scavenging due to moist processes. Three-dimensional precipitation fields can be derived from two possible sources: Assimilated precipitation fields from GEOS-DAS or fields obtained from the Special Sensor Microwave/imager (SSM/1). Of interest is the sensitivity of the modeled aerosol fields to the different precipitation fields driving the moist scavenging module.

Weaver, Clark J.↗

TH-CD-201-01: BEST IN PHYSICS (THERAPY): Assessment of Thin-Film CdTe as An Effective Detector for Microdosimety at Gold-Tissue Interface

Purpose: Presence of interfaces between high and low atomic number materials, often encountered in diagnostic imaging and radiation therapy, leads to radiation dose perturbation. This phenomenon is characterized by a very narrow region of sharp dose enhancement at the interface. The rapid fall-off of the dose enhancement over a very short distance from the interface makes the experimental dosimetry nontrivial. We use an in-house-built inexpensive thin-film Cadmium Telluride (CdTe) photodetector to study this effect at the gold-tissue interface and verify our experimental results with Monte Carlo (MC) modeling. Methods: Three micron thick CdTe photodetectors were fabricated in our lab. One, ten or one hundred micron thick gold foils placed in a tissue-equivalent-phantom were irradiated with a clinical Ir-192 high dose rate source and current measured with a CdTe detector in each case was compared against the current measured for all uniform tissue-equivalent phantom. Percentage signal enhancement (PSE) due to each gold foil was compared against MC modeled percentage dose enhancement (PDE), obtained from the geometry mimicking the experimental setup. Results: The experiment based PSEs due to 1, 10, and 100 micron thick gold foils at the closest measured distance of measurement (12.5 micron) from the interface were 42.6 ± 10.8, 137.0 ± 11.9 and 203.0 ± 15.4 respectively. The corresponding MC modeled PDEs were 38.1 ± 1, 164 ± 1 and 249 ± 1 respectively. The experimental and MC modeled values showed a closer agreement at the larger distances from the interface. Conclusion: The dose enhancement near the gold-tissue interface was measured using an in-house-built high-resolution CdTe-based photodetector and validated with MC simulations. A close agreement of the experimental results with the corresponding MC modeled results shows that CdTe detector can be utilized for mapping interface dose distribution encountered in the application of ionizing radiation.

60 APPLIED LIFE SCIENCES↗

Atomistic simulations of nanowelding of single-crystal and amorphous gold nanowires

The mechanism and quality of the welding of single-crystal (SC) and amorphous gold nanowires (NWs) with head-to-head contact are studied using molecular dynamics simulations based on the second-moment approximation of the many-body tight-binding potential. The results are discussed in terms of atomic trajectories, slip vectors, stress, and radial distribution function. Simulation results show that the alignment for the amorphous NWs during welding is easier than that for the SC NWs due to the former's relatively stable geometry. A few dislocations nucleate and propagate on the (111) close-packed plane (slip plane) inside the SC NWs during the welding and stretching processes. During welding, an incomplete jointing area first forms through the interactions of the van der Waals attractive force, and the jointing area increases with increasing extent of contact between the two NWs. A crystallization transition region forms in the jointing area for the welding of SC-amorphous or amorphous-SC NWs. With increasing interference, an amorphous gold NW shortens more than does a SC gold NW due to the former's relatively poor strength. The pressure required for welding decreases with increasing temperature.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗