Starspot evolution, differential rotation, and correlation between chromospheric and photospheric activities on Kepler-411
Explore the source record for details and available documents.
Fukun Xu↗
SEARCH · Search NASA
Results for “differential evolution”
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.
LA17-V-InverseRadTransportProblem-PD2Kb. Final Project Report
In the inverse radiation transport problem, a radiation signal emanating from an unknown object is used to reconstruct geometric and material parameters of the object , e.g., the mass of special nuclear material (SNM), the isotopic composition of the SNM (e.g., the enrichment of 235U), the thickness of shielding, etc. The inverse radiation transport problem is one of the most fundamental and important in the emergency response mission, and has been addressed in previous NA-22 and NA-42/84 investments by using derivative-free optimization methods. Our project was aimed at (i) consolidating previous work by applying polynomial chaos expansion (PCE) for reducing significantly the number of forward model simulations required by derivative-free optimization methods; (ii) developing a general-purpose framework and stand-alone software for solving the inverse radiation transport problem by applying a new “predictive modeling” methodology (predictive modeling of coupled multi-physics systems, PM_CMPS) which enables the fusing of experimental and computational information; and (iii) developing several benchmarks (for unscattered gamma rays) for demonstrating the quantification of non-Gaussian features (e.g., skewness and kurtosis, which characterize asymmetries and long tails) of response distributions by deterministically and efficiently propagating parameter uncertainties (covariances and higher-order moments, where available) using high-order sensitivities computed by the recently developed comprehensive adjoint sensitivity analysis methodology (C-ASAM). With regard to (i), application of PCE, the most significant accomplishment of this project was an implementation of PCE into INVERSE that significantly reduced computation time for the DiffeRential Evolution Adaptive Metropolis (DREAM) and Differential Evolution (DE) methods. For inverse transport problems with measurements of neutron leakage, neutron multiplication, and neutron-induced gamma-rays, the solution times for the DREAM and DE methods were reduced by factors of 30-300, reducing run times from hours to minutes while maintaining the same accuracy as the standard (non-PCE) DREAM and DE. Speedups were further improved by exploring sparse grid quadrature sets using the Toolkit for Adaptive Stochastic Modeling and Non-Intrusive ApproximatioN (Tasmanian) toolset. An initial exploration of sparse grids indicated that further speedups factors of 5-10 are obtainable. The implementation of PCE in INVERSE is not in a form that users can easily access. With regard to (ii), solving the inverse problem using a predictive modeling paradigm, the most significant accomplishment of this project were the derivation and implementation of adjoint-based second derivatives in neutron and photon transport theory. There are two codes that compute these derivatives. SENSPG computes them for uncollided gamma rays and SENSMG computes them for neutron leakage. SENSMG is open-source code available on github. Other significant accomplishments included a stand-alone code, MULTI-PRED, that implements the predictive modeling paradigm, along with a detailed manual. The code was implemented in the existing INVERSE code for inverse problems, but a head-to-head comparison of inverse methods was not achieved. We also derived and implemented the capability to compute sensitivities for multiplicity measurements into the SENSMG code, and this capability was used in a data assimilation application, but it was not used for solving inverse problems. With regard to (iii), developing benchmarks to demonstrate the importance of high-order sensitivity and uncertainty analysis, the most significant accomplishment of this project were two test problems that showed how important second-order sensitivities are.
07 ISOTOPE AND RADIATION SOURCES↗
Parallel Variable Population Multi-Objective Optimizer (pvpmoo) v1.0
This is a parallel variable population multi-objective optimizer with an adaptive unified differential evolution algorithm or a genetic algorithm. It can also be used for single objective optimization. Some features of this code include: 1) The population size varies from generation to generation to save the total # of objective function evaluations. 2) The population is uniformly distributed to a number of parallel processors for simultaneous objective function evaluation. 3) The objective function evaluation can be attained from an external simulation program with control variables in its input file and objectives calculated from its output files. 4) The optimizer includes an adaptive unified differential evolution algorithm and a real value genetic algorithm. The parameters in the unified differential evolution algorithm can be chosen to attain any mutation schemes in the published literature.
Qiang, Ji↗
A parallel variable population multi-objective optimizer for accelerator beam dynamics optimization
The simultaneous optimization of multiple objective functions is needed in many particle accelerator applications. In this paper, we present a parallel evolution based multi-objective optimizer that uses a variable population from generation to generation and an external storage to save good solutions. Two heuristic optimization methods, one uses the unified differential evolution and the other uses the real-coded genetic algorithm, are included in the optimizer to generate next generation candidate solutions, and are compared in the test examples. Finally, as an application, we applied this optimizer to the beam dynamics design optimization of a photoinjector and attained the optimal front solutions after 200 generations with the unified differential evolution offspring production scheme.
46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
Application of Markov Chain Monte Carlo Methods for Uncertainty Quantification in Inverse Transport Problems
Determination of the components of a radioactive source/shield system using the system’s radiation signature is of great importance in homeland security, material safeguards, and waste management. Although significant progress has been made toward solving this inverse transport problem in recent years, work remains to be done to quantify the uncertainty in reconstructed results. In this article we apply two Markov chain Monte Carlo (MCMC) approaches, the delayed rejection adaptive metropolis (DRAM) and differential evolution adaptive metropolis (DREAM) methods, to solve inverse problems and quantify uncertainty. The DRAM method uses delayed rejection combined with global adaptation of the proposal covariance matrix. Furthermore, the DREAM method hybridizes MCMC sampling with the differential evolution (DE) algorithm. In numerical test cases, the DRAM and DREAM methods are shown to be superior to a first-order inverse Hessian approach for problems with noisy data and multiple unknown quantities, with DREAM converging to the posterior distribution more quickly than DRAM. The DREAM and DRAM results indicate that a full posterior distribution is required to quantify uncertainty in many inverse transport problems.
98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗
Temperature-dependent mechanical properties and crystal plasticity parameters for additively manufactured Haynes-214 alloy: Experiments and numerical modeling
Our experimental mechanical testing data demonstrated that the additively manufactured (AM) laser powder bed fusion (L-PBF) Haynes-214 alloy exhibits non-linear mechanical properties as the temperature rises from ambient to 870 °C. Crystal plasticity (CP) simulations provide an effective approach to gaining deeper insights into microstructure-property linkages under thermomechanical loading. This method can reduce the need for costly high-temperature mechanical testing while accounting for the effects of crystallographic texture and grain morphology on the mechanical behavior of AM materials. However, calibrating a CP model is time-consuming because individual simulations are computationally expensive and hundreds (or more) of iterations over parameter sets may be required. To address this issue, we have designed a machine learning-differential evolution (ML-DE) CP framework that can accurately interpolate the tensile properties of AM L-PBF Haynes-214 alloy across a wide temperature range from ambient to 870 °C, with minimal reliance on experimental data. The framework uses electron backscatter diffraction (EBSD) measurements to generate statistically equivalent microstructural volume elements to serve as inputs to the CP modeling framework. Stress–strain curves were generated from 1000 CP simulations, which serve as the training data set for the three ML regression algorithms explored: linear, extra-trees, and multi-layer perceptron. These three regression models were independently evaluated to compare their efficiency and identify the most suitable algorithm for the given problem. Results revealed that the extra-trees ML regressor outperforms the other models in both qualitative and quantitative aspects with an R 2 of 0.98. Subsequently, the differential evolution optimization approach is employed to calibrate the ML-based CP material parameters with experimental results obtained at various temperatures. Finally, temperature-dependent CP material parameters are formulated. The effectiveness and efficiency of the designed framework are validated through comparison with experimental results, demonstrating a high degree of agreement. These calibrated parametric constitutive equations enable further use of the CP model to study the deformation behavior of this alloy under a wide range of thermo-mechanical loading conditions.
36 MATERIALS SCIENCE↗
Grainsize evolution and differential comminution in an experimental regolith
The comminution of planetary surfaces by exposure to continuous meteorite bombardment was simulated by impacting the same fragmental gabbro target 200 times. The role of comminution and in situ gardening of planetary regoliths was addressed. Mean grain size continuously decreased with increasing shot number. Initially it decreased linearly with accumulated energy, but at some stage comminution efficiency started to decrease gradually. Point counting techniques, aided by the electron microprobe for mineral identification, were performed on a number of comminution products. Bulk chemical analyses of specific grain size fractions were also carried out. The finest sizes ( 10 microns) display generally the strongest enrichment/depletion factors. Similar, if not exactly identical, trends are reported from lunar soils. It is, therefore, not necessarily correct to explain the chemical characteristics of various grain sizes via different admixtures of materials from distant source terrains. Differential comminution of local source rocks may be the dominating factor.
Horz, F.↗
Chemical differentiation, thermal evolution, and catastrophic overturn on Venus: Predictions and geologic observations
Observations from Magellan show that: (1) the surface of Venus is generally geologically young, (2) there is no evidence for widespread recent crustal spreading or subduction, (3) the crater population permits the hypothesis that the surface is in production, and (4) relatively few impact craters appear to be embayed by volcanic deposits suggesting that the volcanic flux has drastically decreased as a function of time. These observations have led to consideration of hypotheses suggesting that the geological history of Venus may have changed dramatically as a function of time due to general thermal evolution, and/or thermal and chemical evolution of a depleted mantle layer, perhaps punctuated by catastrophic overturn of upper layers or episodic plate tectonics. We have previously examined the geological implications of some of these models, and here we review the predictions associated with two periods of Venus history. Stationary thick lithosphere and depleted mantle layer, and development of regional to global development of regional to global instabilities, and compare these predictions to the geological characteristics of Venus revealed by Magellan.
Head, James W.↗
Re-Os Isotopic Constraints on the Chemical Evolution and Differentiation of the Martian Mantle
The (187)Re-187Os isotopic systematics of SNC meteorites, thought to be from Mars, provide valuable information regarding the chemical processes that affected the Martian mantle, particularly with regard to the relative abundances of highly siderophile elements (HSE). Previously published data (Birck and Allegre 1994, Brandon et al. 2000), and new data obtained since these studies, indicate that the HSE and Os isotopic composition of the Martian mantle was primarily set in its earliest differentiation history. If so, then these meteorites provide key constraints on the processes that lead to variation in HSE observed in not only Mars, but also Earth, the Moon and other rocky bodies in the Solar System. Processes that likely have an effect on the HSE budgets of terrestrial mantles include core formation, magma ocean crystallization, development of juvenile crust, and the addition of a late veneer. Each of these processes will result in different HSE variation and the isotopic composition of mantle materials and mantle derived lavas. Two observations on the SNC data to present provide a framework for which to test the importance of each of these processes. First, the concentrations of Re and Os in SNC meteorites indicate that they are derived from a mantle that has similar concentrations to the Earth's mantle. Such an observation is consistent with a model where a chondritic late veneer replenished the Earth and Martian mantles subsequent to core formation on each planet. Alternative models to explain this observation do exist, but will require additional data to test the limitations of each. Second, Re-Os isotopic results from Brandon et al. (2000) and new data presented here, show that initial yos correlates with variations in the short-lived systems of (182)Hf- (182)W and (142)Sm-142Nd in the SNC meteorites (epsilon(sub W) and epsilon(sub 142Nd)). These systematics require an isolation of mantle reservoirs during the earliest differentiation history of Mars, and subsequent inefficient mixing between these reservoirs. These data show that models for the origin of isotopic variation for SNC meteorites require at least two long-lived mantle reservoirs, and possibly three. The range in the projected present day gamma(sub Os) of these reservoirs is from -5.4+/-2.6, to +4+/-1. The isotopic systematics of these reservoirs may be linked to development of cumulate crystal piles in a Martian magma ocean and variable amounts of late stage intercumulus melt. In this model, fractional crystallization of olivine and possibly other phases with slightly subchondritic Re/Os, from a solidifying magma ocean, resulted in a lower Re/Os ratio early cumulates, and a resultant low gamma(sub Os). Later cumulates or evolved melts crystallized with higher Re/Os ratios to produce the mantle reservoir(s) with consequent higher gamma(sub Os). Crystallization of the Martian magma ocean followed earliest core formation, as indicated by the correlation of epsilon(sub W) with epsilon(sub 142Nd) and initial gamma(sub Os).
Brandon, Alan D.↗
Advances in global optimization of high brightness beams
High brightness electron beams play an important role in accelerator-based applications such as driving X-ray free electron laser (FEL) radiation.Here in this paper, we report on advances in global beam dynamics optimization of an accelerator design using start-to-end simulations and a new parallel multi-objective differential evolution optimization method. The global optimization results in significant improvement of the final electron beam brightness.
43 PARTICLE ACCELERATORS↗
Multiclass Reduced-Set Support Vector Machines
There are well-established methods for reducing the number of support vectors in a trained binary support vector machine, often with minimal impact on accuracy. We show how reduced-set methods can be applied to multiclass SVMs made up of several binary SVMs, with significantly better results than reducing each binary SVM independently. Our approach is based on Burges' approach that constructs each reduced-set vector as the pre-image of a vector in kernel space, but we extend this by recomputing the SVM weights and bias optimally using the original SVM objective function. This leads to greater accuracy for a binary reduced-set SVM, and also allows vectors to be 'shared' between multiple binary SVMs for greater multiclass accuracy with fewer reduced-set vectors. We also propose computing pre-images using differential evolution, which we have found to be more robust than gradient descent alone. We show experimental results on a variety of problems and find that this new approach is consistently better than previous multiclass reduced-set methods, sometimes with a dramatic difference.
reduced set methods↗
Q-POP-Thermo: A general-purpose thermodynamics solver for ferroelectric materials
We report that Q-POP-Thermo is a program designed to compute thermodynamic monodomain equilibrium states and their properties for ferroelectric single crystals and thin films based on the Landau-Ginzburg-Devonshire (LGD) Theory. Utilizing symbolic manipulation with the SymPy Library, the governing equations along with appropriate boundary conditions are solved for speedy minimization of the free energy of a crystal. Utilizing the popular Differential Evolution algorithm, with appropriate hybridization, multiple phase diagrams, such as the pressure-temperature phase diagram for bulk single crystals and the common strain-temperature phase diagram for monodomain thin-film systems can be readily generated. Furthermore, a variety of material properties of stable ferroelectric phases, including dielectric, piezoelectric, and electrocaloric properties, can simultaneously be calculated. Validation studies are presented for both thin-film and single crystal systems to test the effectiveness and capability of the open-source program.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Complex modeling for the quantification of nanoscale disorder using genetic algorithms, density functional theory and line-profile analysis
Abstract not provided
Koch, Robert J.↗
Estimate of the Mass and Radial Profile of the Orphan–Chenab Stream's Dwarf-galaxy Progenitor Using MilkyWay@home
We fit the mass and radial profile of the Orphan–Chenab Stream's (OCS) dwarf-galaxy progenitor by using turnoff stars in the Sloan Digital Sky Survey and the Dark Energy Camera to constrain N-body simulations of the OCS progenitor falling into the Milky Way on the 1.5 PetaFLOPS MilkyWay@home distributed supercomputer. We infer the internal structure of the OCS's progenitor under the assumption that it was a spherically symmetric dwarf galaxy composed of a stellar system embedded in an extended dark matter halo. We optimize the evolution time, the baryonic and dark matter scale radii, and the baryonic and dark matter masses of the progenitor using a differential evolution algorithm. The likelihood score for each set of parameters is determined by comparing the simulated tidal stream to the angular distribution of OCS stars observed in the sky. We fit the total mass of the OCS's progenitor to (2.0 ± 0.3) × 10 7 M ⊙ with a mass-to-light ratio of γ = 73.5 ± 10.6 and (1.1 ± 0.2) × 10 6 M ⊙ within 300 pc of its center. Within the progenitor's half-light radius, we estimate a total mass of (4.0 ± 1.0) × 10 5 M ⊙ . We also fit the current sky position of the progenitor's remnant to be (α, δ) = ((166.0 ± 0.9)°, (–11.1 ± 2.5)°) and show that it is gravitationally unbound at the present time. The measured progenitor mass is on the low end of previous measurements and, if confirmed, lowers the mass range of ultrafaint dwarf galaxies. Our optimization assumes a fixed Milky Way potential, OCS orbit, and radial profile for the progenitor, ignoring the impact of the Large Magellanic Cloud.
79 ASTRONOMY AND ASTROPHYSICS↗
The nonlinear evolution of inviscid Goertler vortices in three-dimensional boundary layers
The nonlinear development of inviscid Gortler vortices in a three-dimensional boundary layer is considered. We do not follow the classical approach of weakly nonlinear stability problems and consider a mode which has just become unstable. Instead we extend the method of Blackaby, Dando, and Hall (1992), which considered the closely related nonlinear development of disturbances in stratified shear flows. The Gortler modes we consider are initially fast growing and we assume, following others, that boundary-layer spreading results in them evolving in a linear fashion until they reach a stage where their amplitudes are large enough and their growth rates have diminished sufficiently so that amplitude equations can be derived using weakly nonlinear and non-equilibrium critical-layer theories. From the work of Blackaby, Dando and Hall (1993) is apparent, given the range of parameters for the Gortler problem, that there are three possible nonlinear integro-differential evolution equations for the disturbance amplitude. These are a cubic due to viscous effects, a cubic which corresponds to the novel mechanism investigated in this previous paper, and a quintic. In this paper we shall concentrate on the two cubic integro-differential equations and in particular, on the one due to the novel mechanism as this will be the first to affect a disturbance. It is found that the consideration of a spatial evolution problem as opposed to temporal (as was considered in Blackaby, Dando, and Hall, 1992) causes a number of significant changes to the evolution equations.
Blackaby, Nicholas↗
The fully nonlinear development of Goertler vortices in growing boundary layers
The fully nonlinear development of small wavelength Goertler vortices in a growing boundary layer is investigated using a combination of asymptotic and numerical methods. The starting point for the analysis is the weakly nonlinear theory of Hall (1982b) who discussed the initial development of small amplitude vortices in a neighborhood of the location where they first become linearly unstable. That development is unusual in the context of nonlinear stability theory in that it is not described by the Stuart-Watson approach. In fact the development is governed by a pair of coupled nonlinear partial differential evolution equations for the vortex flow and the mean flow correction. Here the further development of this interaction is considered for vortices so large that the mean flow correction driven by them is as large as the basic state. Surprisingly it is found that such a nonlinear interaction can still be described by asymptotic means. It is shown that the vortices spread out across the boundary layer and effectively drive the boundary layer. In fact the system obtained by writing down the equations for the fundamental component of the vortex generate a differential equation for the basic state. Thus the mean flow adjusts so as to make these large amplitude vortices locally neutral. Moreover in the region where the vortices exist the mean flow has a square-root profile and the vortex velocity field can be written down in closed form. The upper and lower boundaries of the region of vortex activity are determined by a free-boundary problem involving the boundary layer equations. In general it is found that this region ultimately includes almost all of the original boundary layer and much of the free-stream. In this situation the mean flow has essentially no relationship to the flow which exists in the absence of the vortices.
Hall, Philip↗
Inferring fracture dilation and shear slip from surface deformation utilising trained surrogate models
An important task in energy and CO 2 storage (sequestration) in the subsurface is to verify that the surrounding fractures and faults are not activated, acting as leakage pathways. This is achievable through effective and efficient Measurement, Monitoring and Verification (MMV) plans. In this work, two surrogate models are trained to captures dilation (opening) and shear deformation of fractures, and the associated surface deformation. The trained surrogate model, based on conditional Generative-Adversarial Networks (cGAN) receives fracture apertures from dilational fractures together with fracture slips from shear fractures and predicts the combined surface deformation. An inversion algorithm based on Bayesian framework is proposed to identify the geometry of both types of fractures, as well as volume of dilational fractures and deformation moment induced by shear fractures, all from the measured surface deformation data. The inversion algorithm utilises the Differential Evolution (DE) optimisation technique that has the superior performance in finding the global minimum of cost function. The proposed surrogate-assisted inversion successfully inferred the unknown dip, dip direction and the volume of the dilational fractures as well as the induced deformation moment in shear fractures. The model was further tested for the inversion of a field hydraulic fracturing tilt dataset applying different scenarios with varying unknowns to show the model's performance, as well as incorporating shear deformation for better match with the observed data.
Dilation and shear↗
An ion specific continuum model on ionic surfactant's binary phase diagram in aqueous solution
Here, a continuum aggregation model is proposed to account for specific ion effects, enabling accurate phase diagrams calculations for the micellar, cylindrical and lamellar aggregates of sodium/potassium carboxylate surfactants in aqueous solution across a range of temperatures. Three groups of concentrations at distinctive temperatures are fitted to build empirical temperature dependence given the limited availability of relevant experimental measurements. The specific ion effects are manifested in the aggregates' surface tension as well as in the distributions of counter-ions' concentrations in the vicinity of the aggregates. The aggregates' geometric sizes are well-reproduced. The differential evolution algorithm is applied to address boundary conditions of the electrostatic potential, aggregate size optimization as well as the equilibrium of monomers transferring between the aggregate and aqueous region.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗