Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Differential Geometry”

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 847 records · Page 47
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

Performance evaluation of the JPL TOPSAR system: An across-track interferometric SAR system for topographic mapping

An across-track interferometric (XTI) SAR has the capability of determining the three dimensional target location from the recorded radar echos thus providing high spatial resolution both horizontally and vertically. The output of an XTI SAR will typically include both a radar backscatter map and a digital elevation model (DEM). The accuracy of the derived DEM is determined by factors such as data acquisition geometry, signal-to-noise ratio, knowledge of platform position and attitude as well as the accuracy of the processing system. In the summer of 1992 the authors acquired TOPSAR data over the Ft. Irwin area in California, a desert area with significant relief (height standard deviation 150 m). Very accurate DEM's have been derived for this area by the Topographic Engineering Center (TEC) using digital correlation methods on 1:20,000 scale digitized photographs. Corner reflectors were deployed in the area, and their locations were determined to cm accuracies by the Defense Mapping Agency (DMA) using differential GPS techniques. DEM's generated from the acquired radar data have been rotated and translated to overlay the reference DEM's provided by TEC. A detailed description of the errors and their characteristics will be given. The standard deviations of the radar data sets we acquired measured over a 5.6 by 7 km area were 1.9 and 2.3 m. The corresponding numbers for flat areas were 1.1 and 2.0 m and for mountain areas 3.3 and 2.2 m.

Madsen, Soren N.↗
DOE OSTI2013

A GENERALIZED STOCHASTIC COLLOCATION APPROACH TO CONSTRAINED OPTIMIZATION FOR RANDOM DATA IDENTIFICATION PROBLEMS

We present a scalable, parallel mechanism for stochastic identification/control for problems constrained by partial differential equations with random input data. Several identification objectives will be discussed that either minimize the expectation of a tracking cost functional or minimize the difference of desired statistical quantities in the appropriate $L^p$ norm, and the distributed parameters/control can both deterministic or stochastic. Given an objective we prove the existence of an optimal solution, establish the validity of the Lagrange multiplier rule and obtain a stochastic optimality system of equations. The modeling process may describe the solution in terms of high dimensional spaces, particularly in the case when the input data (coefficients, forcing terms, boundary conditions, geometry, etc) are affected by a large amount of uncertainty. For higher accuracy, the computer simulation must increase the number of random variables (dimensions), and expend more effort approximating the quantity of interest in each individual dimension. Hence, we introduce a novel stochastic parameter identification algorithm that integrates an adjoint-based deterministic algorithm with the sparse grid stochastic collocation FEM approach. This allows for decoupled, moderately high dimensional, parameterized computations of the stochastic optimality system, where at each collocation point, deterministic analysis and techniques can be utilized. The advantage of our approach is that it allows for the optimal identification of statistical moments (mean value, variance, covariance, etc.) or even the whole probability distribution of the input random fields, given the probability distribution of some responses of the system (quantities of physical interest). Our rigorously derived error estimates, for the fully discrete problems, will be described and used to compare the efficiency of the method with several other techniques. Numerical examples illustrate the theoretical results and demonstrate the distinctions between the various stochastic identification objectives.

Webster, Clayton G↗
DOE OSTI2013

Quantum-mechanical predictions of electron-induced ionization cross sections of DNA components

Ionization of biomolecules remains still today rarely investigated on both the experimental and the theoretical sides. In this context, the present work appears as one of the first quantum mechanical approaches providing a multi-differential description of the electron-induced ionization process of the main DNA components for impact energies ranging from the target ionization threshold up to about 10 keV. The cross section calculations are here performed within the 1st Born approximation framework in which the ejected electron is described by a Coulomb wave whereas the incident and the scattered electrons are both described by a plane wave. The biological targets of interest, namely, the DNA nucleobases and the sugar-phosphate backbone, are here described by means of the GAUSSIAN 09 system using the restricted Hartree-Fock method with geometry optimization. The theoretical predictions also obtained have shown a reasonable agreement with the experimental total ionization cross sections while huge discrepancies have been pointed out with existing theoretical models, mainly developed within a semi-classical framework.

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

Computational Modeling of Graphite Degradation due to Molten Salt Infiltration and Wear

Molten-salt reactors (MSRs) represent a promising next-generation reactor design, with graphite serving as a moderator and/or reflector in several designs. However, due to limited experimental data and operational experience, a technical understanding of the structural integrity of graphite in molten salt environments remains incomplete. This report presents a modeling-based evaluation of graphite degradation in MSR environments, focusing on the effects of salt infiltration in fuel salt-based designs and surface wear in pebble bed reactor designs. The objective of this study is to enhance understanding of the structural integrity challenges posed by these degradation mechanisms and to provide a framework for assessing graphite behavior in MSRs. The first part of the report investigates the phenomenon of molten salt infiltration into graphite. This infiltration occurs when molten salt permeates the interconnected pore structure of the graphite moderator, driven by factors such as pressure differentials and the physical properties of both the salt and graphite. The infiltration process is influenced by characteristics of the pore structure, viscosity of the molten salt, and the interfacial energies between the graphite, salt, and the atmosphere within the graphite pore. Utilizing a coupled multiphysics modeling approach with Grizzly software, the study evaluates the stress induced by internal heat sources due to infiltration, which can lead to structural concerns. This evaluation is crucial for understanding how infiltration affects the mechanical integrity of graphite components in MSRs. The study considers the Molten-Salt Reactor Experiment (MSRE) graphite stringer geometry due to the availability of relevant data. Through detailed finite element analysis, the study examines stress distributions at varying infiltration percentages, revealing that stress levels increase with higher amounts of infiltration. Rare-event simulations, using the parallel subset simulation (PSS) framework, further quantify the failure probabilities under input uncertainties, with a user-specified failure metric. The PSS framework also identifies critical input parameters that significantly affect the stress values, including infiltration amount, thermal conductivity, and power density. Additionally, considering realistic reactor scenarios, the analysis was performed to account for the combined effects of radiation and infiltration, and modeling strategies on how to analyze new reactor designs or new graphite grades are discussed. The second part of the report focuses on wear mechanisms in pebble bed-based MSRs. As graphite fuel pebbles interact with the graphite reflector block, wear can result in material loss and the formation of surface defects, which may act as stress concentrators. A similar multiphysics modeling framework is employed to assess the impact of wear on the structural integrity of graphite components. This study considers a generic fluoride-cooled high-temperature reactor (gFHR) design due to the availability of comprehensive data. Worst-case scenario dimensions of the reflector blocks were analyzed under thermal and radiation conditions. Subsequently, wear in the form of idealized pits and grooves is modeled on the inner surface of the graphite block, with the maximum stress from previous simulations. The simulations show that groove-type defects are more detrimental than pits, leading to higher stress concentrations. Considering worst-case simulation scenarios and experimental wear rates, it was determined that the formation of a surface defect critical enough to affect the stress may not be possible in a gFHR design. Overall, the findings of this research contribute to the development of robust modeling tools for predicting graphite behavior under various operational conditions in MSRs.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS1992

Evaluation and application of the Baldwin-Lomax turbulence model in two-dimensional, unsteady, compressible boundary layers with and without separation in engine inlets

There is a practical need to model high speed flows that exist in jet engine inlets. The boundary layers that form in these inlets may be turbulent or laminar and either separated or attached. Also, unsteady supersonic inlets may be subject to frequent changes in operating conditions. Some changes in the operating conditions of the inlets may include varying the inlet geometry, bleeds and bypasses, and rotating or translating the centerbody. In addition, the inlet may be either started or unstarted. Therefore, a CFD code, used to model these inlets, may have to run for several different cases. Also, since the flow conditions through an unsteady inlet may be continually fluctuating, the CFD code which models these flows may have to be run over many time steps. Therefore, it would be beneficial that the code run quickly. Many turbulence models, however, are cumbersome to implement and require a lot of computer time to run, since they add to the number of differential equations to be solved to model a flow. The Baldwin-Lomax turbulence model is a popular model. It is an algebraic, eddy viscosity model. The Baldwin-Lomax model is used in many CFD codes because it is quick and easy to implement. In this paper, we will discuss implementing the Baldwin-Lomax turbulence model for both steady and unsteady compressible flows. In addition, these flows may be either separated or attached. In order to apply this turbulence model to flows which may be subjected to these conditions, certain modifications should be made to the original Baldwin-Lomax model. We will discuss these modifications and determine whether the Baldwin-Lomax model is a viable turbulence model that produces reasonably accurate results for high speed flows that can be found in engine inlets.

Sakowski, Barbara↗
NASA NTRS1993

Two wavelength satellite laser ranging using SPAD

When ranging to satellites with lasers, there are several principal contributions to the error budget: from the laser ranging system on the ground, from the satellite retroarray geometry, and from the atmosphere. Using a single wavelength, we have routinely achieved a ranging precision of 8 millimeters when ranging to the ERS-1 and Starlette satellites. The systematic error of the atmosphere, assuming the existing dispersion models, is expected to be of the order of 1 cm. Multiple wavelengths ranging might contribute to the refinement of the existing models. Taking into account the energy balance, the existing picosecond lasers and the existing receiver and detection technology, several pairs or multiple wavelengths may be considered. To be able to improve the atmospheric models to the subcentimeter accuracy level, the differential time interval (DTI) has to be determined within a few picoseconds depending on the selected wavelength pair. There exist several projects based on picosecond lasers as transmitters and on two types of detection techniques: one is based on photodetectors, like photomultipliers or photodiodes connected to the time interval meters. Another technique is based on the use of a streak camera as an echo signal detector, temporal analyzer, and time interval vernier. The temporal analysis at a single wavelength using the streak camera showed the complexity of the problem.

Prochazka, Ivan↗
DOE OSTI2015

Preparation and characterization of novel foamed porous glass-ceramics

Foamed glass-ceramics without using foaming agent have been synthesized in a novel glass system of SrO-CaO-Al{sub 2}O{sub 3}-TiO{sub 2}-B{sub 2}O{sub 3}-SiO{sub 2}-P{sub 2}O{sub 5}-M{sub x}O{sub y} (where M = Ba, Mg, La, Ce and Ni) by a simple process of powder sintering. The glass and glass-ceramics are characterized by dilatometry, differential scanning calorimetry (DSC), heating stage microscopy (HSM), X-ray diffraction (XRD), field emission scanning electron microscopy (FESEM), optical microscopy and Fourier transformed infrared spectroscopy (FTIR). All the glasses formed are amorphous and the glass transition temperature and dilatometric softening temperature of these glasses are found to be in the range 673–678 °C and 706–728 °C respectively. The glasses are highly stable as indicated by the DSC evaluated glass stability parameters of the range 195–240 °C. Quantitative sintering study of glass powder compacts revealed swelling in the samples with NiO and CeO{sub 2} corresponding to a geometry change of 75 and 108% around 900 °C respectively. With reference to this finding the glass powder compacts are heated to 900 °C and the foamed glass-ceramics are obtained. Characteristic crystalline silicate phases have been identified in the XRD studies and their microstructures are recorded by FESEM. Optical microscope study of the foamed samples revealed formation of bigger foamed cavity with residual pores in samples with NiO and CeO{sub 2} in comparison to samples with BaO, MgO and La{sub 2}O{sub 3}. The mean pore diameters of the samples with NiO and CeO{sub 2} are determined to be 43 and 32 μm, and their respective porosities are 2.34 and 1.82 cm{sup 3}/g respectively. Thus NiO and CeO{sub 2} are found to be very effective to obtain foamed glass-ceramics without using foaming agent by the viscous flow sintering of fine glass powder compacts along with the reduction of the respective polyvalent ions. - Highlights: • Synthesis of foamed porous glass-ceramics without foaming agent by sintering method • Only powder compact yielded foamed porous glass-ceramics but bulk glass did not. • Glasses containing NiO and CeO{sub 2} exhibited significant foaming efficiency. • Bloating of entrapped gas during viscous flow sintering is the origin of foaming. • Residual void created pores in the sintered glass-ceramics as evidenced in FESEM.

36 MATERIALS SCIENCE↗
DOE OSTI2014

Nonlinear physics of shear Alfvén waves

Shear Alfvén waves (SAW) play fundamental roles in thermonuclear plasmas of fusion interest, since they are readily excited by energetic particles in the MeV range as well as by the thermal plasma components. Thus, understanding fluctuation induced transport in burning plasmas requires understanding nonlinear SAW physics. There exist two possible routes to nonlinear SAW physics: (i) wave-wave interactions and the resultant spectral energy transfer; (ii) nonlinear wave-particle interactions of SAW instabilities with energetic particles. Within the first route, it is advantageous to understand and describe nonlinear processes in term of proximity of the system to the Alfvénic state, where wave-wave interactions are minimized due to the cancellation of Reynolds and Maxwell stresses. Here, various wave-wave nonlinear dynamics are elucidated in terms of how they break the Alfvénic state. In particular, we discuss the qualitative and quantitative modification of the SAW parametric decay process due to finite ion compressibility and finite ion Larmor radius. We also show that toroidal geometry plays a crucial role in the nonlinear excitation of zonal structures by Alfvén eigenmodes. Within the second route, the coherent nonlinear dynamics of structures in the energetic particle phase space, by which secular resonant particle transport can occur on meso- and macro-scales, must be addressed and understood. These 'nonlinear equilibria' or 'phase-space zonal structures' dynamically evolve on characteristic (fluctuation induced) turbulent transport time scales, which are generally of the same order of the nonlinear time scale of the underlying fluctuations. In this work, we introduce the general structure of nonlinear Schrödinger equations with complex integro-differential nonlinear terms, which govern these physical processes. To elucidate all these aspects, theoretical analyses are presented together with numerical simulation results.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1997

Constrained Multipoint Aerodynamic Shape Optimization Using an Adjoint Formulation and Parallel Computers

An aerodynamic shape optimization method that treats the design of complex aircraft configurations subject to high fidelity computational fluid dynamics (CFD), geometric constraints and multiple design points is described. The design process will be greatly accelerated through the use of both control theory and distributed memory computer architectures. Control theory is employed to derive the adjoint differential equations whose solution allows for the evaluation of design gradient information at a fraction of the computational cost required by previous design methods. The resulting problem is implemented on parallel distributed memory architectures using a domain decomposition approach, an optimized communication schedule, and the MPI (Message Passing Interface) standard for portability and efficiency. The final result achieves very rapid aerodynamic design based on a higher order CFD method. In order to facilitate the integration of these high fidelity CFD approaches into future multi-disciplinary optimization (NW) applications, new methods must be developed which are capable of simultaneously addressing complex geometries, multiple objective functions, and geometric design constraints. In our earlier studies, we coupled the adjoint based design formulations with unconstrained optimization algorithms and showed that the approach was effective for the aerodynamic design of airfoils, wings, wing-bodies, and complex aircraft configurations. In many of the results presented in these earlier works, geometric constraints were satisfied either by a projection into feasible space or by posing the design space parameterization such that it automatically satisfied constraints. Furthermore, with the exception of reference 9 where the second author initially explored the use of multipoint design in conjunction with adjoint formulations, our earlier works have focused on single point design efforts. Here we demonstrate that the same methodology may be extended to treat complete configuration designs subject to multiple design points and geometric constraints. Examples are presented for both transonic and supersonic configurations ranging from wing alone designs to complex configuration designs involving wing, fuselage, nacelles and pylons.

Reuther, James↗
DOE OSTI2020

VALVE TIMING AND AREA PROFILE SELECTION FOR HYDRAULIC PUMPS AND MOTORS

A significant amount of research has been conducted to select valve timing and area profiles that create efficient and quiet hydraulic pumps and motors. Numerous active valve architectures have been modeled and optimized, but the rationale for the final solution is often unclear. The solution is usually highly dependent on the modeled valve geometry constraints and the duty cycle of the pump or motor for which the valve was optimized. This paper first presents an ideal valve area profile and what compromises are made to design for a working pressure and speed range. The ideal valve area profile is resolved using a piston-cylinder model with a pressure dependent bulk modulus and variable valve openings. The piston-cylinder is first modeled as a pump with passive (check) valves and a specified pressure drop across the valves to determine the ideal area profile. A generalization of the ideal timing is demonstrated for positive pressure differential and positive rotation direction, also known as the first quadrant. The effects on valve timing when varying pressure are shown for the first quadrant, but the trends can be extrapolated for all quadrants of operation. The ideal valve profile is implemented as fixed valve timing, meaning the valve area is only a function of the timing angle, in a pump-motor. The pump-motor is simulated in one rotation direction through a pressure range. Performance is high in pumping operation, but when the pressure differential is reversed, cylinder pressure spikes ensue. Two strategies to modify an ideal valve area profile are presented: timing grooves and a pressure shifted valve timing. Timing grooves reduce pressure spikes and cavitation in the cylinder but generally increase throttling losses. A pressure shifted valve timing has lower throttling energy losses, making it the favored solution.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2025

Cartesian equivariant representations for learning and understanding molecular orbitals

Qualitative and quantitative orbital properties such as bonding/antibonding character, localization, and orbital energies are critical to how chemists understand reactivity, catalysis, and excited-state behavior. Despite this, representations of orbitals in deep learning models have been very underdeveloped relative to representations of molecular geometries and Hamiltonians. Here, we apply state-of-the-art equivariant deep learning architectures to the task of assigning global labels to orbitals, namely energies characterizations, given the molecular coefficients from Hartree–Fock or density functional theory. The architecture we have developed, the Cartesian Equivariant Orbital Network (CEONET), shows how molecular orbital coefficients are readily featurized as equivariant node features common to all graph-based machine-learned potentials. We find that CEONET performs well at predicting difficult quantitative labels such as the orbital energy and orbital entropy. Furthermore, we find that the CEONET representation provides an intuitive latent space for differentiating orbital character for the qualitative assignment of e.g. bonding or antibonding character. In addition to providing a useful representation for further integrating deep learning with electronic structure theory, we expect CEONET to be useful for automatizing and interpreting the results of advanced electronic structure methods such as complete active space self-consistent field theory. In particular, the ability of CEONET to infer multireference character via the orbital entropy paves the way toward the machine-learned selection of active spaces.

chemical reactions↗
DOE OSTI2024

AMReX v2024

The software framework, AMReX, supports the development of block-structured adaptive mesh refinement (AMR) algorithms for solving systems of partial differential equations. AMR reduces the computational cost and memory footprint compared to a uniform mesh while preserving the essential local descriptions of different physical processes in complex multiphysics algorithms. AMR uses a hierarchical representation of the solution at multiple levels of resolution where the solution on each level is defined on the union of data containers at that resolution. These data containers, which represent the solution over a logically rectangular subregion of the domain, can contain field data defined on a mesh, Lagrangian particles or combinations of both. In addition to these basic data types, AMReX supports a multilevel embedded boundary representation of complex geometry; linear solvers for cell-centered and nodal data; asynchronous I/O in a native format readable by ParaView, VisIt and yt; and interfaces to hypre and PETSc solvers. AMReX enables applications to run on distributed memory architectures with multicore CPUs and with GPU accelerators. AMReX uses a lightweight abstraction layer that effectively hides the details of the architecture from the application. The framework currently supports CUDA, HIP and SYCL for GPU acceleration and OpenMP for multi-core CPU architectures.

Almgren, Ann↗
DOE OSTI2016

A pencil beam dose calculation model for CyberKnife system

Purpose: CyberKnife system is initially equipped with fixed circular cones for stereotactic radiosurgery. Two dose calculation algorithms, Ray-Tracing and Monte Carlo, are available in the supplied treatment planning system. A multileaf collimator system was recently introduced in the latest generation of system, capable of arbitrarily shaped treatment field. The purpose of this study is to develop a model based dose calculation algorithm to better handle the lateral scatter in an irregularly shaped small field for the CyberKnife system. Methods: A pencil beam dose calculation algorithm widely used in linac based treatment planning system was modified. The kernel parameters and intensity profile were systematically determined by fitting to the commissioning data. The model was tuned using only a subset of measured data (4 out of 12 cones) and applied to all fixed circular cones for evaluation. The root mean square (RMS) of the difference between the measured and calculated tissue-phantom-ratios (TPRs) and off-center-ratio (OCR) was compared. Three cone size correction techniques were developed to better fit the OCRs at the penumbra region, which are further evaluated by the output factors (OFs). The pencil beam model was further validated against measurement data on the variable dodecagon-shaped Iris collimators and a half-beam blocked field. Comparison with Ray-Tracing and Monte Carlo methods was also performed on a lung SBRT case. Results: The RMS between the measured and calculated TPRs is 0.7% averaged for all cones, with the descending region at 0.5%. The RMSs of OCR at infield and outfield regions are both at 0.5%. The distance to agreement (DTA) at the OCR penumbra region is 0.2 mm. All three cone size correction models achieve the same improvement in OCR agreement, with the effective source shift model (SSM) preferred, due to their ability to predict more accurately the OF variations with the source to axis distance (SAD). In noncircular field validation, the pencil beam calculated results agreed well with the film measurement of both Iris collimators and the half-beam blocked field, fared much better than the Ray-Tracing calculation. Conclusions: The authors have developed a pencil beam dose calculation model for the CyberKnife system. The dose calculation accuracy is better than the standard linac based system because the model parameters were specifically tuned to the CyberKnife system and geometry correction factors. The model handles better the lateral scatter and has the potential to be used for the irregularly shaped fields. Comprehensive validations on MLC equipped system are necessary for its clinical implementation. It is reasonably fast enough to be used during plan optimization.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2018

Continuous Particulate Filter State of Health Monitoring Using Radio Frequency Sensing

Reliable means for on-board detection of particulate filter failures or malfunctions are needed to meet diagnostics (OBD) requirements. Detecting these failures, which result in tailpipe particulate matter (PM) emissions exceeding the OBD limit, over all operating conditions is challenging. Current approaches employ differential pressure sensors and downstream PM sensors, in combination with particulate filter and engine-out soot models. These conventional monitors typically operate over narrowly-defined time windows and do not provide a direct measure of the filter’s state of health. In contrast, radio frequency (RF) sensors, which transmit a wireless signal through the filter substrate provide a direct means for interrogating the condition of the filter itself. Here, this study investigated the use of RF sensors for the continuous measurement of filter trapping efficiency, which was compared to downstream measurements with an AVL Microsoot Sensor, and a PM sampling probe simulating the geometry and installation configuration of a conventional PM sensor. The study included several particulate filter failure modes, both above and below the OBD threshold. Finally, the results confirmed the use of RF sensors to provide a direct and continuous measure of the particulate filter’s state of health over a range of typical in-use operating conditions, thereby significantly increasing the time window over which filter failures may be detected.

42 ENGINEERING↗
NASA NTRS1983

Destabilization of a 650 km chemical boundary layer and its bearing on the evolution of the continental crust

A mechanism for the production of a chemical change in the mantle, from primordial silicate compositions above the 650 km discontinuity to differentiated compositions below, is reviewed. Some consequences of this are the stabilization of two layer convection with a temperature contrast between the anhydrous mantle solidus and the geotherm which, at 650 km depth, is lower than any other location in the mantle. With thermal contributions from the concentration of the heat producing elements U, Th and K below the 650 km mantle boundary layer and the higher geotherms in the past, widespread or catastrophic melting may have taken place at this location. An episodic breach of this boundary layer by extensive heat and mass transport may have periodically destroyed any simple two-layer convection geometry in the mantle. Such episodic injections of mass and energy into the upper mantle from below may have been the mechanism responsible for episodes of enhanced surface tectonism and thermal activity which appear to be recorded in apparent polar wandering paths and radiometric ages of continental rocks.

Herzberg, C. T.↗
DOE OSTI2021

Fracture Network Localization Preceding Catastrophic Failure in Triaxial Compression Experiments on Rocks

We quantify the spatial distribution of fracture networks throughout six in situ X-ray tomography triaxial compression experiments on crystalline rocks at confining stresses of 5–35 MPa in order to quantify how fracture development controls the final macroscopic failure of the rock, a process analogous to those that control geohazards such as earthquakes and landslides. Tracking the proportion of the cumulative volume of fractures with volumes >90th percentile to the total fracture volume, ∑ v 90 / v t o t indicates that the fracture networks tend to increase in localization toward these largest fractures for up to 80% of the applied differential stress. The evolution of this metric also matches the evolution of the Gini coefficient, which measures the deviation of a population from uniformity. These results are consistent with observations of localizing low magnitude seismicity before large earthquakes in southern California. In both this analysis and the present work, phases of delocalization interrupt the general increase in localization preceding catastrophic failure, indicating that delocalization does not necessarily indicate a reduction of seismic hazard. However, the proportion of the maximum fracture volume to the total fracture volume does not increase monotonically. Experiments with higher confining stress tend to experience greater localization. To further quantify localization, we compare the geometry of the largest fractures, with volumes >90th percentile, to the best fit plane through these fractures immediately preceding failure. The r 2 scores and the mean distance of the fractures to the plane indicate greater localization in monzonite than in granite. The smaller mean mineral diameter and lower confining stress in the granite experiments may contribute to this result. Tracking these various metrics of localization reveals a close association between macroscopic yielding and the acceleration of fracture network localization. Near yielding, ∑ v 90 / v t o t and the Gini coefficient increase while the mean distance to the final failure plane decreases. Macroscopic yielding thus occurs when the rate of fracture network localization increases.

58 GEOSCIENCES↗
NASA NTRS1988

Nonconforming mortar element methods: Application to spectral discretizations

Spectral element methods are p-type weighted residual techniques for partial differential equations that combine the generality of finite element methods with the accuracy of spectral methods. Presented here is a new nonconforming discretization which greatly improves the flexibility of the spectral element approach as regards automatic mesh generation and non-propagating local mesh refinement. The method is based on the introduction of an auxiliary mortar trace space, and constitutes a new approach to discretization-driven domain decomposition characterized by a clean decoupling of the local, structure-preserving residual evaluations and the transmission of boundary and continuity conditions. The flexibility of the mortar method is illustrated by several nonconforming adaptive Navier-Stokes calculations in complex geometry.

Maday, Yvon↗
NASA NTRS2018

Simulation of Charge Collection to Spacecraft Surfaces: Freja Satellite

As spacecraft travel through plasma, spacecraft surfaces become charged by the collection of charged particles. This process is referred to as Surface Charging. These charges can be detrimental to the vehicle's electronic subsystems as they present a threat of electrostatic discharge (ESD) to onboard circuitry. The process of Surface Charging is complex and is affected by many elements. The charging of each surface is unique. The potential of an individual surface is dependent upon many variables including but not limited to the surface's geometry, material and its location. Each surface also has unique interactions with the surrounding plasma. Other factors that play large roles in the charging process is the density and temperature of plasma ions and electrons. Using Nascap-2k, a model of the Freja satellite has been constructed, and its auroral plasma environment has been imitated to simulate surface charging characteristics. The charging process of the Freja satellite has been modeled iteratively with incremental changes in both the Maxwellian electron temperature (eV) as well as the Gaussian electron energy (eV). This study provides an analysis of the sensitivity between spacecraft surface charging and these two primary variables of electron differential flux.

Wolfley, Olivia H.↗