SEARCH · Search NASA

Results for “boundary value problems”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,441 records · Page 80

All sources NASA NTRS DOE OSTI

DOE OSTI2017

The structure and propagation of laminar flames under autoignitive conditions

The laminar flame speed s l is an important reference quantity for characterising and modelling combustion. Experimental measurements of laminar flame speed require the residence time of the fuel/air mixture (τ f ) to be shorter than the autoignition delay time (τ). This presents a considerable challenge for conditions where autoignition occurs rapidly, such as in compression ignition engines. As a result, experimental measurements in typical compression ignition engine conditions do not exist. Simulations of freely propagating premixed flames, where the burning velocity is found as an eigenvalue of the solution, are also not well posed in such conditions, since the mixture ahead of the flame can autoignite, leading to the so called “cold boundary problem”. In this paper, a numerical method for estimating a reference flame speed, s R , is proposed that is valid for laminar flame propagation at autoignitive conditions. Two isomer fuels are considered to test this method: ethanol, which in the considered conditions is a single-stage ignition fuel; and dimethyl ether, which has a temperature-dependent single- or two-stage ignition and a negative temperature coefficient regime for τ. Calculations are performed for the flame position in a one-dimensional computational domain with inflow-outflow boundary conditions, as a function of the inlet velocity U I and for stoichiometric fuel–air premixtures. The response of the flame position, L F , to U I shows distinct stabilisation regimes. For single-stage ignition fuels, at low U I the flame speed exceeds U I and the flame becomes attached to the inlet. Above a critical U I value, the flame detaches from the inlet and L f becomes extremely sensitive to U I until, for sufficiently high U I , the sensitivity decreases and L f corresponds to the location expected from a purely autoignition stabilised flame. The transition from the attached to the autoignition regimes has a corresponding peak dL f /dU I value which is proposed to be a unique reference flame speed s R for single-stage ignition fuels. For two-stage ignition fuels, there is an additional stable regime where a high-temperature flame propagates into a pool of combustion intermediates generated by the first stage of autoignition. This results in two peaks in dL f /dU I and therefore two reference flame speed values. The lower value corresponds to the definition of s R for single-stage ignition fuels, while the higher value exists only for two-stage ignition fuels and corresponds to a high temperature flame propagating into the first stage of autoignition and is denoted s R ' . Finally, a transport budget analysis for low- and high-temperature radical species is also performed, which confirms that the flame structures at U I = s R and U I = s R ' do indeed correspond to premixed flames (deflagrations), as opposed to spontaneous ignition fronts which do not have a unique propagation speed.

33 ADVANCED PROPULSION SYSTEMS↗

NASA NTRS1989

Confidence set inference with a prior quadratic bound

In the uniqueness part of a geophysical inverse problem, the observer wants to predict all likely values of P unknown numerical properties z=(z sub 1,...,z sub p) of the earth from measurement of D other numerical properties y (sup 0) = (y (sub 1) (sup 0), ..., y (sub D (sup 0)), using full or partial knowledge of the statistical distribution of the random errors in y (sup 0). The data space Y containing y(sup 0) is D-dimensional, so when the model space X is infinite-dimensional the linear uniqueness problem usually is insoluble without prior information about the correct earth model x. If that information is a quadratic bound on x, Bayesian inference (BI) and stochastic inversion (SI) inject spurious structure into x, implied by neither the data nor the quadratic bound. Confidence set inference (CSI) provides an alternative inversion technique free of this objection. Confidence set inference is illustrated in the problem of estimating the geomagnetic field B at the core-mantle boundary (CMB) from components of B measured on or above the earth's surface.

Backus, George E.↗

NASA NTRS1992

Airborne lidar and radiometric observations of PBL- and low clouds

Boundary layer- and low altitude clouds over open ocean and continent areas have been studied during several field campaigns since mid-1990 using the French airborne backscatter lidar LEANDRE in conjunction with on-board IR and visible radiometers. LEANDRE is an automatic system, and a modification of the instrumental parameters, when airborne, is computer controlled through an operator keyboard. The vertical range squared lidar signals and instrument status are displayed in real time on two dedicated monitors. The lidar is used either down- or up-looking while the aircraft is flying above or below clouds. A switching of the viewing configuration takes about a minute. The lidar measurements provide a high resolution description of cloud morphology and holes in cloud layers. The flights were conducted during various meteorological conditions on single or multilayer stratocumulus and cumulus decks. Analysis on a single shot basis of cloud top (or bottom) altitude and a plot of the corresponding histogram allows one to determine a probability density function (PDF). The preliminary results show the PDFs for cloud top are not Gaussian and symmetric about the mean value. The skewness varies with atmospheric conditions. An example of results recorded over the Atlantic ocean near Biarritz is displayed, showing: (1) the range squared lidar signals as a function of time (here 100 s corresponds to about 8 km, 60 shots are averaged on horizontal); the Planetary Boundary Layer (PBL) - up to 600 m - is observed at the beginning of the leg as well as on surface returns, giving an indication of the porosity; (2) the cloud top altitude variation between 2.4 to 2.8 km during the 150 to 320 s section; and (3) the corresponding PDF. Similar results are obtained on stratocumulus over land. Single shot measurements can be used also to determine an optical porosity at a small scale as well as a fractional cloudiness at a larger scale. A comparison of cloud top altitude retrieved from lidar and narrowbeam IR radiometer is conducted to study the scale integration problem. A good agreement within less than 100 m relies on spatial uniformity and an optically thick layer. In the presence of holes, a discrepancy is observed. This is illustrated in figure 2, displaying as a function of time (1) the lidar signals; (2) the target temperature (either clouds or sea surface) retreived from a narrowbeam IR radiometer, 17 C is the sea surface temperature on that day; and (3) the visible flux, linked to cloud albedo, measured by a pyranometer. In preparation of ASTEX, down- and up-looking measurements where conducted on stratocumulus clouds over the Atlantic Ocean near Quimper in Brittany. Depending on the flight pattern orientation with respect to the wind, the top and bottom cloud morphologies are different. Preliminary results are given on cloud morphology, cloud top PDFs, optical porosity, fractional cloudiness, and comparison of lidar and radiometric measurements.

Flamant, P. H.↗

DOE OSTI2018

Minimal area surfaces in AdS n+1 and Wilson loops

The AdS/CFT correspondence relates the expectation value of Wilson loops in N = 4 SYM to the area of minimal surfaces in AdS 5 . In this paper we consider minimal area surfaces in generic Euclidean AdS n+1 using the Pohlmeyer reduction in a similar way as we did previously in Euclidean AdS 3 . As in that case, the main obstacle is to find the correct parameterization of the curve in terms of a conformal parameter. Once that is done, the boundary conditions for the Pohlmeyer fields are obtained in terms of conformal invariants of the curve. After solving the Pohlmeyer equations, the area can be expressed as a boundary integral involving a generalization of the conformal arc-length, curvature and torsion of the curve. Furthermore, one can introduce the λ-deformation symmetry of the contours by a simple change in the conformal invariants. This determines the λ-deformed contours in terms of the solution of a boundary linear problem. In fact the condition that all λ deformed contours are periodic can be used as an alternative to solving the Pohlmeyer equations and is equivalent to imposing the vanishing of an infinite set of conserved charges derived from integrability.

't Hooft and Polyakov loops↗

NASA NTRS1997

Multidisciplinary Optimization of Tilt Rotor Blades Using Comprehensive Composite Modeling Technique

An optimization procedure is developed for addressing the design of composite tilt rotor blades. A comprehensive technique, based on a higher-order laminate theory, is developed for the analysis of the thick composite load-carrying sections, modeled as box beams, in the blade. The theory, which is based on a refined displacement field, is a three-dimensional model which approximates the elasticity solution so that the beam cross-sectional properties are not reduced to one-dimensional beam parameters. Both inplane and out-of-plane warping are included automatically in the formulation. The model can accurately capture the transverse shear stresses through the thickness of each wall while satisfying stress free boundary conditions on the inner and outer surfaces of the beam. The aerodynamic loads on the blade are calculated using the classical blade element momentum theory. Analytical expressions for the lift and drag are obtained based on the blade planform with corrections for the high lift capability of rotor blades. The aerodynamic analysis is coupled with the structural model to formulate the complete coupled equations of motion for aeroelastic analyses. Finally, a multidisciplinary optimization procedure is developed to improve the aerodynamic, structural and aeroelastic performance of the tilt rotor aircraft. The objective functions include the figure of merit in hover and the high speed cruise propulsive efficiency. Structural, aerodynamic and aeroelastic stability criteria are imposed as constraints on the problem. The Kreisselmeier-Steinhauser function is used to formulate the multiobjective function problem. The search direction is determined by the Broyden-Fletcher-Goldfarb-Shanno algorithm. The optimum results are compared with the baseline values and show significant improvements in the overall performance of the tilt rotor blade.

Chattopadhyay, Aditi↗

DOE OSTI2022

FROSch Preconditioners for Land Ice Simulations of Greenland and Antarctica

Numerical simulations of Greenland and Antarctic ice sheets involve the solution of large-scale highly nonlinear systems of equations on complex shallow geometries. This work is concerned with the construction of Schwarz preconditioners for the solution of the associated tangent problems, which are challenging for solvers mainly because of the strong anisotropy of the meshes and wildly changing boundary conditions that can lead to poorly constrained problems on large portions of the domain. In this study, two-level generalized Dryja--Smith--Widlund (GDSW)--type Schwarz preconditioners are applied to different land ice problems, i.e., a velocity problem, a temperature problem, as well as the coupling of the former two problems. We employ the message passing interface (MPI)--parallel implementation of multilevel Schwarz preconditioners provided by the package FROSch (fast and robust Schwarz) from the Trilinos library. The strength of the proposed preconditioner is that it yields out-of-the-box scalable and robust preconditioners for the single physics problems. To the best of our knowledge, this is the first time two-level Schwarz preconditioners have been applied to the ice sheet problem and a scalable preconditioner has been used for the coupled problem. The preconditioner for the coupled problem differs from previous monolithic GDSW preconditioners in the sense that decoupled extension operators are used to compute the values in the interior of the subdomains. Additionally, several approaches for improving the performance, such as reuse strategies and shared memory OpenMP parallelization, are explored. In our numerical study we target both uniform meshes of varying resolution for the Antarctic ice sheet as well as nonuniform meshes for the Greenland ice sheet. We present several weak and strong scaling studies confirming the robustness of the approach and the parallel scalability of the FROSch implementation. Among the highlights of the numerical results are a weak scaling study for up to 32,K processor cores (8,K MPI ranks and 4 OpenMP threads) and 566,M degrees of freedom for the velocity problem as well as a strong scaling study for up to 4,K processor cores (and MPI ranks) and 68,M degrees of freedom for the coupled problem.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2014

Adjustable magnetoelectric effect of self-assembled vertical multiferroic nanocomposite films by the in-plane misfit strain and ferromagnetic volume fraction

The strain-mediated magnetoelectric (ME) property of self-assembled vertical multiferroic nanocomposite films epitaxially grown on cubic substrates was calculated by a nonlinear thermodynamic theory combined with the elastic theory. The dependent relations of phase state of ferroelectric films with the in-plane misfit strain, out-of-plane misfit strain, temperature, and volume fraction of ferromagnetic phase were confirmed. The effects of in-plane misfit strain and ferromagnetic volume fraction on the polarization and dielectric constant of ferroelectric films at room temperature were elaborately analyzed for the vertical BaTiO{sub 3}-CoFe{sub 2}O{sub 4} and PbTiO{sub 3}-CoFe{sub 2}O{sub 4} nanocomposite films. Our calculated results confirmed the relationship among ME effect and in-plane misfit strain and ferromagnetic volume fraction in the nanocomposite films. The ME voltage coefficients of vertical BaTiO{sub 3}-CoFe{sub 2}O{sub 4} and PbTiO{sub 3}-CoFe{sub 2}O{sub 4} nanocomposite films displayed various maximums and abrupt points at special phases and phase transition boundaries. The ME voltage coefficients of lead-free BaTiO{sub 3}-CoFe{sub 2}O{sub 4} nanocomposite films epitaxially grown on different substrates could reach a comparative value of ∼2 V·cm{sup −1}·Oe{sup −1} under the controllable in-plane misfit strain induced by substrate clamping. Our results provided an available method for the optimal design of vertical multiferroic nanocomposites with adjustable ME effect by optimizing the ferromagnetic volume fraction and substrate type.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2017

Optical properties of magnetic photonic crystals with an arbitrary magnetization orientation

We have studied the peculiarities of diffraction of light in magnetic photonic crystals at large values of magnetooptical activity parameter and modulation depth. We have considered the case of an arbitrary angle between the directions of the external static magnetic field and the normal to the layer. The problem has been solved by the modified Ambartsumyan layer summation method. It has been shown that the given system is nonreciprocal with respect to not only circular, but linear polarizations also. In this case, a new type of nonreciprocity is observed (namely, the relation R(α) ≠ R(–α) holds, where R is the reflection coefficient and α is the angle of incidence). It has been demonstrated that in the case of oblique incidence, there appears a new photonic forbidden band that is not selective relative to the polarization of incident light. We have detected strong dependences of reflectance, absorbance, transmittance nonreciprocity, and other characteristics on the angle between the direction of the external static magnetic field and the normal to the layer boundary. Such a system can be used as a controllable polarization filter and a mirror, as well as a source of circular (elliptic) polarization, a controllable optical diode, and so on.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

NASA NTRS1996

Heat and Momentum Transfer Studies in High Reynolds Number Wavy Films at Normal and Reduced Gravity Conditions

We examined the effect of the gas flow on the liquid film when the gas flows in the countercurrent direction in a vertical pipe at normal gravity conditions. The most dramatic effect of the simultaneous flow of gas and liquid in pipes is the greatly increased transport rates of heat, mass, and momentum. In practical situations this enhancement can be a benefit or it can result in serious operational problems. For example, gas-liquid flow always results in substantially higher pressure drop and this is usually undesirable. However, much higher heat transfer coefficients can be expected and this can obviously be of benefit for purposes of design. Unfortunately, designers know so little of the behavior of such two phase systems and as a result these advantages are not utilized. Due to the complexity of the second order boundary model as well as the fact that the pressure variation across the film is small compared to the imposed gas phase pressure, the countercurrent gas flow affect was studied for the standard boundary layer model. A different stream function that can compensate the shear stress affect was developed and this stream function also can predict periodic solutions. The discretized model equations were transformed to a traveling wave coordinate system. A stability analysis of these sets of equations showed the presence of a Hopf bifurcation for certain values of the traveling wave velocity and the shear stress. The Hopf celerity was increased due to the countercurrent shear. For low flow rate the increases of celerity are more than for the high flow rate, which was also observed in experiments. Numerical integration of a traveling wave simplification of the model also predicts the existence of chaotic large amplitude, nonperiodic waves as observed in the experiments. The film thickness was increased by the shear.

Balakotaiah, V.↗

DOE OSTI2023

Sensitivity analysis of a layered piezoelectric system using ZFEM

The complex variable finite element method (ZFEM) is a numerical technique which aims to find the partial derivatives of the independent variables with respect to variation in dependent parameters declared in the physics. This is done by combining the complex Taylor series expansion within the weak formulation of the governing equation in a coupled system of linear equations forming a complex valued block matrix given by the Cauchy–Riemann matrix representation. In this work, two-dimensional linear first-order elements have been implemented in ZFEM to predict the design derivatives of the mechanical displacement field and the voltage potential field for a layered piezoelectric system in a steady-state study with Dirichlet boundary condition applied at the top and bottom edges of the geometry. This approach allows the standard FEM solution to quantify the sensitivity of the mechanical displacement and voltage potential fields with respect to small variations in the material properties through the information obtained from the computation of the derivatives. The domain is formed by a layered body with PZT-4 and PZT-5 stacked together. For result verification, the numerical solution obtained with ZFEM was compared to results from a commercial FEM package and the solution from the imaginary part was compared to the exact solution of a well-known benchmark problem. In conclusion, comparison of the results showed good agreement for both the real and imaginary parts of the solution and the largest sensitivities were found in PZT-5 specifically in C 13 , C 33 , and ε 33 .

42 ENGINEERING↗

NASA NTRS2016

Avoiding Stair-Step Artifacts in Image Registration for GOES-R Navigation and Registration Assessment

In developing software for independent verification and validation (IVV) of the Image Navigation and Registration (INR) capability for the Geostationary Operational Environmental Satellite R Series (GOES-R) Advanced Baseline Imager (ABI), we have encountered an image registration artifact which limits the accuracy of image offset estimation at the subpixel scale using image correlation. Where the two images to be registered have the same pixel size, subpixel image registration preferentially selects registration values where the image pixel boundaries are close to lined up. Because of the shape of a curve plotting input displacement to estimated offset, we call this a stair-step artifact. When one image is at a higher resolution than the other, the stair-step artifact is minimized by correlating at the higher resolution. For validating ABI image navigation, GOES-R images are correlated with Landsat-based ground truth maps. To create the ground truth map, the Landsat image is first transformed to the perspective seen from the GOES-R satellite, and then is scaled to an appropriate pixel size. Minimizing processing time motivates choosing the map pixels to be the same size as the GOES-R pixels. At this pixel size image processing of the shift estimate is efficient, but the stair-step artifact is present. If the map pixel is very small, stair-step is not a problem, but image correlation is computation-intensive. This paper describes simulation-based selection of the scale for truth maps for registering GOES-R ABI images.

stair-step artifact↗

NASA NTRS2016

Avoiding Stair-Step Artifacts in Image Registration for GOES-R Navigation and Registration Assessment

image registration↗

DOE OSTI2018

A novel numerical treatment of the near-wall regions in the k – ω class of RANS models

In this study, we discuss a novel approach to modeling the near-wall region in the class of k – ω to models. The proposed methodology obviates the need for ad hoc boundary conditions of to omega on the wall as typically required in the k – ω model. The primary motivation of this work is to provide a formulation equivalent to the standard k – ω to and k – ω to SST models, but which at the same time overcomes their limitations in the context of their implementation in high-order methods. This is achieved by subtracting the asymptotically known singular behavior of omega at walls. Imposing a grid-dependent value of omega at walls in high-order codes is not straightforward as is demonstrated below and it causes instability as well as accuracy issues. The mathematical formulation of the two novel approaches, termed as the "regularized k – ω model" and the "regularized k – ω OJ SST model", is discussed in detail. A consistency and verification study for these two approaches is performed by proving that the regularized models recover the results of the standard models in various canonical problems, such as turbulent channel and pipe flows, and a systematic investigation of convergence, using both p- (polynomial order) as well as h- (grid) refinement is reported. Furthermore, comparisons highlighting the performance of the proposed methods in more complex configurations such as flow over a backward facing step and the turbulent mixing of fluid streams of different temperatures in a T-junction are also presented.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2017

GeneratorSE: A Sizing Tool for Variable-Speed Wind Turbine Generators

This report documents a set of analytical models employed by the optimization algorithms within the GeneratorSE framework. The initial values and boundary conditions employed for the generation of the various designs and initial estimates for basic design dimensions, masses, and efficiency for the four different models of generators are presented and compared with empirical data collected from previous studies and some existing commercial turbines. These models include designs applicable for variable-speed, high-torque application featuring direct-drive synchronous generators and low-torque application featuring induction generators. In all of the four models presented, the main focus of optimization is electromagnetic design with the exception of permanent-magnet and wire-wound synchronous generators, wherein the structural design is also optimized. Thermal design is accommodated in GeneratorSE as a secondary attribute by limiting the winding current densities to acceptable limits. A preliminary validation of electromagnetic design was carried out by comparing the optimized magnetic loading against those predicted by numerical simulation in FEMM4.2, a finite-element software for analyzing electromagnetic and thermal physics problems for electrical machines. For direct-drive synchronous generators, the analytical models for the structural design are validated by static structural analysis in ANSYS.

17 WIND ENERGY↗

DOE OSTI2015

THE INFLUENCE OF SPATIAL RESOLUTION ON NONLINEAR FORCE-FREE MODELING

The nonlinear force-free field (NLFFF) model is often used to describe the solar coronal magnetic field, however a series of earlier studies revealed difficulties in the numerical solution of the model in application to photospheric boundary data. We investigate the sensitivity of the modeling to the spatial resolution of the boundary data, by applying multiple codes that numerically solve the NLFFF model to a sequence of vector magnetogram data at different resolutions, prepared from a single Hinode/Solar Optical Telescope Spectro-Polarimeter scan of NOAA Active Region 10978 on 2007 December 13. We analyze the resulting energies and relative magnetic helicities, employ a Helmholtz decomposition to characterize divergence errors, and quantify changes made by the codes to the vector magnetogram boundary data in order to be compatible with the force-free model. This study shows that NLFFF modeling results depend quantitatively on the spatial resolution of the input boundary data, and that using more highly resolved boundary data yields more self-consistent results. The free energies of the resulting solutions generally trend higher with increasing resolution, while relative magnetic helicity values vary significantly between resolutions for all methods. All methods require changing the horizontal components, and for some methods also the vertical components, of the vector magnetogram boundary field in excess of nominal uncertainties in the data. The solutions produced by the various methods are significantly different at each resolution level. We continue to recommend verifying agreement between the modeled field lines and corresponding coronal loop images before any NLFFF model is used in a scientific setting.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2021

Deep Learning Approaches to Surrogates for Solving the Diffusion Equation for Mechanistic Real-World Simulations

In many mechanistic medical, biological, physical, and engineered spatiotemporal dynamic models the numerical solution of partial differential equations (PDEs), especially for diffusion, fluid flow and mechanical relaxation, can make simulations impractically slow. Biological models of tissues and organs often require the simultaneous calculation of the spatial variation of concentration of dozens of diffusing chemical species. One clinical example where rapid calculation of a diffusing field is of use is the estimation of oxygen gradients in the retina, based on imaging of the retinal vasculature, to guide surgical interventions in diabetic retinopathy. Furthermore, the ability to predict blood perfusion and oxygenation may one day guide clinical interventions in diverse settings, i.e., from stent placement in treating heart disease to BOLD fMRI interpretation in evaluating cognitive function (Xie et al., 2019; Lee et al., 2020). Since the quasi-steady-state solutions required for fast-diffusing chemical species like oxygen are particularly computationally costly, we consider the use of a neural network to provide an approximate solution to the steady-state diffusion equation. Machine learning surrogates, neural networks trained to provide approximate solutions to such complicated numerical problems, can often provide speed-ups of several orders of magnitude compared to direct calculation. Surrogates of PDEs could enable use of larger and more detailed models than are possible with direct calculation and can make including such simulations in real-time or near-real time workflows practical. Creating a surrogate requires running the direct calculation tens of thousands of times to generate training data and then training the neural network, both of which are computationally expensive. Often the practical applications of such models require thousands to millions of replica simulations, for example for parameter identification and uncertainty quantification, each of which gains speed from surrogate use and rapidly recovers the up-front costs of surrogate generation. We use a Convolutional Neural Network to approximate the stationary solution to the diffusion equation in the case of two equal-diameter, circular, constant-value sources located at random positions in a two-dimensional square domain with absorbing boundary conditions. Such a configuration caricatures the chemical concentration field of a fast-diffusing species like oxygen in a tissue with two parallel blood vessels in a cross section perpendicular to the two blood vessels. To improve convergence during training, we apply a training approach that uses roll-back to reject stochastic changes to the network that increase the loss function. The trained neural network approximation is about 1000 times faster than the direct calculation for individual replicas. Because different applications will have different criteria for acceptable approximation accuracy, we discuss a variety of loss functions and accuracy estimators that can help select the best network for a particular application. We briefly discuss some of the issues we encountered with overfitting, mismapping of the field values and the geometrical conditions that lead to large absolute and relative errors in the approximate solution.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2015

Blockage of domain growth by nanoscale heterogeneities in a relaxor ferroelectric Sr{sub 0.61}Ba{sub 0.39}Nb{sub 2}O{sub 6}

The growth of localized subsurface domains in a relaxor ferroelectric Sr{sub 0.61}Ba{sub 0.39}Nb{sub 2}O{sub 6} is studied using the technique of piezoresponse force microscopy (PFM). Ferroelectric domains are created by applying moderate voltages of 10–50 V to the conductive tip of a scanning force microscope brought into contact with a nonpolar face of a Sr{sub 0.61}Ba{sub 0.39}Nb{sub 2}O{sub 6} crystal. PFM images of written domains are acquired and analyzed quantitatively to determine the domain length along the polar axis and its width in the transverse direction. The dependences of domain sizes on the applied voltage, pulse duration, and the time passed after completion of the voltage pulse are reported and analyzed theoretically. It is shown that the observed kinetics of domain growth can be explained by the creep of domain boundaries occurring in the presence of random electric fields inherent in Sr{sub 0.61}Ba{sub 0.39}Nb{sub 2}O{sub 6}. The comparison of measured domain sizes with their equilibrium values calculated with the aid of the thermodynamic theory demonstrates that the growth of subsurface domains in Sr{sub 0.61}Ba{sub 0.39}Nb{sub 2}O{sub 6} is blocked by nanoscale heterogeneities characteristic of this relaxor ferroelectric. These results may have important implications for the development of nonlinear optical devices based on nanoheterogeneous ferroelectrics.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2017

Lithium Depletion in Solar-like Stars: Effect of Overshooting Based on Realistic Multi-dimensional Simulations

We study lithium depletion in low-mass and solar-like stars as a function of time, using a new diffusion coefficient describing extra-mixing taking place at the bottom of a convective envelope. This new form is motivated by multi-dimensional fully compressible, time-implicit hydrodynamic simulations performed with the MUSIC code. Intermittent convective mixing at the convective boundary in a star can be modeled using extreme value theory, a statistical analysis frequently used for finance, meteorology, and environmental science. In this Letter, we implement this statistical diffusion coefficient in a one-dimensional stellar evolution code, using parameters calibrated from multi-dimensional hydrodynamic simulations of a young low-mass star. We propose a new scenario that can explain observations of the surface abundance of lithium in the Sun and in clusters covering a wide range of ages, from ∼50 Myr to ∼4 Gyr. Because it relies on our physical model of convective penetration, this scenario has a limited number of assumptions. It can explain the observed trend between rotation and depletion, based on a single additional assumption, namely, that rotation affects the mixing efficiency at the convective boundary. We suggest the existence of a threshold in stellar rotation rate above which rotation strongly prevents the vertical penetration of plumes and below which rotation has small effects. In addition to providing a possible explanation for the long-standing problem of lithium depletion in pre-main-sequence and main-sequence stars, the strength of our scenario is that its basic assumptions can be tested by future hydrodynamic simulations.

79 ASTRONOMY AND ASTROPHYSICS↗