Search NASA⌕ Search

SEARCH · Search NASA

Results for “error fields”

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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Beam distortion effects upon focusing an ultrashort petawatt laser pulse to greater than 10 22 W/cm 2

When an ultrashort laser pulse is tightly focused to a size approaching its central wavelength, the properties of the focused spot diverge from the diffraction-limited case. Here, we report on this change in behavior of a tightly focused petawatt-class laser beam by an $f/1$ off-axis parabolic mirror (OAPM). Considering the effects of residual aberration, the spatial profile of the near-field, and pointing error, we estimate the deviation in peak intensities of the focused spot from the ideal case. Here, we verify that the estimated peak intensity values are within an acceptable error range of the measured values. With the added uncertainties in target alignment, we extend the estimation to infer on-target peak intensities of ≥10 22 W/cm 2 for a target at the focal plane of this $f/1$ OAPM.

47 OTHER INSTRUMENTATION↗
DOE OSTI2014

SU-E-T-226: Junction Free Craniospinal Irradiation in Linear Accelerator Using Volumetric Modulated Arc Therapy : A Novel Technique Using Dose Tapering

Purpose: Spatially separated fields are required for craniospinal irradiation due to field size limitation in linear accelerator. Field junction shits are conventionally done to avoid hot or cold spots. Our study was aimed to demonstrate the feasibility of junction free irradiation plan of craniospinal irradiation (CSI) for Meduloblastoma cases treated in linear accelerator using Volumetric modulated arc therapy (VMAT) technique. Methods: VMAT was planned using multiple isocenters in Monaco V 3.3.0 and delivered in Elekta Synergy linear accelerator. A full arc brain and 40° posterior arc spine fields were planned using two isocentre for short (<1.3 meter height ) and 3 isocentres for taller patients. Unrestricted jaw movement was used in superior-inferior direction. Prescribed dose to PTV was achieved by partial contribution from adjacent beams. A very low dose gradient was generated to taper the isodoses over a long length (>10 cm) at the conventional field junction. Results: In this primary study five patients were planned and three patients were delivered using this novel technique. As the dose contribution from the adjacent beams were varied (gradient) to create a complete dose distribution, therefore there is no specific junction exists in the plan. The junction were extended from 10–14 cm depending on treatment plan. Dose gradient were 9.6±2.3% per cm for brain and 7.9±1.7 % per cm for spine field respectively. Dose delivery error due to positional inaccuracy was calculated for brain and spine field for ±1mm, ±2mm, ±3mm and ±5 mm were 1%–0.8%, 2%–1.6%, 2.8%–2.4% and 4.3%–4% respectively. Conclusion: Dose tapering in junction free CSI do not require a junction shift. Therefore daily imaging for all the field is also not essential. Due to inverse planning dose to organ at risk like thyroid kidney, heart and testis can be reduced significantly. VMAT gives a quicker delivery than Step and shoot or dynamic IMRT.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2023

Accurate simulation of direct laser acceleration in a laser wakefield accelerator

In a laser wakefield accelerator (LWFA), an intense laser pulse excites a plasma wave that traps and accelerates electrons to relativistic energies. When the pulse overlaps the accelerated electrons, it can enhance the energy gain through direct laser acceleration (DLA) by resonantly driving the betatron oscillations of the electrons in the plasma wave. The traditional particle-in-cell (PIC) algorithm, although often the tool of choice to study DLA, contains inherent errors due to numerical dispersion and the time staggering of the electric and magnetic fields. Furthermore, conventional PIC implementations cannot reliably disentangle the fields of the plasma wave and laser pulse, which obscures interpretation of the dominant acceleration mechanism. Here, a customized field solver that reduces errors from both numerical dispersion and time staggering is used in conjunction with a field decomposition into azimuthal modes to perform PIC simulations of DLA in an LWFA. Comparisons with traditional PIC methods, model equations, and experimental data show improved accuracy with the customized solver and convergence with an order-of-magnitude fewer cells. Furthermore, the azimuthal-mode decomposition reveals that the most energetic electrons receive comparable energy from DLA and LWFA.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Human performance benefits gained by dynamic instructions compared to smart PDFs

Human performance benefits gained by dynamic instructions compared to smart PDFs In recent years multiple nuclear utilities have deployed an electronic work management solution to replace the traditional paper-based processes. There are many benefits of electronic work management solutions compared to paper processes, such as a more streamlined review and approval process of work packages and a better way to track and resolve any potential bottlenecks causing delays in the review process. These benefits have been recognized where electronic work management solutions have been successfully deployed. However the utilities also hoped the new solution would help the field worker conduct his/her job more efficiently, correctly, and safely. Unfortunately, the industry is now realizing that the deployed solutions did not fully help them reach that goal. The electronic work management solutions present the work instructions as an electronic copy (i.e., a portable document format (PDF)) of the paper instruction. The PDF versions used in these electronic work management solutions have some ability for data input, notes, and mark-ups, but are unable to provide many additional benefits compared to using paper. In fact, conducting work using a PDF on a tablet can be more cumbersome than with a paper instruction. For example, flipping back and forth in a paper instruction can be time consuming, but it is not at all as frustrating as scrolling back and forth on a tablet or continuously having to navigate between tabs. The researchers at Idaho National Laboratory (INL) approached the next step beyond smart PDFs from a human factors perspective. In close collaboration with multiple utilities they investigated how to present instruction content and how to design the interaction with the instruction to improve the workers performance and efficiency. Their solution is a dynamic presentation of the instruction that guides the worker through the correct path as well as ensuring the worker is on the correct component. The solution also automatically validates input and provides guidance if unexpected conditions are identified. A well-designed solution has the opportunity to vastly reduce the risk of human errors in the field, such as preventing actions on incorrect component or unintentionally skipping steps. The human factors researchers at INL leveraged their expertise, existing human performance tools, and a great amount of involvement from field workers in multiple utilities to develop and evaluate design concepts that will provide the human performance and efficiency benefits that the electronic work management solutions fall short to deliver. This paper will describe in detail how these new design concepts improve human performance.

47 OTHER INSTRUMENTATION↗
DOE OSTI2013

Magnetic field calculation for a 10 MeV positron emission tomography cyclotron

The magnetic field calculation and correction for a 10 MeV positron emission tomography cyclotron is presented. 3D TOSCA analysis results are compared with the measured data, and the calculation error is used to calibrate the B-H curve to obtain a very precise finite element method estimator, which is used to predict the correction of the magnet pole for achieving the isochronous field. The isochronous field error is approximated with the effects of a set of standard patches. On the assumption that the effect of each small patch is proportional to its surface, the correction of the magnet pole is found by solving a system of equations using the least square scheme. The magnet shimming is performed and the measured magnetic field is found in good agreement with the prediction, with an error less than 2 G.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2021

A Brief Account of Steven Weinberg’s Legacy in ab initio Many-Body Theory

In this contribution to the special issue “Celebrating 30 years of Steven Weinberg’s papers on Nuclear Forces from Chiral Lagrangians,” we emphasize the important role chiral effective field theory has played in leading nuclear physics into a precision era. To this end, we share our perspective on a few of the recent advances made in ab initio calculations of nuclear structure and nuclear matter observables, as well as Bayesian uncertainty quantification of effective field theory truncation errors.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Quantifying Progress Across Different CMIP Phases With the ESMValTool

Abstract More than 40 model groups worldwide are participating in the Coupled Model Intercomparison Project Phase 6 (CMIP6), providing a new and rich source of information to better understand past, present, and future climate change. Here, we use the Earth System Model Evaluation Tool (ESMValTool) to assess the performance of the CMIP6 ensemble compared to the previous generations CMIP3 and CMIP5. While CMIP5 models did not capture the observed pause in the increase in global mean surface temperature between 1998 and 2013, the historical CMIP6 simulations agree well with the observed recent temperature increase, but some models have difficulties in reproducing the observed global mean surface temperature record of the second half of the twentieth century. While systematic biases in annual mean surface temperature and precipitation remain in the CMIP6 multimodel mean, individual models and high‐resolution versions of the models show significant reductions in many long‐standing biases. Some improvements are also found in the vertical temperature, water vapor, and zonal wind speed distributions, and root‐mean‐square errors for selected fields are generally smaller with reduced intermodel spread and higher average skill in the correlation patterns relative to observations. An emerging property of the CMIP6 ensemble is a higher effective climate sensitivity with an increased range between 2.3 and 5.6 K. A possible reason for this increase in some models is improvements in cloud representation resulting in stronger shortwave cloud feedbacks than in their predecessor versions.

Bock, L.↗
DOE OSTI2023

Distributed Ledger Technology for Fault Tolerant Distribution Grid Operations

This paper explores the potential of distributed ledger technology (DLT) to improve fault-tolerant grid operations by leveraging its core features as an immutable, decentralized ledger, a distributed, consensus-based agreement process, and a distributed state-replication engine. Distribution power systems deliver electricity to millions of customers; however, they are susceptible to various threats that can result in customer interruptions. These include faults caused by adverse weather conditions, natural disasters, vegetation growth, equipment failure, and malicious attacks. To minimize the effects of these faults, fault-handling approaches rely on network knowledge to isolate affected areas and reconnect unaffected areas, reducing the number of affected customers while maintaining safety. Here, we present a trusted data-sharing architecture that enables independent, distributed actors to reconstruct the pre-fault system state by enabling distributed resources to make appropriate decisions with limited network/system information. Although the process requires some data sharing between switch-delimited areas, the approach limits the amount of private information shared, preserving customers' privacy and business-sensitive information. We include three use cases that form a foundation for third parties to develop functional solutions that can eventually be deployed in the field. The gross error detection method used within switch-delimited areas can identify sensor errors and accurately detect circuit breaker states. The evaluation of possible reconnection while preserving data ownership resulted in a voltage magnitude difference smaller than 0.001% from the OpenDSS power flow solution that has full system knowledge, which is below the expected power flow tolerance. The approach offers a promising opportunity for improving fault-tolerant distribution grid operations.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2023

Q-BEEP: Quantum Bayesian Error Mitigation Employing Poisson Modeling over the Hamming Spectrum

Quantum computing technology has grown rapidly in recent years, with new technologies being explored, error rates being reduced, and quantum processor’s qubit capacity growing. However, near-term quantum algorithms are still unable to be induced without compounding consequential levels of noise, leading to non-trivial erroneous results. Quantum Error Correction (in-situ error mitigation) and Quantum Error Mitigation (post-induction error mitigation) are promising fields of research within the quantum algorithm scene, aiming to alleviate quantum errors, increasing the overall fidelity and hence the overall quality of circuit induction. Earlier this year, a pioneering work, namely HAMMER, published in ASPLOS-22 demonstrated the existence of a latent structure regarding post-circuit induction errors when mapping to the Hamming spectrum. However, they intuitively assumed that errors occur in local clusters, and that at higher average Hamming distances this structure falls away. In this work, we show that such a correlation structure is not only local but extends certain non-local clustering patterns which can be precisely described by a Poisson distribution model taking the input circuit, the device run time status (i.e., calibration statistics) and qubit topology into consideration. Using this quantum error characterizing model, we developed an iterative algorithm over the generated Bayesian network state-graph for post-induction error mitigation. Thanks to more precise modeling of the error distribution latent structure and the new iterative method, our Q-Beep approach provides state of the art performance and can boost circuit execution fidelity by up to 234.6% on Bernstein-Vazirani circuits and on average 71.0% on QAOA solution quality, using 16 practical IBMQ quantum processors. For other benchmarks such as those in QASMBench, the fidelity improvement is up to 17.8%. Q-Beep is a light-weight post-processing technique that can be performed offline and remotely, making it a useful tool for quantum vendors to integrate and provide more reliable circuit induction results.

Stein, Samuel A.↗
DOE OSTI2024

Prediction of laser beam spatial profiles in a high-energy laser facility by use of deep learning

We adapt the significant advances achieved recently in the field of generative artificial intelligence/machine-learning to laser performance modeling in multipass, high-energy laser systems with application to high-shot-rate facilities relevant to inertial fusion energy. Advantages of neural-network architectures include rapid prediction capability, data-driven processing, and the possibility to implement such architectures within future low-latency, low-power consumption photonic networks. Four models were investigated that differed in their generator loss functions and utilized the U-Net encoder/decoder architecture with either a reconstruction loss alone or combined with an adversarial network loss. We achieved inference times of 1.3 ms for a 256 × 256 pixel near-field beam with errors in predicted energy of the order of 1% over most of the energy range. It is shown that prediction errors are significantly reduced by ensemble averaging the models with different weight initializations. These results suggest that including the temporal dimension in such models may provide accurate, real-time spatiotemporal predictions of laser performance in high-shot-rate laser systems.

47 OTHER INSTRUMENTATION↗
DOE OSTI2020

Predicting secondary organic aerosol phase state and viscosity and its effect on multiphase chemistry in a regional-scale air quality model

Atmospheric aerosols are a significant public health hazard and have substantial impacts on the climate. Secondary organic aerosols (SOAs) have been shown to phase separate into a highly viscous organic outer layer surrounding an aqueous core. This phase separation can decrease the partitioning of semi-volatile and low-volatile species to the organic phase and alter the extent of acid-catalyzed reactions in the aqueous core. A new algorithm that can determine SOA phase separation based on their glass transition temperature (T g ), oxygen to carbon (O:C) ratio and organic mass to sulfate ratio, and meteorological conditions was implemented into the Community Multiscale Air Quality Modeling (CMAQ) system version 5.2.1 and was used to simulate the conditions in the continental United States for the summer of 2013. SOA formed at the ground/surface level was predicted to be phase separated with core–shell morphology, i.e., aqueous inorganic core surrounded by organic coating 65.4 % of the time during the 2013 Southern Oxidant and Aerosol Study (SOAS) on average in the isoprene-rich southeastern United States. Our estimate is in proximity to the previously reported ~70 % in literature. The phase states of organic coatings switched between semi-solid and liquid states, depending on the environmental conditions. The semi-solid shell occurring with lower aerosol liquid water content (western United States and at higher altitudes) has a viscosity that was predicted to be 10 2 –10 12 Pa s, which resulted in organic mass being decreased due to diffusion limitation. Organic aerosol was primarily liquid where aerosol liquid water was dominant (eastern United States and at the surface), with a viscosity <10 2 Pa s. Phase separation while in a liquid phase state, i.e., liquid–liquid phase separation (LLPS), also reduces reactive uptake rates relative to homogeneous internally mixed liquid morphology but was lower than aerosols with a thick viscous organic shell. The sensitivity cases performed with different phase-separation parameterization and dissolution rate of isoprene epoxydiol (IEPOX) into the particle phase in CMAQ can have varying impact on fine particulate matter (PM 2.5 ) organic mass, in terms of bias and error compared to field data collected during the 2013 SOAS. This highlights the need to better constrain the parameters that govern phase state and morphology of SOA, as well as expand mechanistic representation of multiphase chemistry for non-IEPOX SOA formation in models aided by novel experimental insights.

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

MATCHED-INDEX-OF-REFRACTION FLOW FACILITY FOR FUNDAMENTAL AND APPLIED RESEARCH

Significant challenges face reactor designers with regard to thermal hydraulic design and associated modeling for advanced reactor concepts. Computational thermal hydraulic codes solve only a piece of the core. There is a need for a whole core dynamics system code with local resolution to investigate and understand flow behavior with all the relevant physics and thermo-mechanics. The matched index of refraction (MIR) flow facility at Idaho National Laboratory (INL) has a unique capability to contribute to the development of validated computational fluid dynamics (CFD) codes through the use of state-of-the-art optical measurement techniques, such as Laser Doppler Velocimetry (LDV) and Particle Image Velocimetry (PIV). PIV is a non-intrusive velocity measurement technique that tracks flow by imaging the movement of small tracer particles within a fluid. At the heart of a PIV calculation is the cross correlation algorithm, which is used to estimate the displacement of particles in some small part of the image over the time span between two images. Generally, the displacement is indicated by the location of the largest peak. To quantify these measurements accurately, sophisticated processing algorithms correlate the locations of particles within the image to estimate the velocity (Ref. 1). Prior to use with reactor deign, the CFD codes have to be experimentally validated, which requires rigorous experimental measurements to produce high quality, multi-dimensional flow field data with error quantification methodologies. Computational thermal hydraulic codes solve only a piece of the core. There is a need for a whole core dynamics system code with local resolution to investigate and understand flow behavior with all the relevant physics and thermo-mechanics. Computational techniques with supporting test data may be needed to address the heat transfer from the fuel to the coolant during the transition from turbulent to laminar flow, including the possibility of an early laminarization of the flow (Refs. 2 and 3) (laminarization is caused when the coolant velocity is theoretically in the turbulent regime, but the heat transfer properties are indicative of the coolant velocity being in the laminar regime). Such studies are complicated enough that computational fluid dynamics (CFD) models may not converge to the same conclusion. Thus, experimentally scaled thermal hydraulic data with uncertainties should be developed to support modeling and simulation for verification and validation activities. The fluid/solid index of refraction matching technique allows optical access in and around geometries that would otherwise be impossible while the large test section of the INL system provides better spatial and temporal resolution than comparable facilities. Benchmark data for assessing computational fluid dynamics can be acquired for external flows, internal flows, and coupled internal/external flows for better understanding of physical phenomena of interest. The core objective of this study is to describe MIR and its capabilities, and mention current development areas for uncertainty quantification, mainly the uncertainty surface method and cross-correlation method. Using these methods, it is anticipated to establish a suitable approach to quantify PIV uncertainty for experiments performed in the MIR.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

Towards Feature-Preserving 2D and 3D Vector Field Compression

The objective of this work is to develop error-bounded lossy compression methods to preserve topological features in 2D and 3D vector fields. Specifically, we explore the preservation of critical points in piecewise linear vector fields. We define the preservation of critical points as, without any false positive, false negative, or false type change in the decompressed data, (1) keeping each critical point in its original cell and (2) retaining the type of each critical point (e.g., saddle and attracting node). The key to our method is to adapt a vertex-wise error bound for each grid point and to compress input data together with the error bound field using a modified lossy compressor. Our compression algorithm can be also embarrassingly parallelized for large data handling and in situ processing. We benchmark our method by comparing it with existing lossy compressors in terms of false positive/negative/type rates, compression ratio, and various vector field visualizations with several scientific applications.

critical points↗
DOE OSTI2011

EVOLUTION OF GALAXY LUMINOSITY FUNCTION USING PHOTOMETRIC REDSHIFTS

We examine the impact of using photometric redshifts for studying the evolution of both the global galaxy luminosity function (LF) and that for different galaxy types. To this end, we compare the LFs obtained using photometric redshifts from the Canada-France-Hawaii Telescope Legacy Survey (CFHTLS) D1 field with those from the spectroscopic survey VIMOS VLT Deep Survey (VVDS) comprising {approx}4800 galaxies. We find that for z {<=} 2.0, in the interval of magnitudes considered by this survey, the LFs obtained using photometric and spectroscopic redshifts show a remarkable agreement. This good agreement led us to use all four Deep fields of the CFHTLS comprising {approx}386,000 galaxies to compute the LF of the combined fields and directly estimate the error in the parameters based on the field-to-field variation. We find that the characteristic absolute magnitude M* of Schechter fits fades by {approx}0.7 mag from z {approx} 1.8 to z {approx} 0.3, while the characteristic density {phi}* increases by a factor of {approx}4 in the same redshift interval. We use the galaxy classification provided by the template fitting program used to compute photometric redshifts and split the sample into galaxy types. We find that these Schechter parameters evolve differently for each galaxy type, an indication that their evolution is a combination of several effects: galaxy merging, star formation quenching, and mass assembly. All these results are compatible with those obtained by different spectroscopic surveys such as VVDS, DEEP2, and zCosmos, which reinforces the fact that photometric redshifts can be used to study galaxy evolution, at least for the redshift bins adopted so far. This is of great interest since future very large imaging surveys containing hundreds of millions of galaxies will allow us to obtain important precise measurements to constrain the evolution of the LF and to explore the dependence of this evolution on morphology and/or color helping constrain the mechanisms of galaxy evolution.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

A deep learning-based workflow for fast prediction of 3D state variables in geological carbon storage: A dimension reduction approach

Deep learning (DL) models are extensively used as surrogate models for high-fidelity simulations of multiphase fluid flow in porous media at large scales, enabling fast forecasts of the spatial–temporal evolution of three-dimensional (3D) state variables in geological carbon storage (GCS). However, training these models in high-dimensional space remains computationally demanding and prone to overfitting because of limited training data. This paper presents a novel workflow to address these challenges by integrating dimension reduction (DR) methods. Here, the proposed workflow employed pre-trained DR models to extract the latent variables of geological models and state variables and utilized the multi-layer perceptron (MLP) for constructing mapping functions between the input and output variables in latent spaces. Subsequently, the pre-trained reconstruction models converted the MLP-predicted latent state variables to their original high-dimensional form. Furthermore, we proposed a novel strategy for the DR and reconstruction of 3D saturation fields to account for the unique data characteristics of sparsity, nonuniformity, and discontinuity. The proposed strategy applied PCA and inverse PCA for 2D average saturation fields and developed a DL-based 3D reconstruction model, leveraging three 2D average saturation fields as input to produce a 3D saturation field as output. The pre-training of DR and reconstruction models and training of MLP models were conducted on 84 Gulf of Mexico (GoM) simulations and evaluated on 12 testing simulations. Each simulation contained 720 monthly time steps, with the first 360 months as the injection period and the rest as the post-injection period. The proposed workflow, incorporating DR and DL models, accurately predicts the normalized 3D pressure fields, achieving mean square error (MSE) of 2.92 × 10 -7 compared to the ground truth obtained from a full-physics simulator. Furthermore, the proposed strategy outperformed PCA and convolutional autoencoder (CAE) models on 3D saturation fields, resulting in minor workflow prediction errors with an MSE of 2.93 × 10 -5 . The results suggest the proposed workflow provides sufficient predictive fidelity across temporal and spatial scales, and enables a speedup of 160 times compared to the full-physics simulator, facilitating improved decision-making and risk assessment for large-scale GCS management in real-time scenarios.

3D reconstruction model↗
DOE OSTI2017

Checking Questionable Entry of Personally Identifiable Information Encrypted by One-Way Hash Transformation

Background: As one of the several effective solutions for personal privacy protection, a global unique identifier (GUID) is linked with hash codes that are generated from combinations of personally identifiable information (PII) by a one-way hash algorithm. On the GUID server, no PII is permitted to be stored, and only GUID and hash codes are allowed. The quality of PII entry is critical to the GUID system. Objective: The goal of our study was to explore a method of checking questionable entry of PII in this context without using or sending any portion of PII while registering a subject. Methods: According to the principle of GUID system, all possible combination patterns of PII fields were analyzed and used to generate hash codes, which were stored on the GUID server. Based on the matching rules of the GUID system, an error-checking algorithm was developed using set theory to check PII entry errors. We selected 200,000 simulated individuals with randomly-planted errors to evaluate the proposed algorithm. These errors were placed in the required PII fields or optional PII fields. The performance of the proposed algorithm was also tested in the registering system of study subjects. Results: There are 127,700 error-planted subjects, of which 114,464 (89.64%) can still be identified as the previous one and remaining 13,236 (10.36%, 13,236/127,700) are discriminated as new subjects. As expected, 100% of nonidentified subjects had errors within the required PII fields. The possibility that a subject is identified is related to the count and the type of incorrect PII field. For all identified subjects, their errors can be found by the proposed algorithm. The scope of questionable PII fields is also associated with the count and the type of the incorrect PII field. The best situation is to precisely find the exact incorrect PII fields, and the worst situation is to shrink the questionable scope only to a set of 13 PII fields. In the application, the proposed algorithm can give a hint of questionable PII entry and perform as an effective tool. Conclusions: The GUID system has high error tolerance and may correctly identify and associate a subject even with few PII field errors. Correct data entry, especially required PII fields, is critical to avoiding false splits. In the context of one-way hash transformation, the questionable input of PII may be identified by applying set theory operators based on the hash codes. The count and the type of incorrect PII fields play an important role in identifying a subject and locating questionable PII fields.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2017

Magnus: A New Resistive MHD Code with Heat Flow Terms

We present a new magnetohydrodynamic (MHD) code for the simulation of wave propagation in the solar atmosphere, under the effects of electrical resistivity—but not dominant—and heat transference in a uniform 3D grid. The code is based on the finite-volume method combined with the HLLE and HLLC approximate Riemann solvers, which use different slope limiters like MINMOD, MC, and WENO5. In order to control the growth of the divergence of the magnetic field, due to numerical errors, we apply the Flux Constrained Transport method, which is described in detail to understand how the resistive terms are included in the algorithm. In our results, it is verified that this method preserves the divergence of the magnetic fields within the machine round-off error (∼1×10{sup −12}). For the validation of the accuracy and efficiency of the schemes implemented in the code, we present some numerical tests in 1D and 2D for the ideal MHD. Later, we show one test for the resistivity in a magnetic reconnection process and one for the thermal conduction, where the temperature is advected by the magnetic field lines. Moreover, we display two numerical problems associated with the MHD wave propagation. The first one corresponds to a 3D evolution of a vertical velocity pulse at the photosphere–transition–corona region, while the second one consists of a 2D simulation of a transverse velocity pulse in a coronal loop.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

The magnetic field of active region 11158 during the 2011 February 12-17 flares: Differences between photospheric extrapolation and coronal forward-fitting methods

We developed a coronal nonlinear force-free field (COR-NLFFF) forward-fitting code that fits an approximate nonlinear force-free field (NLFFF) solution to the observed geometry of automatically traced coronal loops. In contrast to photospheric NLFFF codes, which calculate a magnetic field solution from the constraints of the transverse photospheric field, this new code uses coronal constraints instead, and this way provides important information on systematic errors of each magnetic field calculation method, as well as on the non-force-freeness in the lower chromosphere. In this study we applied the COR-NLFFF code to NOAA Active Region 11158, during the time interval of 2011 February 12-17, which includes an X2.2 GOES-class flare plus 35 M- and C-class flares. We calculated the free magnetic energy with a 6 minute cadence over 5 days. We find good agreement between the two types of codes for the total nonpotential E{sub N} and potential energy E{sub P} but find up to a factor of 4 discrepancy in the free energy E {sub free} = E{sub N} – E{sub P} and up to a factor of 10 discrepancy in the decrease of the free energy ΔE {sub free} during flares. The coronal NLFFF code exhibits a larger time variability and yields a decrease of free energy during the flare that is sufficient to satisfy the flare energy budget, while the photospheric NLFFF code shows much less time variability and an order of magnitude less free-energy decrease during flares. The discrepancy may partly be due to the preprocessing of photospheric vector data but more likely is due to the non-force-freeness in the lower chromosphere. We conclude that the coronal field cannot be correctly calculated on the basis of photospheric data alone and requires additional information on coronal loop geometries.

79 ASTRONOMY AND ASTROPHYSICS↗