Energy Mapping Method Research Articles

Simulation of Explosion Crater of Earth Penetrating Nuclear Bomb Based on Equivalency to TNT Mass


 
 
 The structure size of the crater formed by the earth-penetrating nuclear bomb explosion is one of the important parameters for evaluating the earth-penetrating nuclear bomb damage power. Obtaining the structure size of the crater formed by the earth-penetrating nuclear bomb explosion with different yields is great significance for the evaluation and design of the nuclear bomb damage power. In this study, considering the contradictory relationship between the structure size of the earth-penetrating nuclear bomb and the structure size of the equivalent TNT mass, we propose to use the equivalent energy mapping method to realize the finite element numerical simulation of the earth-penetrating nuclear bomb exploding into craters analyzed and compared the simulation results with the structure size of the crater formed by the ESS nuclear bomb explosion in the United States. The analysis results show that the error between the simulated crater radius and the real crater radius is 3.26%, and the error between the simulated crater depth and the real crater depth is 28.57 %. It meets the calculation accuracy error range of crater formation from nuclear explosion to chemical explosion. Therefore, this method provides an effective numerical simulation method and a means of obtaining the structural size data of the explosion crater for the earth-penetrating nuclear bomb cratering.
 
 

Development of a multi-layer quality assurance program to evaluate the uncertainty of deformable dose accumulation in adaptive radiotherapy.

Deformable dose accumulation (DDA) has uncertainties which impede the implementation of DDA-based adaptive radiotherapy (ART) in clinic. The purpose of this study is to develop a multi-layer quality assurance (MLQA) program to evaluate uncertainties in DDA. A computer program is developed to generate a pseudo-inverse displacement vector field (DVF) for each deformable image registration (DIR) performed in Accuray's PreciseART. The pseudo-inverse DVF is first used to calculate a pseudo-inverse consistency error (PICE) and then implemented in an energy and mass congruent mapping (EMCM) method to reconstruct a deformed dose. The PICE is taken as a metric to estimate DIR uncertainties. A pseudo-inverse dose agreement rate (PIDAR) is used to evaluate the consequence of the DIR uncertainties in DDA and the principle of energy conservation is used to validate the integrity of dose mappings. The developed MLQA program was tested using the data collected from five representative cancer patients treated with tomotherapy. DIRs were performed in PreciseART to generate primary DVFs for the five patients. The fidelity index and PICE of these DVFs on average are equal to 0.028mm and 0.169mm, respectively. With the criteria of 3mm/3% and 5mm/5%, the PIDARs of the PreciseART-reconstructed doses are 73.9±4.4% and 87.2±3.3%, respectively. The PreciseART and EMCM-based dose reconstructions have their deposited energy changed by 5.6±3.9% and 2.6±1.5% in five GTVs, and by 9.2±7.8% and 4.7±3.6% in 30 OARs, respectively. A pseudo-inverse map-based EMCM program has been developed to evaluate DIR and dose mapping uncertainties. This program could also be used as a sanity check tool for DDA-based ART.

Heat Mapping, a Method for Enhancing the Sustainability of the Smart District Heat Networks

District heating (DH) has a major potential to increase the efficiency, security, and sustainability of energy management at the community scale. However, there is a huge challenge for decision makers due to the lack of knowledge about thermal energy demand during a year. Thermal energy demand is strongly dependent on the outdoor temperature, building area, and activities. In this context, this paper presents an innovative monthly thermal energy mapping method to calculate and visualize heat demand accurately for various types of buildings. The method includes three consecutive phases: (i) calculating energy loss, (ii) completing a dataset that includes energy and building information, and (iii) generating the monthly heat demand maps for the community. Determining the amount of demand and the best location for energy generators from the perspective of energy efficiency in a DH system in an urban context is one of the important applications of heat maps. Exploring heat demand characteristics and visualizing them on maps is the foundation of smart DHs.

Magnetic interactions in BiFeO3 : A first-principles study

First-principles calculations, in combination with the four-state energy mapping method, are performed to extract the magnetic interaction parameters of multiferroic BiFeO$_3$. Such parameters include the symmetric exchange (SE) couplings and the Dzyaloshinskii-Moriya (DM) interactions up to second nearest neighbors, as well as the single ion anisotropy (SIA). All magnetic parameters are obtained not only for the $R3c$ structural ground state, but also for the $R3m$ and $R\bar{3}c$ phases in order to determine the effects of ferroelectricity and antiferrodistortion distortions, respectively, on these magnetic parameters. In particular, two different second-nearest neighbor couplings are identified and their origins are discussed in details. Moreover, Monte-Carlo (MC) simulations using a magnetic Hamiltonian incorporating these first-principles-derived interaction parameters are further performed. They result (i) not only in the accurate prediction of the spin-canted G-type antiferromagnetic structure and of the known magnetic cycloid propagating along a $<$1$\bar{1}$0$>$ direction, as well as their unusual characteristics (such as a weak magnetization and spin-density-waves, respectively); (ii) but also in the finding of another cycloidal state of low-energy and that awaits to be experimentally confirmed. Turning on and off the different magnetic interaction parameters in the MC simulations also reveal the precise role of each of them on magnetism.

Electronic structure and surface properties ofMgB2(0001) upon oxygen adsorption

We use density-functional theory to investigate the bulk and surface properties of ${\mathrm{MgB}}_{2}$. The unique bonding structure of ${\mathrm{MgB}}_{2}$ is investigated by Bader's atoms-in-molecules, charge density difference, and occupancy projected band structure analyses. Oxygen adsorption on the charge-depleted surfaces of ${\mathrm{MgB}}_{2}$ is studied by a surface potential energy mapping method, reporting a complete map including low-symmetry binding sites. The B-terminated ${\mathrm{MgB}}_{2}$(0001) demonstrates reconstruction of the graphenelike B layer, and the reconstructed geometry exposes a threefold site of the subsurface Mg, making it accessible from the surface. Detailed reconstruction mechanisms are studied by simulated annealing method based on ab initio molecular dynamics and nudged elastic band calculations. The surface clustering of B atoms significantly modifies the B $2p$ states to occupy low energy valence states. The present paper emphasizes that a thorough understanding of the surface phase may explain an apparent inconsistency in the experimental surface characterization of ${\mathrm{MgB}}_{2}$. Furthermore, these results suggest that the surface passivation can be an important technical challenge when it comes to development of a superconducting device using ${\mathrm{MgB}}_{2}$.

Aircraft integrated structural health monitoring using lasers, piezoelectricity, and fiber optics

Various structural health monitoring (SHM) systems that have been developed based on laser ultrasonics and fiber optics are introduced in this paper. The systems are used to realize the new SHM paradigm for ground SHM, called the Smart Hangar. Guided wave ultrasonic propagation imaging (G-UPI) technology is implemented in the Smart Hangar in the form of built-in and mobile G-UPI systems. The laser-induced guided waves generated by pulsed beam scanning in the wings and fuselage of an aircraft are captured by a fiber optic, piezoelectric, or laser ultrasonic sensor, and their propagation is visualized. The wave propagation video is further processed to visualize damage in the presence of multiple sources using a multi-time-frame ultrasonic energy mapping (mUEM) method. For in-flight SHM, laser scanning-based smartification of the structure with sensors is presented, and an event localization method based on fiber optics and piezoelectric sensing is also introduced. Optic sensors are also utilized to reconstruct the wing deformation from the measured strain. The wing deformation and impact localization information is transferred to a ground pilot in the case of unmanned aerial vehicles (UAV), and the ground pilot can feel the wing deformation and impact by using a pilot arm-wearable haptic interface, which makes it possible to achieve human-UAV interactive decision making.

FPGA-based ultrasonic energy mapping with source removal method for damage visualization in composite structures

This paper presents the implementation of an ultrasonic energy mapping (UEM) method in a field programmable gate array (FPGA)-based processing board. This method is integrated into an ultrasonic propagation imaging (UPI) system for real-field application. UEM method involves two-dimensional Fourier transform and wavenumber filtering to capture the energy scattering waves. Performing multi-dimensional Fourier transform in FPGA is a challenging task because it cannot fit into a single FPGA chip, and therefore, external memory is required to handle a large amount of data. The row–column decomposition method is used to perform 2D Fourier transform using two DDR3 SDRAM memories. Additionally, a source removal function is added into the system to highlight the damage when high energy at the source location occurs. The UPI system was applied for inspection of a composite wing box with impact damages. The UEM results showed that a damage of 20 × 20 mm2 was detected successfully in the presence of the source inside the scan area. Moreover, the UEM result is produced in less than 5 s for 57,600 points, which is practical for a real-world, non-destructive testing application.

Freeform illumination lens design combining energy and intensity mapping

A composite method combining energy and intensity mapping is proposed to address the issue of surface error caused by the irregular sampling phenomenon in freeform illumination lens design. In the combined method, the central region of the freeform lens is designed by the intensity mapping method, whereas the peripheral region is designed by the energy-mapping method. Furthermore, an iterative feedback optimization is added to scale out the application in extended light sources. As an evaluating example, a freeform lens with a 120-deg viewing angle, fitted with an appropriate number of points is designed by the proposed combined method. Compared with that designed by the energy method, the lens forms a more uniform illumination on the target surface without the appearance of a hot spot in the central region. The proposed method also exhibits superiority in extended-source design, where only a one-time optimization is needed to achieve the preset uniformity with the proper choice of coefficients.

Giant biquadratic interaction-induced magnetic anisotropy in the iron-based superconductorAxFe2−ySe2

The emergence of the electron-pocket only iron-based superconductor AxFe2-ySe2 (A = alkali metal) challenges the Fermi-surface nesting picture established in iron-pnictides. It was widely believed that magnetism is correlated with the superconductivity in AxFe2-ySe2. Unfortunately, the highly anisotropic exchange parameters and the disagreement between theoretical calculations and experimental results triggered a fierce debate on the nature of magnetism in AxFe2-ySe2. Here we find that the strong magnetic anisotropy is from the anisotropic biquadratic interaction. In order to accurately obtain the magnetic interaction parameters, we propose a universal method, which does not need including other high energy configurations as did in conventional energy mapping method. We show that our model successfully captures the magnetic interactions in AxFe2-ySe2 and correctly predicts the spin wave spectrum, in quantitative agreement with the experimental observation. These results suggest that the local moment picture, including the biquadratic term, can describe accurately the magnetic properties and spin excitations in AxFe2-ySe2, which sheds new light on the future study of the high-Tc iron-based superconductors.

Freeform lens design for light-emitting diode uniform illumination by using a method of source-target luminous intensity mapping.

We present a freeform lens for application to LED uniform illumination. This lens, which is designed with a method of simple source-target luminous intensity mapping, can produce irradiance uniformity of greater than 0.8 and optical efficiency above 90% with an arbitrary half-beam angle greater than 45deg. Typically, as compared with a conventional source-target energy mapping method, this design method can achieve better optical performance of lenses for general LED lighting. When a non-Lambertian-type light source is employed, for example, the chip on board LED, the use of the method can result in a compact LED lens without losing the optical performances of high irradiance uniformity and high optical efficiency as yielded by lenses for Lambertian-type LED light sources.

Design a light pattern of multiple concentric circles for LED fishing lamps using Fourier series and an energy mapping method.

Fourier series and an energy mapping method were used in this study to design a lens that produces a light pattern of multiple concentric circles (LPMCC) for a light-emitting diode (LED) fishing lamp. Fourier series were used to represent the light intensity distribution curve (LIDC) of the LPMCC light pattern. Energy mapping involves performing angular energy mapping based on the LIDCs of an LED light source and LPMCC to design a freeform lens. Type I and Type II LPMCC lenses were designed according to the phototaxis behavior of fish to create a LPMCC light pattern of interleaving light-dark zones that attracts fish shoals to stay in an area for a long period. The experimental results indicated that, in comparing the LIDCs of the Type I and II lenses with the respective simulation values, the normalized cross-correlation (NCC) value reached 96%. According to a 24-hour observation of the phototaxis of Poecilia reticulata to evaluate the effectiveness of the proposed light pattern to attract fish, when a fish shoal was habituated to a light source that emitted constant illumination light, it gradually moved away from the intense light zone and hovered around the junction of the light and dark zones. In the future, the design used in this study can be applied to LED fishing lamps to replace traditional fishing lamps.

Investigation of voxel warping and energy mapping approaches for fast 4D Monte Carlo dose calculations in deformed geometries using VMC++

A new deformable geometry class for the VMC++ Monte Carlo code was implemented based on the voxel warping method. Alternative geometries which use tetrahedral sub-elements were implemented and efficiency improvements investigated. A new energy mapping method, based on calculating the volume overlap between deformed reference dose grid and the target dose grid, was also developed. Dose calculations using both the voxel warping and energy mapping methods were compared in simple phantoms as well as a patient geometry. The new deformed geometry implementation in VMC++ increased calculation times by approximately a factor of 6 compared to standard VMC++ calculations in rectilinear geometries. However, the tetrahedron-based geometries were found to improve computational efficiency, relative to the dodecahedron-based geometry, by a factor of 2. When an exact transformation between the reference and target geometries was provided, the voxel and energy warping methods produced identical results. However, when the transformation is not exact, there were discrepancies in the energy deposited on the target geometry which lead to significant differences in the dose calculated by the two methods. Preliminary investigations indicate that these energy differences may correlate with registration errors; however, further work is needed to determine the usefulness of this metric for quantifying registration accuracy.

CT-Based Attenuation Correction on the FLEX Triumph Preclinical PET/CT Scanner

Positron Emission Tomography (PET) has emerged as a valuable molecular imaging modality for quantitative measurement of biochemical processes in vivo in the clinical and preclinical imaging domains. However, PET imaging suffers from various physical degrading factors including photon attenuation, which can be corrected using CT-based attenuation correction (CTAC) on combined PET/CT scanners. The attenuation map is calculated by converting CT numbers derived from low-energy polyenergetic x-ray spectra to linear attenuation coefficients at 511 keV. Generation of accurate attenuation maps is crucial for reliable attenuation correction of PET data and hence is a prerequisite for accurate quantification of biological processes. In this study, we implemented the CTAC procedure on the FLEX Triumph™ preclinical PET/CT scanner and evaluated tube voltage dependence for different kVps (40, 50, 60, 70, and 80). The quantitative impact of both bilinear and quadratic based energy-mapping methods on linear attenuation coefficients, attenuation maps and corrected PET images was assessed at different CT tube voltages. Attenuation maps were calculated from CT images of a cylindrical polyethylene phantom containing different concentrations of K2HPO4 in water. Correlation coefficients and best regression fit equations were calculated for both methods. Phantom and rodent PET/CT images were used to assess improvements in image quality and quantitative accuracy. It was observed that the slopes of the bilinear calibration curves for CT numbers greater than 0 HU increase with increasing tube voltage. In addition, higher correlation coefficients were obtained for the quadratic compared to the bilinear energy-mapping method. Tube voltage of 70 kVp produced the smallest relative error and higher correlation coefficient compared to other tube voltages. For low concentrations of K2HPO4, the mean relative difference (in %) between theoretical and calculated attenuation coefficients when using bilinear and quadratic energy-mapping methods are 1.39 ± 1.9 and 1.33 ± 1.8, respectively. They are 2.78 ± 1.3 and 2.5 ± 1.3, respectively, for high concentrations of K2HPO4 . As expected, higher activity concentrations were obtained for PET after attenuation correction. The increased PET signal for mouse tissues ranged between 21 and 31% for bilinear energy-mapping and between 21.8 and 35% for quadratic energy-mapping, whereas these varied from 40 to 51% and from 41 to 56%, respectively, for rat tissues. For biological tissues having a high atomic number such as bone, the quadratic energy-mapping method produced slightly improved results compared to the bilinear energy-mapping method. Phantom and rodent PET studies were successfully corrected for photon attenuation using the developed CTAC procedure.

Monte Carlo dose mapping on deforming anatomy

This paper proposes a Monte Carlo-based energy and mass congruent mapping (EMCM) method to calculate the dose on deforming anatomy. Different from dose interpolation methods, EMCM separately maps each voxel's deposited energy and mass from a source image to a reference image with a displacement vector field (DVF) generated by deformable image registration (DIR). EMCM was compared with other dose mapping methods: energy-based dose interpolation (EBDI) and trilinear dose interpolation (TDI). These methods were implemented in EGSnrc/DOSXYZnrc, validated using a numerical deformable phantom and compared for clinical CT images. On the numerical phantom with an analytically invertible deformation map, EMCM mapped the dose exactly the same as its analytic solution, while EBDI and TDI had average dose errors of 2.5% and 6.0%. For a lung patient's IMRT treatment plan, EBDI and TDI differed from EMCM by 1.96% and 7.3% in the lung patient's entire dose region, respectively. As a 4D Monte Carlo dose calculation technique, EMCM is accurate and its speed is comparable to 3D Monte Carlo simulation. This method may serve as a valuable tool for accurate dose accumulation as well as for 4D dosimetry QA.