Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “consequence”

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,531 records · Page 85
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Prediction of Metabolite Concentrations, Rate Constants and Post-Translational Regulation Using Maximum Entropy-Based Simulations with Application to Central Metabolism of Neurospora crassa

We report the application of a recently proposed approach for modeling biological systems using a maximum entropy production rate principle in lieu of having in vivo rate constants. The method is applied in four steps: (1) a new ordinary differential equation (ODE) based optimization approach based on Marcelin’s 1910 mass action equation is used to obtain the maximum entropy distribution; (2) the predicted metabolite concentrations are compared to those generally expected from experiments using a loss function from which post-translational regulation of enzymes is inferred; (3) the system is re-optimized with the inferred regulation from which rate constants are determined from the metabolite concentrations and reaction fluxes; and finally (4) a full ODE-based, mass action simulation with rate parameters and allosteric regulation is obtained. From the last step, the power characteristics and resistance of each reaction can be determined. The method is applied to the central metabolism of Neurospora crassa and the flow of material through the three competing pathways of upper glycolysis, the non-oxidative pentose phosphate pathway, and the oxidative pentose phosphate pathway are evaluated as a function of the NADP/NADPH ratio. It is predicted that regulation of phosphofructokinase (PFK) and flow through the pentose phosphate pathway are essential for preventing an extreme level of fructose 1,6-bisphophate accumulation. Such an extreme level of fructose 1,6-bisphophate would otherwise result in a glassy cytoplasm with limited diffusion, dramatically decreasing the entropy and energy production rate and, consequently, biological competitiveness.

Cannon, William R.↗
DOE OSTI2020

Connecting Risk and Resilience for a Power System Using the Portland Hills Fault Case Study

Active seismic faults in the Pacific Northwest area have encouraged electric utilities in the region to deeply contemplate and proactively intervene to support grid resilience. To further this effort this research introduces Monte Carlo (MC)-based power system modeling as a means to inform the Performance Based Earthquake Engineering method and simulates 100,000 sample earthquakes of a 6.8 magnitude (M6.8) Portland Hills Fault (PHF) scenario in the Portland General Electric (PGE) service territory as a proof of concept. This paper also proposes the resilience metric Seismic Load Recovery Factor (SLRF) to quantify the recovery of a downed power system and thus can be used to quantify earthquake economic risk. Using MC results, the SLRF was evaluated to be 19.7 h and the expected economic consequence cost of a M6.8 PHF event was found to be $\$180$ million with an annualized risk of $\$90,000$ given the event’s 1 in 2000 year probability of occurrence. The MC results also identified the eight most consequential substations in the PGE system—i.e., those that contributed to maximum load loss. This paper concludes that retrofitting these substations reduced the expected consequence cost of a M6.8 PHF event to $117 million.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2018

Implications of Emerging Vehicle Technologies on Rare Earth Supply and Demand in the United States

In this article, we explore the long-term demand and supply potentials of rare earth elements in alternative energy vehicles (AEVs) in the United States until 2050. Using a stock-flow model, we compare a baseline scenario with scenarios that incorporate an exemplary technological innovation: a novel aluminum–cerium–magnesium alloy. We find that the introduction of the novel alloy demonstrates that even low penetration rates can exceed domestic cerium production capacity, illustrating possible consequences of technological innovations to material supply and demand. End-of-life vehicles can, however, overtake domestic mining as a source of materials, calling for proper technologies and policies to utilize this emerging source. The long-term importing of critical materials in manufactured and semi-manufactured products shifts the location of material stocks and hence future secondary supply of high-value materials, culminating in a double benefit to the importing country. This modeling approach is adaptable to the study of varied scenarios and materials, linking technologies with supply and demand dynamics in order to understand their potential economic and environmental consequences.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2025

The Principle of Maximum Conformality Correctly Resolves the Renormalization-Scheme-Dependence Problem

In this paper, we clarify a serious misinterpretation and consequent misuse of the Principle of Maximum Conformality (PMC), which also can serve as a mini-review of PMC. In a recently published article, P. M. Stevenson has claimed that “the PMC is ineffective and does nothing to resolve the renormalization-scheme-dependence problem”, concluding incorrectly that the success of PMC predictions is due to the PMC being a “laborious, ad hoc, and back-door” version of the Principle of Minimal Sensitivity (PMS). We show that such conclusions are incorrect, deriving from a misinterpretation of the PMC and an overestimation of the applicability of the PMS. The purpose of the PMC is to achieve precise fixed-order pQCD predictions, free from conventional renormalization schemes and scale ambiguities. We demonstrate that the PMC predictions satisfy all the self-consistency conditions of the renormalization group and standard renormalization-group invariance; the PMC predictions are thus independent of any initial choice of renormalization scheme and scale. The scheme independence of the PMC is also ensured by commensurate scale relations, which relate different observables to each other. Moreover, in the Abelian limit, the PMC dovetails into the well-known Gell-Mann–Low framework, a method universally revered for its precision in QED calculations. Due to the elimination of factorially divergent renormalon terms, the PMC series not only attains a convergence behavior far superior to that of its conventional counterparts but also deftly curtails any residual scale dependence caused by the unknown higher-order terms. This refined convergence, coupled with its robust suppression of residual uncertainties, furnishes a sound and reliable foundation for estimating the contributions from unknown higher-order terms. Anchored in the bedrock of standard renormalization-group invariance, the PMC simultaneously eradicates the factorial divergences and eliminates superfluous systematic errors, which inversely provides a good foundation for achieving high-precision pQCD predictions. Consequently, owing to its rigorous theoretical underpinnings, the PMC is eminently applicable to virtually all high-energy hadronic processes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

On the role of the ΩΓ limit in the formation of population III massive stars

We explore the role of the modified Eddington limit due to rapid rotation (the so-called ΩΓ limit) in the formation of Population III stars. We performed one-dimensional stellar evolution simulations of zero-metallicity protostars accreting mass at a very high rate ( M-dot ∼10{sup −3} M{sub ⊙} yr{sup −1}) and dealt with stellar rotation as a separate post-process. The protostar would reach the Keplerian rotation very soon after the onset of mass accretion, but mass accretion would continue as stellar angular momentum is transferred outward to the accretion disk by viscous stress. The envelope of the protostar expands rapidly when the stellar mass reaches ∼5–7 M{sub ⊙} and the Eddington factor increases sharply. This makes the protostar rotate critically at a rate that is significantly below the Keplerian value (i.e., the ΩΓ limit). The resultant positive gradient of the angular velocity in the boundary layer between the protostar and the Keplerian disk prohibits angular momentum transport from the star to the disk, and consequently further rapid mass accretion. This would prevent the protostar from growing significantly beyond 20–40 M{sub ⊙}. Another important consequence of the ΩΓ limit is that the protostar can remain fairly compact (R ≲ 50 R{sub ⊙}) and avoid a fluffy structure (R ≳ 500 R{sub ⊙}) that is usually found with a very high rate of mass accretion. This effect would make the protostar less prone to binary interactions during the protostar phase. Although our analysis is based on Population III protostar models, this role of the ΩΓ limit would be universal in the formation process of massive stars, regardless of metallicity.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

KINETIC STUDY OF RADIATION-REACTION-LIMITED PARTICLE ACCELERATION DURING THE RELAXATION OF UNSTABLE FORCE-FREE EQUILIBRIA

Many powerful and variable gamma-ray sources, including pulsar wind nebulae, active galactic nuclei and gamma-ray bursts, seem capable of accelerating particles to gamma-ray emitting energies efficiently over very short timescales. These are likely due to the rapid dissipation of electromagnetic energy in a highly magnetized, relativistic plasma. In order to understand the generic features of such processes, we have investigated simple models based on the relaxation of unstable force-free magnetostatic equilibria. In this work, we make the connection between the corresponding plasma dynamics and the expected radiation signal, using 2D particle-in-cell simulations that self-consistently include synchrotron radiation reactions. We focus on the lowest order unstable force-free equilibrium in a 2D periodic box. We find that rapid variability, with modest apparent radiation efficiency as perceived by a fixed observer, can be produced during the evolution of the instability. The “flares” are accompanied by an increased polarization degree in the high energy band, with rapid variation in the polarization angle. Furthermore, the separation between the acceleration sites and the synchrotron radiation sites for the highest energy particles facilitates acceleration beyond the synchrotron radiation reaction limit. We also discuss the dynamical consequences of the radiation reaction, and some astrophysical applications of this model. Our current simulations with numerically tractable parameters are not yet able to reproduce the most dramatic gamma-ray flares, e.g., from the Crab Nebula. Higher magnetization studies are promising and will be carried out in the future.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

SHOULD ONE USE THE RAY-BY-RAY APPROXIMATION IN CORE-COLLAPSE SUPERNOVA SIMULATIONS?

We perform the first self-consistent, time-dependent, multi-group calculations in two dimensions (2D) to address the consequences of using the ray-by-ray+ transport simplification in core-collapse supernova simulations. Such a dimensional reduction is employed by many researchers to facilitate their resource-intensive calculations. Our new code (Fornax) implements multi-D transport, and can, by zeroing out transverse flux terms, emulate the ray-by-ray+ scheme. Using the same microphysics, initial models, resolution, and code, we compare the results of simulating 12, 15, 20, and 25 M {sub ⊙} progenitor models using these two transport methods. Our findings call into question the wisdom of the pervasive use of the ray-by-ray+ approach. Employing it leads to maximum post-bounce/pre-explosion shock radii that are almost universally larger by tens of kilometers than those derived using the more accurate scheme, typically leaving the post-bounce matter less bound and artificially more “explodable.” In fact, for our 25 M {sub ⊙} progenitor, the ray-by-ray+ model explodes, while the corresponding multi-D transport model does not. Therefore, in two dimensions, the combination of ray-by-ray+ with the axial sloshing hydrodynamics that is a feature of 2D supernova dynamics can result in quantitatively, and perhaps qualitatively, incorrect results.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

BABCOCK–LEIGHTON SOLAR DYNAMO: THE ROLE OF DOWNWARD PUMPING AND THE EQUATORWARD PROPAGATION OF ACTIVITY

The key elements of the Babcock–Leighton dynamos are the generation of poloidal field through decay and the dispersal of tilted bipolar active regions and the generation of toroidal field through the observed differential rotation. These models are traditionally known as flux transport dynamo models as the equatorward propagations of the butterfly wings in these models are produced due to an equatorward flow at the bottom of the convection zone. Here we investigate the role of downward magnetic pumping near the surface using a kinematic Babcock–Leighton model. We find that the pumping causes the poloidal field to become predominately radial in the near-surface shear layer, which allows the negative radial shear to effectively act on the radial field to produce a toroidal field. We observe a clear equatorward migration of the toroidal field at low latitudes as a consequence of the dynamo wave even when there is no meridional flow in the deep convection zone. Both the dynamo wave and the flux transport type solutions are thus able to reproduce some of the observed features of the solar cycle including the 11-year periodicity. The main difference between the two types of solutions is the strength of the Babcock–Leighton source required to produce the dynamo action. A second consequence of the magnetic pumping is that it suppresses the diffusion of fields through the surface, which helps to allow an 11-year cycle at (moderately) larger values of magnetic diffusivity than have previously been used.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Observation and characterization of a cosmic muon neutrino flux from the northern hemisphere using six years of IceCube data

The IceCube Collaboration has previously discovered a high-energy astrophysical neutrino flux using neutrino events with interaction vertices contained within the instrumented volume of the IceCube detector. We present a complementary measurement using charged current muon neutrino events where the interaction vertex can be outside this volume. As a consequence of the large muon range the effective area is significantly larger but the field of view is restricted to the Northern Hemisphere. IceCube data from 2009 through 2015 have been analyzed using a likelihood approach based on the reconstructed muon energy and zenith angle. At the highest neutrino energies between 194 TeV and 7.8 PeV a significant astrophysical contribution is observed, excluding a purely atmospheric origin of these events at 5.6σ significance. The data are well described by an isotropic, unbroken power-law flux with a normalization at 100 TeV neutrino energy of (0.90{sub −0.27}{sup +0.30})×10{sup −18} GeV{sup −1} cm{sup −2} s{sup −1} sr{sup −1} and a hard spectral index of γ=2.13±0.13. The observed spectrum is harder in comparison to previous IceCube analyses with lower energy thresholds which may indicate a break in the astrophysical neutrino spectrum of unknown origin. The highest-energy event observed has a reconstructed muon energy of (4.5±1.2) PeV which implies a probability of less than 0.005% for this event to be of atmospheric origin. Analyzing the arrival directions of all events with reconstructed muon energies above 200 TeV no correlation with known γ-ray sources was found. Using the high statistics of atmospheric neutrinos we report the current best constraints on a prompt atmospheric muon neutrino flux originating from charmed meson decays which is below 1.06 in units of the flux normalization of the model in Enberg et al.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

The Gemini Planet Imager View of the HD 32297 Debris Disk

We present new H-band scattered light images of the HD 32297 edge-on debris disk obtained with the Gemini Planet Imager. The disk is detected in total and polarized intensity down to a projected angular separation of 0.″15, or 20 au. On the other hand, the large-scale swept-back halo remains undetected, likely a consequence of its markedly blue color relative to the parent body belt. We analyze the curvature of the disk spine and estimate a radius of ≈100 au for the parent body belt, smaller than past scattered light studies but consistent with thermal emission maps of the system. We employ three different flux-preserving post-processing methods to suppress the residual starlight and evaluate the surface brightness and polarization profile along the disk spine. Unlike past studies of the system, our high-fidelity images reveal the disk to be highly symmetric and devoid of morphological and surface brightness perturbations. We find the dust scattering properties of the system to be consistent with those observed in other debris disks, with the exception of HR 4796. Finally, we find no direct evidence for the presence of a planetary-mass object in the system.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

THE ROLE OF KELVIN–HELMHOLTZ INSTABILITY FOR PRODUCING LOOP-TOP HARD X-RAY SOURCES IN SOLAR FLARES

We propose a model for the formation of loop-top hard X-ray (HXR) sources in solar flares through the inverse Compton mechanism, scattering the surrounding soft X-ray (SXR) photons to higher energy HXR photons. We simulate the consequences of a flare-driven energy deposit in the upper chromosphere in the impulsive phase of single loop flares. The consequent chromosphere evaporation flows from both footpoints reach speeds up to hundreds of kilometers per second, and we demonstrate how this triggers Kelvin–Helmholtz instability (KHI) in the loop top, under mildly asymmetric conditions, or more toward the loop flank for strongly asymmetric cases. The KHI vortices further fragment the magnetic topology into multiple magnetic islands and current sheets, and the hot plasma within leads to a bright loop-top SXR source region. We argue that the magnetohydrodynamic turbulence that appears at the loop apex could be an efficient accelerator of non-thermal particles, which the island structures can trap at the loop-top. These accelerated non-thermal particles can upscatter the surrounding thermal SXR photons emitted by the extremely hot evaporated plasma to HXR photons.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Spectral evolution in high redshift quasars from the final Baryon Oscillation Spectroscopic Survey sample

We report on the diversity in quasar spectra from the Baryon Oscillation Spectroscopic Survey. After filtering the spectra to mitigate selection effects and Malmquist bias associated with a nearly flux-limited sample, we create high signal-to-noise ratio composite spectra from 58,656 quasars (2.1⩽z⩽3.5), binned by luminosity, spectral index, and redshift. With these composite spectra, we confirm the traditional Baldwin effect (BE, i.e., the anti-correlation of C iv equivalent width (EW) and luminosity) that follows the relation W{sub λ}∝L{sup β{sub w}} with slope β{sub w}=−0.35±0.004, −0.35 ± 0.005, and −0.41 ± 0.005 for z = 2.25, 2.46, and 2.84, respectively. In addition to the redshift evolution in the slope of the BE, we find redshift evolution in average quasar spectral features at fixed luminosity. The spectroscopic signature of the redshift evolution is correlated at 98% with the signature of varying luminosity, indicating that they arise from the same physical mechanism. At a fixed luminosity, the average C iv FWHM decreases with increasing redshift and is anti-correlated with C iv EW. The spectroscopic signature associated with C iv FWHM suggests that the trends in luminosity and redshift are likely caused by a superposition of effects that are related to black hole mass and Eddington ratio. The redshift evolution is the consequence of a changing balance between these two quantities as quasars evolve toward a population with lower typical accretion rates at a given black hole mass.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Early GRB afterglows from reverse shocks in ultra-relativistic, long-lasting winds

We develop a model of early Gamma-Ray Burst (GRB) afterglows with dominant X-ray contribution from the reverse shock (RS) propagating in highly relativistic (Lorentz factor γ{sub w} ∼ 10{sup 6}) magnetized wind of a long-lasting central engine. The model reproduces, in a fairly natural way, the overall trends and yet allows for variations in the temporal and spectral evolution of early optical and X-ray afterglows. The high energy and the optical synchrotron emission from the RS particles occurs in the fast cooling regime; the resulting synchrotron power L{sub s} is a large fraction of the wind luminosity, L{sub s}≈L{sub w}/√(1+σ{sub w}) (L{sub w} and σ{sub w} are wind power and magnetization). Thus, plateaus—parts of afterglow light curves that show slowly decreasing spectral power—are a natural consequence of the RS emission. Contribution from the forward shock (FS) is negligible in the X-rays, but in the optical both FS and RS contribute similarly: FS optical emission is in the slow cooling regime, producing smooth components, while RS optical emission is in the fast cooling regime, and thus can both produce optical plateaus and account for fast optical variability correlated with the X-rays, e.g., due to changes in the wind properties. We discuss how the RS emission in the X-rays and combined FS and RS emission in the optical can explain many puzzling properties of early GRB afterglows.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

On the Detectability of the X-tilde {sup 2}A{sup ′′} HSS, HSO, and HOS Radicals in the Interstellar Medium

X-tilde {sup 2}A{sup ′′} HSS has yet to be observed in the gas phase in the interstellar medium (ISM). HSS has been observed in cometary material and in high abundance. However, its agglomeration to such bodies or dispersal from them has not been observed. Similarly, HSO and HOS have not been observed in the ISM, either, even though models support their formation from reactions of known sulfur monoxide and hydrogen molecules, among other pathways. Consequently, this work provides high-level, quantum chemical rovibrational spectroscopic constants and vibrational frequencies in order to assist in interstellar searches for these radical molecules. Furthermore, the HSO−HOS isomerization energy is determined to be 3.63 kcal mol{sup −1}, in line with previous work, and the dipole moment of HOS is 36% larger at 3.87 D than HSO, making the less stable isomer more rotationally intense. Finally, the S−S bond strength in HSS is shown to be relatively weak at 30% of the typical disulfide bond energy. Consequently, HSS may degrade into SH and sulfur atoms, making any ISM abundance of HSS likely fairly low, as recent interstellar surveys have observed.

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

The Formation of Magnetic Depletions and Flux Annihilation Due to Reconnection in the Heliosheath

The misalignment of the solar rotation axis and the magnetic axis of the Sun produces a periodic reversal of the Parker spiral magnetic field and the sectored solar wind. The compression of the sectors is expected to lead to reconnection in the heliosheath (HS). We present particle-in-cell simulations of the sectored HS that reflect the plasma environment along the Voyager 1 and 2 trajectories, specifically including unequal positive and negative azimuthal magnetic flux as seen in the Voyager data. Reconnection proceeds on individual current sheets until islands on adjacent current layers merge. At late time, bands of the dominant flux survive, separated by bands of deep magnetic field depletion. The ambient plasma pressure supports the strong magnetic pressure variation so that pressure is anticorrelated with magnetic field strength. There is little variation in the magnetic field direction across the boundaries of the magnetic depressions. At irregular intervals within the magnetic depressions are long-lived pairs of magnetic islands where the magnetic field direction reverses so that spacecraft data would reveal sharp magnetic field depressions with only occasional crossings with jumps in magnetic field direction. This is typical of the magnetic field data from the Voyager spacecraft. Voyager 2 data reveal that fluctuations in the density and magnetic field strength are anticorrelated in the sector zone, as expected from reconnection, but not in unipolar regions. The consequence of the annihilation of subdominant flux is a sharp reduction in the number of sectors and a loss in magnetic flux, as documented from the Voyager 1 magnetic field and flow data.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

A Panchromatic View of Brown Dwarf Aurorae

Stellar coronal activity has been shown to persist into the low-mass star regime, down to late M-dwarf spectral types. However, there is now an accumulation of evidence suggesting that at the end of the main sequence, there is a transition in the nature of the magnetic activity from chromospheric and coronal to planet-like and auroral, from local impulsive heating via flares and MHD wave dissipation to energy dissipation from strong large-scale magnetospheric current systems. We examine this transition and the prevalence of auroral activity in brown dwarfs through a compilation of multiwavelength surveys of magnetic activity, including radio, X-ray, and optical. We compile the results of those surveys and place their conclusions in the context of auroral emission as a consequence of large-scale magnetospheric current systems that accelerate energetic electron beams and drive the particles to impact the cool atmospheric gas. We explore the different manifestations of auroral phenomena, like H α , in brown dwarf atmospheres and define their distinguishing characteristics. We conclude that large-amplitude photometric variability in the near-infrared is most likely a consequence of clouds in brown dwarf atmospheres, but that auroral activity may be responsible for long-lived stable surface features. We report a connection between auroral H α emission and quiescent radio emission in electron cyclotron maser instability pulsing brown dwarfs, suggesting a potential underlying physical connection between quiescent and auroral emissions. We also discuss the electrodynamic engines powering brown dwarf aurorae and the possible role of satellites around these systems both to power the aurorae and seed the magnetosphere with plasma.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Apex Dips of Experimental Flux Ropes: Helix or Cusp?

We present a new theory for the presence of apex dips in certain experimental flux ropes. Previously such dips were thought to be projections of a helical loop axis generated by the kink instability. However, new evidence from experiments and simulations suggest that the feature is a 2D cusp rather than a 3D helix. The proposed mechanism for cusp formation is a density pileup region generated by nonlinear interaction of neutral gas cones emitted from fast-gas nozzles. The results indicate that density perturbations can result in large distortions of an erupting flux rope, even in the absence of significant pressure or gravitational forces. The density pileup at the apex also suppresses the m = 1 kink mode by acting as a stationary node. Consequently, more accurate density profiles should be considered when attempting to model the stability and shape of solar and astrophysical flux ropes.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Sensor Distortion Effects in Photon Monte Carlo Simulations

We present a detailed method to simulating sensor distortions using a photon and electron Monte Carlo method. We use three-dimensional electrostatic simulations to parameterize the perturbed electric field profile for nonideal sensor details. We follow the conversion of simulated photons, and the subsequent response of the converted electrons to the electric field pattern. These nonideal sensor details can be implemented efficiently in a Monte Carlo approach. We demonstrate that the nonideal sensor distortions have a variety of observable consequence including the modification of the astrometric pattern; the distortion of the electron diffusion size and shape; and the distortion of flats. We show analytic validation of the diffusion physics, reproduce two kinds of edge distortion, and show qualitative validation of field-free regions, lithography errors, and fringing. We also demonstrate that there are two related effects of doping variation having different observable consequences. We show that field distortions from accumulated electrons lead to intensity-dependent point-spread functions and the sublinear variance in flats. The method is implemented in the Photon Simulator and the code is publicly available.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗