Search NASA⌕ Search

SEARCH · Search NASA

Results for “EM Algorithms”

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 109 records · Page 6

PROSPECTS FOR JOINT GRAVITATIONAL-WAVE AND ELECTROMAGNETIC OBSERVATIONS OF NEUTRON-STAR-BLACK-HOLE COALESCING BINARIES

Coalescing neutron-star-black-hole (NS-BH) binaries are a promising source of gravitational-wave (GW) signals detectable with large-scale laser interferometers such as the Advanced Laser Interferometer Gravitational-Wave Observatory and Virgo. They are also one of the main short gamma-ray burst (SGRB) progenitor candidates. If the black hole (BH) tidally disrupts its companion, an SGRB may be ignited when a sufficiently massive accretion disk forms around the remnant BH. Detecting an NS-BH coalescence both in the GW and electromagnetic (EM) spectrum offers a wealth of information about the nature of the source. How much can actually be inferred from a joint detection is unclear, however, as a mass/spin degeneracy may reduce the GW measurement accuracy. To shed light on this problem and on the potential of joint EM+GW observations, we here combine recent semi-analytical predictions for the remnant disk mass with estimates of the parameter-space portion that is selected by a GW detection. We identify cases in which an SGRB ignition is supported, others in which it can be excluded, and finally others in which the outcome depends on the chosen model for the currently unknown NS equation of state. We pinpoint a range of systems that would allow us to place lower bounds on the equation of state stiffness if both the GW emission and its EM counterpart are observed. The methods we develop can broaden the scope of existing GW detection and parameter-estimation algorithms and could allow us to disregard about half of the templates in an NS-BH search following an SGRB trigger, increasing its speed and sensitivity.

79 ASTRONOMY AND ASTROPHYSICS↗

Time-resolved atomic-resolution Brownian tomography of single nanocrystals reveals size-dependent dynamics

Atomic-resolution structure identification of nanocrystals by graphene liquid cell electron microscopy (GLC-EM) has revealed that small, solubilized platinum nanocrystals consist of an ordered crystalline core surrounded by mobile surface atoms, which dissociate during oxidative etching, resulting in distinct temporal structural states. Requirements imposed by the 3D reconstruction algorithm limit the number of structural states that can be resolved. We introduce a regularized 3D reconstruction algorithm that exploits the redundancy inherent in the experimental data, allowing us to improve the time resolution. Our developments provide a comprehensive molecular picture at unprecedented spatial and temporal resolution of the nonlinear, linear, and fluctuating dynamic phenomena that single nanocrystals undergo during the GLC-EM experiment. We determined atomic structures of 66 temporal structural states, extracted from 15 time trajectories of individual nanocrystals. Large (478 to 698 atoms) and small (<300 atoms) nanocrystals show etching that preserves a stable core, whereas mid-sized (351 to 571 atoms) nanocrystals present dynamics that change the coordination of the core.

Meana-Pañeda, Rubén↗

Fully implicit Particle-in-cell algorithms for multiscale plasma simulation

The outline of the paper is as follows: Particle-in-cell (PIC) methods for fully ionized collisionless plasmas, explicit vs. implicit PIC, 1D ES implicit PIC (charge and energy conservation, moment-based acceleration), and generalization to Multi-D EM PIC: Vlasov-Darwin model (review and motivation for Darwin model, conservation properties (energy, charge, and canonical momenta), and numerical benchmarks). The author demonstrates a fully implicit, fully nonlinear, multidimensional PIC formulation that features exact local charge conservation (via a novel particle mover strategy), exact global energy conservation (no particle self-heating or self-cooling), adaptive particle orbit integrator to control errors in momentum conservation, and canonical momenta (EM-PIC only, reduced dimensionality). The approach is free of numerical instabilities: ω pe Δt >> 1, and Δx >> λ D . It requires many fewer dofs (vs. explicit PIC) for comparable accuracy in challenging problems. Significant CPU gains (vs explicit PIC) have been demonstrated. The method has much potential for efficiency gains vs. explicit in long-time-scale applications. Moment-based acceleration is effective in minimizing N FE , leading to an optimal algorithm.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

SU-F-J-10: Sliding Mode Control of a SMA Actuated Active Flexible Needle for Medical Procedures

Purpose: In medical interventional procedures such as brachytherapy, ablative therapies and biopsy precise steering and accurate placement of needles are very important for anatomical obstacle avoidance and accurate targeting. This study presents the efficacy of a sliding mode controller for Shape Memory Alloy (SMA) actuated flexible needle for medical procedures. Methods: Second order system dynamics of the SMA actuated active flexible needle was used for deriving the sliding mode control equations. Both proportional-integral-derivative (PID) and adaptive PID sliding mode control (APIDSMC) algorithms were developed and implemented. The flexible needle was attached at the end of a 6 DOF robotic system. Through LabView programming environment, the control commands were generated using the PID and APIDSMC algorithms. Experiments with artificial tissue mimicking phantom were performed to evaluate the performance of the controller. The actual needle tip position was obtained using an electromagnetic (EM) tracking sensor (Aurora, NDI, waterloo, Canada) at a sampling period of 1ms. During experiment, external disturbances were created applying force and thermal shock to investigate the robustness of the controllers. Results: The root mean square error (RMSE) values for APIDSMC and PID controllers were 0.75 mm and 0.92 mm, respectively, for sinusoidal reference input. In the presence of external disturbances, the APIDSMC controller showed much smoother and less overshooting response compared to that of the PID controller. Conclusion: Performance of the APIDSMC was superior to the PID controller. The APIDSMC was proved to be more effective controller in compensating the SMA uncertainties and external disturbances with clinically acceptable thresholds.

60 APPLIED LIFE SCIENCES↗

Algoritmo regional para estimativa do material particulado em suspensão no baixo Amazonas

The goal of this study was to assess which spectral intervals are better suitable to estimate the suspended particulate matter concentration ([SPM]) as well as test two regional algorithms to estimate [SPM] at the Lower Amazon. Radiometric in situ measurements and surface water sampling to determine the [SPM] were carried during 2014-2017. The best spectral interval to estimate [SPM] was 655-778 nm when all Lower Amazon samples are considered (R2 > 0.70, p<0.005). Among the tested algorithms, a single wavelength non-linear algorithm with the calibrated coefficients for the study area had the best performance (R2 = 0.81; RMSE = 15.3 mg L-1).

Amazon↗

Networked Microgrids for Self-healing Power Systems

This paper proposes a transformative architecture for the normal operation and self-healing of networked microgrids (MGs). MGs can support and interchange electricity with each other in the proposed infrastructure. The networked MGs are connected by a physical common bus and a designed two-layer cyber communication network. The lower layer is within each MG where the energy management system (EMS) schedules the MG operation; the upper layer links a number of EMSs for global optimization and communication. In the normal operation mode, the objective is to schedule dispatchable distributed generators (DGs), energy storage systems (ESs) and controllable loads to minimize the operation costs and maximize the supply adequacy of each MG. When a generation deficiency or fault happens in a MG, the model switches to the self-healing mode and the local generation capacities of other MGs can be used to support the on-emergency portion of the system. A consensus algorithm is used to distribute portions of the desired power support to each individual MG in a decentralized way. The allocated portion corresponds to each MG’s local power exchange target which is used by its EMS to perform the optimal schedule. The resultant aggregated power output of networked MGs will be used to provide the requested power support. Test cases demonstrate the effectiveness of the proposed methodology.

Consensus Algorithm↗

From PINNs to PIKANs: recent advances in physics-informed machine learning

Physics-Informed Neural Networks (PINNs) have emerged as a key tool in Scientific Machine Learning since their introduction in 2017, enabling the efficient solution of ordinary and partial differential equations using sparse measurements. Over the past few years, significant advancements have been made in the training and optimization of PINNs, covering aspects such as network architectures, adaptive refinement, domain decomposition, and the use of adaptive weights and activation functions. A notable recent development is the Physics-Informed Kolmogorov-Arnold Networks (PIKANS), which leverage a representation model originally proposed by Kolmogorov in 1957, offering a promising alternative to traditional PINNs. In this review, we provide a comprehensive overview of the latest advancements in PINNs, focusing on improvements in network design, feature expansion, optimization techniques, uncertainty quantification, and theoretical insights. We also survey key applications across a range of fields, including biomedicine, fluid and solid mechanics, geophysics, dynamical systems, heat transfer, chemical engineering, and beyond. Lastly, we review computational frameworks and software tools developed by both academia and industry to support PINN research and applications.

Kolmogorov-Arnold networks↗

Circuit model of the ITER-like antenna for JET and simulation of its control algorithms

The ITER-like Antenna (ILA) for JET [1] is a 2 toroidal by 2 poloidal array of Resonant Double Loops (RDL) featuring in-vessel matching capacitors feeding RF current straps in conjugate-T manner, a low impedance quarter-wave impedance transformer, a service stub allowing hydraulic actuator and water cooling services to reach the aforementioned capacitors and a 2nd stage phase-shifter-stub matching circuit allowing to correct/choose the conjugate-T working impedance. Toroidally adjacent RDLs are fed from a 3dB hybrid splitter. It has been operated at 33, 42 and 47MHz on plasma (2008-2009) while it presently estimated frequency range is from 29 to 49MHz. At the time of the design (2001-2004) as well as the experiments the circuit models of the ILA were quite basic. The ILA front face and strap array Topica model was relatively crude and failed to correctly represent the poloidal central septum, Faraday Screen attachment as well as the segmented antenna central septum limiter. The ILA matching capacitors, T-junction, Vacuum Transmission Line (VTL) and Service Stubs were represented by lumped circuit elements and simple transmission line models. The assessment of the ILA results carried out to decide on the repair of the ILA identified that achieving routine full array operation requires a better understanding of the RF circuit, a feedback control algorithm for the 2nd stage matching as well as tighter calibrations of RF measurements. The paper presents the progress in modelling of the ILA comprising a more detailed Topica model of the front face for various plasma Scrape Off Layer profiles, a comprehensive HFSS model of the matching capacitors including internal bellows and electrode cylinders, 3D-EM models of the VTL including vacuum ceramic window, Service stub, a transmission line model of the 2nd stage matching circuit and main transmission lines including the 3dB hybrid splitters. A time evolving simulation using the improved circuit model allowed to design and simulate the effectiveness of a feedback control algorithm for the 2nd stage matching and demonstrates the simultaneous matching and control of the 4 RDLs: 11 feedback loops control 21 actuators (8 capacitors, 4 phase shifters and 4 stubs for the 2nd stage matching, 4 main phase shifters controlling of the toroidal phasing and the electronically controlled phase between RF sources feeding top and bottom parts of the array and determines the poloidal phasing of the array which is solved explicitly at each time step) on (simulated) ELMy plasmas.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Proton radiography inversions with source extraction and comparison to mesh methods

In this article, proton radiography is a central diagnostic technique for measuring electromagnetic (EM) fields in high-energy-density, laser-produced plasmas. In this technique, protons traverse the plasma where they accumulate small EM deflections which lead to variations in the proton fluence pattern on a detector. Path-integrated EM fields can then be extracted from the fluence image through an inversion process. In this work, experiments of laser-driven foils were conducted on the OMEGA laser and magnetic field reconstructions were performed using both “fluence-based” techniques and high-fidelity “mesh-based” methods. We implement nonzero boundary conditions into the inversion and show their importance by comparing against mesh measurements. Good agreement between the methods is found only when nonzero boundary conditions are used. We also introduce an approach to determine the unperturbed proton source profile, which is a required input in fluence reconstruction algorithms. In this approach, a fluence inversion is embedded inside of a mesh region, which provides overconstrained magnetic boundary conditions. A source profile is then iteratively optimized to satisfy the boundary information. This method substantially enhances the accuracy in recovering EM fields. Lastly, we propose a scheme to quantify uncertainty in the final inversion that is introduced through errors in the source retrieval.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Reconstruction Performance Studies with MicroBooNE Data in Support of Summer 2018 Analyses

In this note we present a number of reconstruction performance measurements conducted on MicroBooNE data. Results show good agreement between the data and the predictions from simulation, thus validating the version of the reconstruction used by MicroBooNE analyses, which are presented at 2018 conferences. The observed level of performance shows the high level of resolution obtainable with a large-scale LArTPC imaging detector. Tracking and vertexing resolutions at or below the millimeter scale are obtained for neutrino interaction event reconstruction. Energy reconstruction measurements are presented for tracks and EM showers using multiple Coulomb scattering and calorimetry respectively; in the case of EM showers these represent the first result obtained in a fully-automated reconstruction chain. Nevertheless, the present results are expected to improve with further refinements of the employed reconstruction algorithms. The techniques which are presented here have a broader validity than MicroBooNE and could be applied to other LArTPC detectors as they start taking data in the next months or years.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Hierarchical Control of Megawatt-Scale Charging Stations for Electric Trucks with Distributed Energy Resources

Electrifying medium- and heavy-duty trucks is critical to decarbonizing the transportation sector. Energy needs of electric trucks will likely require megawatt-scale charging stations, which could significantly stress the electric distribution grid. Distributed energy resources (DER) can alleviate this stress and reduce charging costs with proper management. To that end, this work develops a hierarchical predictive control algorithm for future multi-port megawatt-scale charging stations that can provide real-time energy management for stations, decide charging rates, dispatch energy storage system (ESS), and provide grid voltage support. We integrate three algorithmic components: (i) an energy management optimization (EMO) that provides supervisory control to DER assets and charging loads at minute scale, (ii) a real-time energy management system (RT-EMS) that heuristically compensates for fast disturbances at sub-second scale, and (iii) a model predictive control (MPC)-based battery management system (BMS) that communicates future charging demands to the EMO, to manage the overall megawatt-scale site. Additionally, validation in a controller hardware-in-the-loop (CHIL) environment shows that the hierarchical controller can reduce the total energy consumption from the grid by approximately 28% compared to an uncontrolled case for the station configuration in this paper, without impacting charging time.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Cryo-electron tomography related radiation-damage parameters for individual-molecule 3D structure determination

To understand the dynamic structure–function relationship of soft- and biomolecules, the determination of the three-dimensional (3D) structure of each individual molecule (nonaveraged structure) in its native state is sought-after. Cryo-electron tomography (cryo-ET) is a unique tool for imaging an individual object from a series of tilted views. However, due to radiation damage from the incident electron beam, the tolerable electron dose limits image contrast and the signal-to-noise ratio (SNR) of the data, preventing the 3D structure determination of individual molecules, especially at high-resolution. Although recently developed technologies and techniques, such as the direct electron detector, phase plate, and computational algorithms, can partially improve image contrast/SNR at the same electron dose, the high-resolution structure, such as tertiary structure of individual molecules, has not yet been resolved. Here, we review the cryo-electron microscopy (cryo-EM) and cryo-ET experimental parameters to discuss how these parameters affect the extent of radiation damage. This discussion can guide us in optimizing the experimental strategy to increase the imaging dose or improve image SNR without increasing the radiation damage. With a higher dose, a higher image contrast/SNR can be achieved, which is crucial for individual-molecule 3D structure. With 3D structures determined from an ensemble of individual molecules in different conformations, the molecular mechanism through their biochemical reactions, such as self-folding or synthesis, can be elucidated in a straightforward manner.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Effects of a nuclear-disturbed environment on electromagnetic wave propagation through the atmosphere

This paper investigates the effects of a nuclear-disturbed environment on the transmission of electromagnetic (EM) waves through the atmosphere. An atmospheric nuclear detonation can produce heightened free electron densities in the surrounding atmosphere that can disrupt EM waves that propagate through the disturbed region. Radiation transport models simulated the ionization and free electron densities created in the atmosphere from a 1 MT detonation at heights of burst of 5 km, 25 km, and 75 km. Recombination rates for the free electrons in the atmosphere were applied, from previous work in the literature, to determine the nuclear-induced electron densities as a function of time and space after the detonation. A ray-tracing algorithm was applied to determine the refraction and reflection of waves propagating in the different nuclear-disturbed environments. The simulation results show that the free electron plasma created from an atmospheric nuclear detonation depend on the height of burst of the weapon, the weapon yield, and the time after detonation. Detonations at higher altitudes produce higher free electron densities for greater durations and over larger ranges. The larger the free electron densities, the greater the impact on EM wavelengths in regards to refraction, reflection, and absorption in the atmosphere. An analysis of modern infrastructure and the effects of nuclear-disturbed atmospheres on different signal wavelengths and systems is discussed.

42 ENGINEERING↗

On coincident loop transient electromagnetic induction logging

Coincident loop transient induction wireline logging is examined as the borehole analog of the well-known land and airborne time-domain electromagnetic (EM) method. The concept of whole-space late-time apparent resistivity is modified from the half-space version commonly used in land and airborne geophysics and applied to the coincident loop voltages produced from various formation, borehole, and invasion models. Given typical tool diameters, off-time measurements with such an instrument must be made on the order of nanoseconds to microseconds — much more rapidly than for surface methods. Departure curves of the apparent resistivity for thin beds, calculated using an algorithm developed to model the transient response of a loop in a multilayered earth, indicate that the depth of investigation scales with the bed thickness. Modeled resistivity logs are comparable in accuracy and resolution with standard frequency-domain focused induction logs. However, if measurement times are longer than a few microseconds, the thicknesses of conductors can be overestimated, whereas resistors are underestimated. Thin-bed resolution characteristics are explained by visualizing snapshots of the EM fields in the formation, where a conductor traps the electric field while two current maxima are produced in the shoulder beds surrounding a resistor. Radial profiling is studied using a concentric cylinder earth model. Results found that true formation resistivity can be determined in the presence of either oil- or water-based mud, although in the latter case, measurements must be taken several orders of magnitude later in time. Lastly, the ability to determine true formation resistivity is governed by the degree that the EM field heals after being distorted by borehole fluid and invasion, a process visualized and particularly evident in the case of conductive water-based mud.

58 GEOSCIENCES↗

On specular reflectivity measurements in high and low-contrast relativistic laser-plasma interactions

Using both experiment and 2D3V particle-in-cell (PIC) simulations, we describe the use of specular reflectivity measurements to study relativistic (Iλ{sup 2 }> 10{sup 18 }W/cm{sup 2}⋅μm{sup 2}) laser-plasma interactions for both high and low-contrast 527 nm laser pulses on initially solid density aluminum targets. In the context of hot-electron generation, studies typically rely on diagnostics which, more-often-than-not, represent indirect processes driven by fast electrons transiting through solid density materials. Specular reflectivity measurements, however, can provide a direct measure of the interaction that is highly sensitive to how the EM fields and plasma profiles, critical input parameters for modeling of hot-electron generation, evolve near the interaction region. While the fields of interest occur near the relativistic critical electron density, experimental reflectivity measurements are obtained centimeters away from the interaction region, well after diffraction has fully manifested itself. Using a combination of PIC simulations with experimentally inspired conditions and an analytic, non-paraxial, pulse propagation algorithm, we calculate reflected pulse properties, both near and far from the interaction region, and compare with specular reflectivity measurements. The experiment results and PIC simulations demonstrate that specular reflectivity measurements are an extremely sensitive qualitative, and partially quantitative, indicator of initial laser/target conditions, ionization effects, and other details of intense laser-matter interactions. The techniques described can provide strong constraints on many systems of importance in ultra-intense laser interactions with matter.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Fractional modeling of viscoelasticity in 3D cerebral arteries and aneurysms

We develop efficient numerical methods for fractional order PDEs, and employ them to investigate viscoelastic constitutive laws for arterial wall mechanics. Recent simulations using one-dimensional models [1] have indicated that fractional order models may offer a more powerful alternative for modeling the arterial wall response, exhibiting reduced sensitivity to parametric uncertainties compared with the integer-calculus-based models. Here, we study three-dimensional (3D) fractional PDEs that naturally model the continuous relaxation properties of soft tissue, and for the first time employ them to simulate flow structure interactions for patient-specific brain aneurysms. To deal with the high memory requirements and in order to accelerate the numerical evaluation of hereditary integrals, we employ a fast convolution method [2] that reduces the memory cost to O(log(N)) and the computational complexity to O(N log(N)). Furthermore, we combine the fast convolution with high-order backward differentiation to achieve third-order time integration accuracy. We confirm that in 3D viscoelastic simulations, the integer order models strongly depends on the relaxation parameters, while the fractional order models are less sensitive. As an application to long-time simulations in complex geometries, we also apply the method to modeling fluid-structure interaction of a 3D patient-specific compliant cerebral artery with an aneurysm. Taken together, our findings demonstrate that fractional calculus can be employed effectively in modeling complex behavior of materials in realistic 3D time-dependent problems if properly designed efficient algorithms are employed to overcome the extra memory requirements and computational complexity associated with the non-local character of fractional derivatives.

Computer Science↗

Advances in Multi-Physics Modeling for Parallel Finite Element Code Suite ACE3P

SLAC developed ACE3P code suite is a comprehensive set of parallel finite-element codes for multi-physics modeling of particle accelerators. Running on massively parallel computer platforms for high fidelity and high accuracy simulation, ACE3P enables rapid virtual prototyping of accelerator and RF component design, optimization and analysis. Furthermore, recent advances of ACE3P have been focused on the implementation of advanced numerical algorithms, enhancement of multi-physics capabilities, and improvement of code performance on state-of-the-art high performance computing (HPC) platforms for large-scale RF modeling in accelerator applications and beyond. In this paper, we will present a nonlinear electromagnetic (EM) eigensolver for damping calculation of resonant modes for cavities coupled with waveguides, the thermal solver for accurate evaluation of heat loads in superconducting RF (SRF) cavity power coupler, the elastic solver for investigation of SRF cavity cryomodule (CM) bowing due to thermal effect, the frequency domain mechanical solver for studying SRF cavity mechanical oscillation modes, the particle tracking code for dark current simulation in SRF CM using a hybrid MPI+OpenMP parallel programing model, the time domain solver for high fidelity modeling of EM wave propagation for medical device in human body.

97 MATHEMATICS AND COMPUTING↗