Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Maxwell method”

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

Sandtank-ML: An Educational Tool at the Interface of Hydrology and Machine Learning

Hydrologists and water managers increasingly face challenges associated with extreme climatic events. At the same time, historic datasets for modeling contemporary and future hydrologic conditions are increasingly inadequate. Machine learning is one promising technological tool for navigating the challenges of understanding and managing contemporary hydrological systems. However, in addition to the technical challenges associated with effectively leveraging ML for understanding subsurface hydrological processes, practitioner skepticism and hesitancy surrounding ML presents a significant barrier to adoption of ML technologies among practitioners. In this paper, we discuss an educational application we have developed—Sandtank-ML—to be used as a training and educational tool aimed at building user confidence and supporting adoption of ML technologies among water managers. We argue that supporting the adoption of ML methods and technologies for subsurface hydrological investigations and management requires not only the development of robust technologic tools and approaches, but educational strategies and tools capable of building confidence among diverse users.

54 ENVIRONMENTAL SCIENCES↗

Expendable Launch Vehicle Studies

Analytical support studies of expendable launch vehicles concentrates on the stability of the dynamics during launch especially during or near the region of maximum dynamic pressure. The in-plane linearized dynamic equations of a generic launch vehicle with multiple flexible bending and fuel sloshing modes are developed. The design of a robust LQR controller based on the reduced order system is accomplished using the parameter perturbation technique. The ELV modeling and analysis team has spent the past three years working on the theoretical development and application of sensitivity analysis to solve eigenvalue problems associated with structural dynamics. Specific application areas include stochastic vibrations, viscously damped vibrations and nonlinear dynamics. In stochastic linear vibrations, sensitivity analysis methods have been developed which determined the expectation and variance of the response for given probability density functions of various stochastic parameters. Examples worked include beam vibrations with up to six stochastic parameters; and the eigensensitivity results differ little from those obtained, with much effect, from Monte Carlo techniques. For viscously damped vibrations, eigensensitivity analysis has given excellent results for homogeneous beams, modeled by either quadratic or quartic eigenvalue equations. Specific applications include Kevin and Maxwell-type viscoelastic beams.

Bainum, P. M.↗

Development of a Maxwell X-57 High Lift Motor Reference Design

NASA's all-electric X-57 airplane will utilize 14 electric motors, of which 12 are exclusively for lift augmentation during takeoff and landing. This report covers the design and development process taken to create an open reference model representative of the 12 lift augmenting motors. A combined worst case scenario was used as the design point, which represents the simultaneously occurring worst case aspects of thermal, static stress, electromagnetic, and rotor dynamic conditions. This work also highlights the tightly coupled nature of aerospace electric motor design, requiring constant iteration between all disciplines involved. Further adding to the uniqueness is the cooling method, which is limited to nacelle skin forced convection cooling only, no internal air flow is permitted. The stator outer diameter limit of 156.45 mm greatly impacts the degree of coupling between the electromagnetic design with the thermal analysis. The permanent magnet synchronous motor developed here operates between 385 V and 538 V, at a peak current of 50 A. Detailed electromagnetic, thermal, static load, and rotordynamic analysis was completed for this electric motor; all of which are required for a full design. The rotordynamic analysis took into consideration the motor housing which is designed specifically for this motor. The final electric motor has a mass of 2.34 kg, produces 24.1 Nm of torque with a specific power of 5.56 kW/kg, and has an efficiency of 96.61% at the combined worst case design point.

Electric Aircraft↗

Evaluation of Advanced Stirling Convertor Net Heat Input Correlation Methods Using a Thermal Standard

Test hardware used to validate net heat prediction models. Problem: Net Heat Input cannot be measured directly during operation. Net heat input is a key parameter needed in prediction of efficiency for convertor performance. Efficiency = Electrical Power Output (Measured) divided by Net Heat Input (Calculated). Efficiency is used to compare convertor designs and trade technology advantages for mission planning.

Briggs, Maxwell↗

Steady State Thermal Analyses of SCEPTOR X-57 Wingtip Propulsion

Electric aircraft concepts enable advanced propulsion airframe integration approaches that promise increased efficiency as well as reduced emissions and noise. NASA's fully electric Maxwell X-57, developed under the SCEPTOR program, features distributed propulsion across a high aspect ratio wing. There are 14 propulsors in all: 12 high lift motor that are only active during take off and climb, and 2 larger motors positioned on the wingtips that operate over the entire mission. The power electronics involved in the wingtip propulsion are temperature sensitive and therefore require thermal management. This work focuses on the high and low fidelity heat transfer analysis methods performed to ensure that the wingtip motor inverters do not reach their temperature limits. It also explores different geometry configurations involved in the X-57 development and any thermal concerns. All analyses presented are performed at steady state under stressful operating conditions, therefore predicting temperatures which are considered the worst-case scenario to remain conservative.

Thermal↗

Characterization of a hard X-ray self-seeding diamond crystal orientation

We present a method to accurately control the photon energies for hard X-ray Self-seeding schemes with a single crystal monochromator in transmissive geometry. The energy calibration is performed by measuring which pairs of the machine pitch and yaw angles for different crystallographic planes reflect the X-ray at the same wavelength. The free parameters of an analytical formula for the self-seeding energies are determined by fitting the observed intersections and the normalized derivative with respect to the pitch and yaw angles in the observed intersections. The method requires a hard X-ray spectrometer, but it does not rely on its absolute energy calibration. Instead, identifying the self-seeded energies above the SASE background or the monochromatic notches within the SASE bandwidth is sufficient for the calibration.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Flow Applications of the Least Squares Finite Element Method

The main thrust of the effort has been towards the development, analysis and implementation of the least-squares finite element method (LSFEM) for fluid dynamics and electromagnetics applications. In the past year, there were four major accomplishments: 1) special treatments in computational fluid dynamics and computational electromagnetics, such as upwinding, numerical dissipation, staggered grid, non-equal order elements, operator splitting and preconditioning, edge elements, and vector potential are unnecessary; 2) the analysis of the LSFEM for most partial differential equations can be based on the bounded inverse theorem; 3) the finite difference and finite volume algorithms solve only two Maxwell equations and ignore the divergence equations; and 4) the first numerical simulation of three-dimensional Marangoni-Benard convection was performed using the LSFEM.

Jiang, Bo-Nan↗

Computer Simulation of Microwave Devices

The accurate simulation of cold-test results including dispersion, on-axis beam interaction impedance, and attenuation of a helix traveling-wave tube (TWT) slow-wave circuit using the three-dimensional code MAFIA (Maxwell's Equations Solved by the Finite Integration Algorithm) was demonstrated for the first time. Obtaining these results is a critical step in the design of TWT's. A well-established procedure to acquire these parameters is to actually build and test a model or a scale model of the circuit. However, this procedure is time-consuming and expensive, and it limits freedom to examine new variations to the basic circuit. These limitations make the need for computational methods crucial since they can lower costs, reduce tube development time, and lessen limitations on novel designs. Computer simulation has been used to accurately obtain cold-test parameters for several slow-wave circuits. Although the helix slow-wave circuit remains the mainstay of the TWT industry because of its exceptionally wide bandwidth, until recently it has been impossible to accurately analyze a helical TWT using its exact dimensions because of the complexity of its geometrical structure. A new computer modeling technique developed at the NASA Lewis Research Center overcomes these difficulties. The MAFIA three-dimensional mesh for a C-band helix slow-wave circuit is shown.

Kory, Carol L.↗

On numerical errors to the fields surrounding a relativistically moving particle in PIC codes

The particle-in-cell (PIC) method is widely used to model the self-consistent interaction between discrete particles and electromagnetic fields. It has been successfully applied to problems across plasma physics including plasma based acceleration, inertial confinement fusion, magnetically confined fusion, space physics, astrophysics, high energy density plasmas. In many cases the physics involves how relativistic particles (those with high relativistic γ factors) are generated and interact with plasmas. However, when relativistic particles stream across the grid, both in “vacuum” and in plasma, many numerical issues may arise which can lead to unphysical results. We present a detailed analysis of how discretized Maxwell solvers used in PIC codes can lead to numerical errors to the fields that surround particles that move at relativistic speeds across the grid. Expressions for the axial electric field as integrals in k space are presented that reveal two types of errors. The first arises from errors to the numerator of the integrand and leads to unphysical fields that are antisymmetric about the particle. Furthermore, the second arises from errors to the denominator of the integrand and lead to Cherenkov like radiation in “vacuum”. These fields are not anti-symmetric, extend behind the particle, and cause the particle to accelerate or decelerate depending on the solver and parameters. The unphysical fields are studied in detail for two representative solvers - the Yee solver and the FFT based solver. Although the Cherenkov fields are absent, the space charge fields are still present in the fundamental Brillouin zone for the FFT based solvers. In addition, the Cherenkov fields are present in higher order zones for the FFT based solvers. Comparison between the analytical solutions and PIC simulation results are presented. A solution for eliminating these unphysical fields by modifying the k operator in the axial direction is also presented. Using a customized finite difference solver, this solution was successfully implemented into OSIRIS [1]. Results from the customized solver are also presented. Additionally, this solution will be useful for a beam of particles that all move in one direction with a small angular divergence.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Using SIMS to decode noisy stratigraphic δ 13 C variations in Ediacaran carbonates

Carbonate carbon isotope (δ 13 C carb ) chemostratigraphy is a valuable tool in Precambrian stratigraphic correlation. The effectiveness of this tool rests on the assumption that δ 13 C carb data record global seawater signals. However, in some cases δ 13 C carb data exhibit rapid and noisy stratigraphic variations that may have been influenced by authigenic or diagenetic carbonate minerals. To assess the contribution of non-primary minerals towards bulk carbonate carbon isotope values, we acquired SIMS (secondary ion mass spectrometry) δ 13 C carb data, electron microscopic data, and CL (cathodoluminescence) microscopic data from two Ediacaran successions—the lower Member II of the Doushantuo Formation in South China and the Mooifontein Member of the Zaris Formation in southern Namibia. The Doushantuo samples came from a stratigraphic interval with noisy meter-scale δ 13 C carb variations of up to 10‰, whereas the Mooifontein Member is characterized by consistent δ13Ccarb values with limited meter-scale variations of less than ~4‰. Our data show that the meter-scale stratigraphic variations in δ 13 C carb are also mirrored in the SIMS data at μm-mm scales in both Doushantuo and Mooifontein samples. In the Doushantuo samples, SIMS δ 13 C carb values of authigenic calcite vary by up to 10‰ over μm-mm scales and can be either higher or lower than those of the co-existing dolomite matrix, which also appears to be affected by authigenic carbonate on the basis of petrographic observation. Bulk-sample δ 13 C carb values measured on powders microdrilled from the same SIMS specimens are within the range of SIMS δ 13 C carb values. Thus, we infer that bulk-sample δ 13 C carb values of Doushantuo samples represent mixtures of different carbonate components. The Doushantuo and Mooifontein SIMS δ 13 C carb data validate the traditional method of using chemostratigraphic consistency to evaluate authigenic/diagenetic alteration, but given that secondary phases may be 13 C-enriched relative to matrix carbonate, caution is warranted when taking maximum values to approximate primary chemostratigraphic trends.

58 GEOSCIENCES↗

Kinetic simulations of collision-less plasmas in open magnetic geometries

Laboratory plasmas in open magnetic geometries can be found in many different applications such as (a) scrape-of-layer (SOL) and divertor regions in toroidal confinement fusion devices, (b) linear divertor simulators, (c) plasma-based thrusters and (d) magnetic mirrors etc. A common feature of these plasma systems is the need to resolve, in addition to velocity space, at least one physical dimension (e.g. along flux lines) to capture the relevant physics. In general, this requires a kinetic treatment. Fully kinetic particle-in-cell (PIC) simulations can be applied but at the expense of large computational effort. A common way to resolve this is to use a hybrid approach: kinetic ions and fluid electrons. In the present work, the development of a hybrid PIC computational tool suitable for open magnetic geometries is described which includes (a) the effect of non-uniform magnetic fields, (b) finite fully-absorbing boundaries for the particles and (c) volumetric particle sources. Analytical expressions for the momentum transport in the paraxial limit are presented with their underlying assumptions and are used to validate the results from the PIC simulations. A general method is described to construct discrete particle distribution functions in a state of mirror-equilibrium. This method is used to obtain the initial state for the PIC simulation. Collisionless simulations in a mirror geometry are performed. The results show that the effect of magnetic compression is correctly described and momentum is conserved. The self-consistent electric field is calculated and is shown to modify the ion velocity distribution function in a manner consistent with analytic theory. Based on this analysis, the ion distribution function is understood in terms of a loss-cone distribution and an isotropic Maxwell-Boltzmann distribution driven by a volumetric plasma source. Finally, the inclusion of a Monte Carlo based Fokker-Planck collision operator is discussed in the context of future work.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Integrating a ponderomotive guiding center algorithm into a quasi-static particle-in-cell code based on azimuthal mode decomposition

High fidelity modeling of plasma based acceleration (PBA) requires the use of three dimensional, fully nonlinear, and kinetic descriptions based on the particle-in-cell (PIC) method. In PBA an intense particle beam or laser (driver) propagates through a tenuous plasma whereby it excites a plasma wave wake. Three-dimensional PIC algorithms based on the quasi-static approximation (QSA) have been successfully applied to efficiently model the interaction between relativistic charged particle beams and plasma. In a QSA PIC algorithm, the plasma response to a charged particle beam or laser driver is calculated based on forces from the driver and self-consistent forces from the QSA form of Maxwell's equations. These fields are then used to advance the charged particle beam or laser forward by a large time step. Since the time step is not limited by the regular Courant-Friedrichs-Lewy (CFL) condition that constrains a standard 3D fully electromagnetic PIC code, a 3D QSA PIC code can achieve orders of magnitude speedup in performance. Recently, a new hybrid QSA PIC algorithm that combines another speedup technique known as an azimuthal Fourier decomposition has been proposed and implemented. This hybrid algorithm decomposes the electromagnetic fields, charge and current density into azimuthal harmonics and only the Fourier coefficients need to be updated, which can reduce the algorithmic complexity of a 3D code to that of a 2D code. Modeling the laser-plasma interaction in a full 3D electromagnetic PIC algorithm is very computationally expensive due the enormous disparity of physical scales to be resolved. In the QSA the laser is modeled using the ponderomotive guiding center (PGC) approach. We describe how to implement a PGC algorithm compatible for the QSA PIC algorithms based on the azimuthal mode expansion. Here this algorithm permits time steps orders of magnitude larger than the cell size and it can be asynchronously parallelized. Details on how this is implemented into the QSA PIC code that utilizes an azimuthal mode expansion, QPAD, are also described. Benchmarks and comparisons between a fully 3D explicit PIC code (OSIRIS), as well as a few examples related to laser wakefield acceleration, are presented.

97 MATHEMATICS AND COMPUTING↗

Documentation of the recording environment

The distortion and degradation of the informational content of satellite sensed data by scene generation and atmospheric propagation are discussed. A method is described for obtaining a system point spread function representing the fractional radiance poisoning each pixel gives and receives. The sensor data output is then reprocessed or restored. The degree of restoration is signal to noise dependent. A technique is given for approximating the impact of the atmosphere on the propagation of scene radiance.

Ritter, M.↗

A global time-dependent model of thunderstorm electricity. I - Mathematical properties of the physical and numerical models

A time-dependent model is introduced that can be used to simulate the interaction of a thunderstorm with its global electrical environment. The model solves the continuity equation of the Maxwell current, which is assumed to be composed of the conduction, displacement, and source currents. Boundary conditions which can be used in conjunction with the continuity equation to form a well-posed initial-boundary value problem are determined. Properties of various components of solutions of the initial-boundary value problem are analytically determined. The results indicate that the problem has two time scales, one determined by the background electrical conductivity and the other by the time variation of the source function. A numerical method for obtaining quantitative results is introduced, and its properties are studied. Some simulation results on the evolution of the displacement and conduction currents during the electrification of a storm are presented.

Browning, G. L.↗

A cradle-to-grave life-cycle-assessment of dry-processed Li-ion batteries for electric vehicles

Maxwell-type dry processing has emerged as a promising manufacturing technology for high-areal-loading Li-ion battery electrodes, offering a significant advantage by eliminating the use of the toxic and costly NMP solvent. This study developed a cradle-to-grave life-cycle-assessment model to evaluate the environmental impacts of this innovative technology, benchmarking the results against the conventional NMP-based technique for a configured 42 kWh NMC622-graphite battery pack. Furthermore, the findings reveal that the dry-processing method merits a 4.8% lower energy consumption and achieves environmental impacts up to 47.5% lower in 12 out of 13 categories, highlighting its environmental benefits compared with the NMP-based manufacturing processes.

Electrode↗

Thermal fluctuations in the Landau-Lifshitz-Bloch model

A formulation for thermal noise in the stochastic form of the Landau-Lifshitz-Bloch equation used for modeling the magnetization dynamics at elevated temperatures is presented here. The diffusion coefficients for thermal fluctuations are obtained via the Fokker-Plank equation using the mean-field approximation of the field in defining the free energy. The presented model leads to a mean magnetization consistent with the equilibrium magnetization for small and large particles. The distribution of the magnetization magnitude is of the Maxwell-Boltzmann type. The presented model was tested by studying the equilibrium magnetization in macrospin particles at high temperatures. The model is appealing for multiscale modeling, such as modeling heat assisted magnetic recording systems and all-optical magnetization reversal.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Maximum likelihood estimation of difference scaling functions for suprathreshold judgments

Maximum likelihood estimation (MLE) has been used to produce perceptual scales from binary judgments of triads and quadruples. This method relies on Thurstone’s theory of a stochastic perceptual process where the perceived difference of two stimuli is the difference in their perceived strengths. It is possible that the perception of a suprathreshold difference is overestimated when adding smaller differences, a phenomenon referred to as diminishing returns. The current approach to construct a perceptual scale using MLE does not account for this phenomenon. We present a way to model the perception of differences using MLE and Thurstone’s theory, adapted to allow the possibility of diminishing returns. This method is validated using Monte Carlo simulated responses to experimental triads and can correctly model diminishing returns, the absence of diminishing returns, and the opposite of diminishing returns both in the cases when a perceptual scale is known and when the true perceived strengths of the stimuli are unknown. Additionally, this method was applied to empirical data sets to determine its feasibility in investigations of perception. Ultimately, it was found that this analysis allows for more accurate modeling of suprathreshold difference judgments, a more complete understanding of the perceptual processes underlying comparisons, and the evaluation of Thurstone’s theory of difference judgments.

97 MATHEMATICS AND COMPUTING↗

Method of fabricating a three-dimensionally patterned mechanical energy absorptive material

A three-dimensionally patterned energy absorptive material and fabrication method having multiple layers of patterned filaments extrusion-formed from a curable pre-cursor material and stacked and cured in a three-dimensionally patterned architecture so that the energy absorptive material produced thereby has an engineered bulk property associated with the three-dimensionally patterned architecture.

Duoss, Eric B.↗