SEARCH · Search NASA

Results for “derivative discontinuity”

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 253 records · Page 14

All sources NASA NTRS DOE OSTI

NASA NTRS2020

Uncertainties in CERES Top-of-Atmosphere Fluxes Caused by Changes in Accompanying Imager

The Clouds and the Earth’s Radiant Energy System (CERES) project provides observations of Earth’s radiation budget using measurements from CERES instruments on board the Terra, Aqua, Suomi National Polar-orbiting Partnership (S-NPP), and NOAA-20 satellites. The CERES top-of-atmosphere (TOA) fluxes are produced by converting radiance measurements using empirical angular distribution models, which are functions of cloud properties that are retrieved from imagers flying with the CERES instruments. As the objective is to create a long-term climate data record, not only calibration consistency of the six CERES instruments needs to be maintained for the entire time period, it is also important to maintain the consistency of other input data sets used to produce this climate data record. In this paper, we address aspects that could potentially affect the CERES TOA flux data quality. Discontinuities in imager calibration can affect cloud retrieval which can lead to erroneous flux trends. When imposing an artificial 0.6 per decade decreasing trend to cloud optical depth, which is similar to the trend difference between CERES Edition 2 and Edition 4 cloud retrievals, the decadal SW flux trend changed from −0.3 5 ± 0.18 Wm −2 to 0.61 ± 0.18 Wm −2 . This indicates that a 13% change in cloud optical depth results in about 1% change in the SW flux. Furthermore, different CERES instruments provide valid fluxes at different viewing zenith angle ranges, and including fluxes derived at the most oblique angels unique to S-NPP (>66 ∘ ) can lead to differences of 0.8 Wm −2 and 0.3 Wm −2 in global monthly mean instantaneous SW flux and LW flux. To ensure continuity, the viewing zenith angle ranges common to all CERES instruments (<66 ∘ ) are used to produce the long-term Earth’s radiation budget climate data record. The consistency of cloud properties retrieved from different imagers also needs to be maintained to ensure the TOA flux consistency.

Wenying Su↗

NASA NTRS1991

Integrated trajectory and control analysis for generic hypersonic vehicles

A tool which supports integrated assessment of air-breathing hypersonic vehicle trajectory management and control strategy has been developed by Hattis and Malchow (1991). Using a two-point boundary value problem solution technique, the tool can construct an entire near-minimum-fuel trajectory and desired control strategy from takeoff to orbit, while satisfying physically derived inequality constraints and while achieving efficient propulsive mode phasing. A strategy is also formulated to construct preliminary trajectory and control history representations with less computational burden than required for the overall flight profile assessment. A tabulated example hypersonic vehicle model has been used to demonstrate the integrated analysis methodology. To assure good numerical behavior when the algorithm uses tabulated model data, a multidimensional cubic spline curve smoothing routine has been developed. A dynamic-pressure-constrained, near-fuel-optimal trajectory was computed to provide examples of information that the methodology can provide. Previously unspecified propulsive discontinuities were located, flight regimes demanding rapid attitude changes were identified, available and closed-loop controller authority was ascertained, and inadequacies in vehicle and specific subsystem model representations were determined.

Hattis, Philip D.↗

NASA NTRS1992

Identification of time-varying structural dynamic systems - An artificial intelligence approach

An application of the artificial intelligence-derived methodologies of heuristic search and object-oriented programming to the problem of identifying the form of the model and the associated parameters of a time-varying structural dynamic system is presented in this paper. Possible model variations due to changes in boundary conditions or configurations of a structure are organized into a taxonomy of models, and a variant of best-first search is used to identify the model whose simulated response best matches that of the current physical structure. Simulated model responses are verified experimentally. An output-error approach is used in a discontinuous model space, and an equation-error approach is used in the parameter space. The advantages of the AI methods used, compared with conventional programming techniques for implementing knowledge structuring and inheritance, are discussed. Convergence conditions and example problems have been discussed. In the example problem, both the time-varying model and its new parameters have been identified when changes occur.

Glass, B. J.↗

NASA NTRS1950

Hydrodynamic Properties of Planing Surfaces and Flying Boats

The study of the hydrodynamic properties of planing bottom of flying boats and seaplane floats is at the present time based exclusively on the curves of towing tests conducted in tanks. In order to provide a rational basis for the test procedure in tanks and practical design data, a theoretical study must be made of the flow at the step and relations derived that show not only qualitatively but quantitatively the inter-relations of the various factors involved. The general solution of the problem of the development of hydrodynamic forces during the motion of the seaplane float or flying boat is very difficult for it is necessary to give a three-dimensional solution, which does not always permit reducing the analysis to the form of workable computation formulas. On the other had, the problem is complicated by the fact that the object of the analysis is concerned with two fluid mediums, namely, air and water, which have a surface of density discontinuity between them. The theoretical and experimental investigations on the hydrodynamics of a ship cannot be completely carried over to the design of floats and flying-boat hulls, because of the difference in the shape of the contour lines of the bodies, and, because of the entirely different flow conditions from the hydrodynamic viewpoint.

Sokolov, N. A.↗

NASA NTRS2004

Sensitivity Equation Derivation for Transient Heat Transfer Problems

The focus of the paper is on the derivation of sensitivity equations for transient heat transfer problems modeled by different discretization processes. Two examples will be used in this study to facilitate the discussion. The first example is a coupled, transient heat transfer problem that simulates the press molding process in fabrication of composite laminates. These state equations are discretized into standard h-version finite elements and solved by a multiple step, predictor-corrector scheme. The sensitivity analysis results based upon the direct and adjoint variable approaches will be presented. The second example is a nonlinear transient heat transfer problem solved by a p-version time-discontinuous Galerkin's Method. The resulting matrix equation of the state equation is simply in the form of Ax = b, representing a single step, time marching scheme. A direct differentiation approach will be used to compute the thermal sensitivities of a sample 2D problem.

Hou, Gene↗

NASA NTRS1985

Use of simultaneous radionuclide and other trace constituent measurements to determine fluxes and atmospheric residence times

Finding a way to use atmospheric concentrations to determine source strengths and atmospheric residence times pose considerable difficulties. An approach is to use simultaneously measured radionuclide concentrations as a means of deriving fluxes and residence times for certain atmospheric trace constituents of interest. The study and use of such methods in a pilot field experiment were among the original goals of this program. Four such approaches were studied or developed; suitable sampling and analysis systems were built and tested or are in the final stages of construction; and a small scale field experiment. The evolution is reviewed of the program and the field experiment is described along with the possibilities for its eventual extension to include other species of interest. The four radionuclide based approaches studied are: Mean vertical profiles; Trajectory correlation; Boundary layer discontinuity; and Lead deposition.

Kritz, Mark↗

NASA NTRS1999

The Analysis and Design of Low Boom Configurations Using CFD and Numerical Optimization Techniques

The use of computational fluid dynamics (CFD) for the analysis of sonic booms generated by aircraft has been shown to increase the accuracy and reliability of predictions. CFD takes into account important three-dimensional and nonlinear effects that are generally neglected by modified linear theory (MLT) methods. Up to the present time, CFD methods have been primarily used for analysis or prediction. Some investigators have used CFD to impact the design of low boom configurations using trial and error methods. One investigator developed a hybrid design method using a combination of Modified Linear Theory (e.g. F-functions) and CFD to provide equivalent area due to lift driven by a numerical optimizer to redesign or modify an existing configuration to achieve a shaped sonic boom signature. A three-dimensional design methodology has not yet been developed that completely uses nonlinear methods or CFD. Constrained numerical optimization techniques have existed for some time. Many of these methods use gradients to search for the minimum of a specified objective function subject to a variety of design variable bounds, linear and nonlinear constraints. Gradient based design optimization methods require the determination of the objective function gradients with respect to each of the design variables. These optimization methods are efficient and work well if the gradients can be obtained analytically. If analytical gradients are not available, the objective gradients or derivatives with respect to the design variables must be obtained numerically. To obtain numerical gradients, say, for 10 design variables, might require anywhere from 10 to 20 objective function evaluations. Typically, 5-10 global iterations of the optimizer are required to minimize the objective function. In terms of using CFD as a design optimization tool, the numerical evaluation of gradients can require anywhere from 100 to 200 CFD computations per design for only 10 design variables. If one CFD computation requires an hour of computational time on a Cray computer, one can see that the use of constrained numerical optimization quickly becomes impractical.Hence, in order to practically couple a numerical design optimization technique with a CFD method, the CFD method must be extremely efficient with running times on the order of only minutes. The CFD Euler code developed under NASA sponsorship and referred to as MIM3D-SB for the most part fulfills these efficiency requirements. Analysis of wing- body configurations can be computed in a matter of a few minutes. The present study will concentrate on the feasibility of the use of this CFD code in conjunction with a numerical design optimization technique for the sonic boom reduction of candidate HSCT configurations. A preliminary supersonic aircraft design system has been established that utilizes the numerical design optimization code NPSOL developed at Stanford University coupled with the supersonic NUM3D-SB CFD code. Many questions still need to be answered in regard to using CFD and numerical optimizers as design tools. There are difficulties related to both the CFD codes and the numerical optimizers. Numerical optimizers can converge to a local minima rather than a global minima. This behavior is largely a function of the initial guess in the design space. The optimizer also is searching for a minimum of the function in terms of its derivative without any regard to the actual function value. Numerically (i.e. CFD) determined gradients can also generate spurious numerical local minima. In addition, for the sonic boom problem, grid fineness will also determine the accuracy of the final design solution. Design optimization methods work well on problems defined by continuous objective functions. The sonic boom signature design problem is not necessarily defined by a continuous objective function. The signature can have a variety of shapes; i.e. from N-wave to multiple shocks. The far-field or ground signature may not transition continuously from one shape to another and hence, may exhibit discontinuous behavior. This is also a source of difficulty in using design optimization methods.In the following sections, several low boom and one reference aircraft configuration will be analyzed to predict their sonic boom signature characteristics. Modifications to some of these designs will also be presented to demonstrate the feasibility of using CFD as a design tool and to demonstrate the feasibility of designing shaped sonic boom signatures. Design modifications to some configurations will be presented to demonstrate the feasibility of achieving shaped signatures with reduced levels and not necessarily to represent realistic or aerodynamically efficient design modifications. Fuselage volume or camber are used as design variables in order to have a minimal effect on the primary wing aerodynamics. The paper will also seek to demonstrate whether a hybrid or ramped signature is feasible to achieve. For the low-boom configurations, the CFD predicted signatures will be compared qualitatively to their MLT design signatures.

Siclari, Michael J.↗

NASA NTRS1988

Lunar composition - A geophysical and petrological synthesis

Lunar compositional constraints are derived on the basis of geophysical data (in particular, the lunar seismic model as revised by Nakamura) and petrological arguments. Only in the case of extreme assumptions can critical aspects of bulk lunar composition be demonstrated to be equivalent to the present-day terrestrial mantle; specifically, the moon has an Mg number that is too low and an alumina abundance that is too high. Over a broad range of crustal densities the presence of a metallic core at least 150 km in radius is necessary to reconcile the upper mantle lunar seismic model with mass and moment of inertia constraints. A significant seismic discontinuity at 500 km depth may mark the lowest extent of early lunar differentiation, possibly representing a transition between highly fractionated upper mantle and less fractionated, perhaps even primordial, middle/lower mantle.

Mueller, Steve↗

NASA NTRS1979

A two phase Mach number description of the equilibrium flow of nitrogen in ducts

Some additional thermodynamic properties of the usual two-phase form which is linear in the moisture fraction are derived which are useful in the analysis of many kinds of duct flow. The method used is based on knowledge of the vapor pressure and Gibbs function as functions of temperature. With these, additional two-phase functions linear in moisture fraction are generated, which ultimately reveal that the squared ratio of mixture specific volume to mixture sound speed depends on liquid mass fraction and temperature in the same manner as do many weighted mean two-phase properties. This leads to a simple method of calculating two-phase Mach numbers for various duct flows. The matching of one- and two-phase flows at a saturated vapor point with discontinuous Mach number is also discussed.

Bursik, J. W.↗

NASA NTRS1991

Air-breathing hypersonic vehicle guidance and control studies; An integrated trajectory/control analysis methodology: Phase 1

A tool which generates optimal trajectory/control histories in an integrated manner is generically adapted to the treatment of single-stage-to-orbit air-breathing hypersonic vehicles. The methodology is implemented as a two point boundary value problem solution technique. Its use permits an assessment of an entire near-minimum-fuel trajectory and desired control strategy from takeoff to orbit while satisfying physically derived inequality constraints and while achieving efficient propulsive mode phasing. A simpler analysis strategy that partitions the trajectory into several boundary condition matched segments is also included to construct preliminary trajectory and control history representations with less computational burden than is required for the overall flight profile assessment. A demonstration was accomplished using a tabulated example (winged-cone accelerator) vehicle model that is combined with a newly developed multidimensional cubic spline data smoothing routine. A constrained near-fuel-optimal trajectory, imposing a dynamic pressure limit of 1000 psf, was developed from horizontal takeoff to 20,000 ft/sec relative air speed while aiming for a polar orbit. Previously unspecified propulsive discontinuities were located. Flight regimes demanding rapid attitude changes were identified, dictating control effector and closed-loop controller authority was ascertained after evaluating effector use for vehicle trim. Also, inadequacies in vehicle model representations and specific subsystem models with insufficient fidelity were determined based on unusual control characteristics and/or excessive sensitivity to uncertainty.

Hattis, Philip D.↗

NASA NTRSDate not supplied

Material Response Analysis of PICA-NuSil at the Hypersonic Materials Environmental Test System

Phenolic impregnated carbon ablator (PICA) is a thermal protection system (TPS) that gained heritage during the Stardust mission and is the baseline TPS material for expeditions to Mars (e.g., Mars Science Laboratory and Mars 2020).1 PICA is friable; therefore, a polysiloxane resin (NuSil CV-1144-0) is deposited onto the surface of flight hardware before assembly, test, and launch operations (ATLO) to mitigate particle shedding. Furthermore, the Mars Science Laboratory (MSL) and Mars 2020 heatshields were instrumented with an array of thermocouple plugs (MEDLI Integrated Sensor Plugs or MISP) designed to measure the in-depth thermal response of PICA during flight. Post-flight analysis of MISP data suggests that recession was lower than predicted, and NuSil is suspected as a reasonable cause. NuSil primarily consists of a copolymer blend of dimethyl and diphenyl-substituted polysiloxane resin. Four decades of research in the field of polymer-derived ceramics demonstrate that polysiloxane resins form oxidation-resistant silicon oxycarbide materials at pyrolysis temperatures as low as 800 °C. -2 Furthermore, silicon oxycarbide materials exhibit excellent thermal stability at relatively high-temperatures (< 1200 °C) but eventually decompose via carbothermal reduction. -3 Therefore, it is reasonable to suggest that NuSil forms an oxidation-resistant thermal barrier coating on the TPS surface and influences the material response during atmospheric entry. Evidence for this notion was observed during prior arc-jet testing at both the Aerodynamic Heating Facility (AHF) and the Panel Test Facility (PTF) at the NASA Ames Research Center. Articles of PICA were subjected to testing in airflow, and a temperature discontinuity was recorded at the surface of coupons coated with NuSil. A glassy coating was observed to form at the surface in shear flow configurations (wedge and panel) under moderate heat fluxes. Furthermore, state-of-the-art material response models have yet to account for the underlying physics and chemistry of the temperature discontinuity. To this end, a pathfinder test campaign was executed at the Hypersonic Materials Environmental Test System (HyMETS) to enhance the fundamental understanding of PICA-NuSil ablation phenomena and provide a preliminary data set to aid in the development of material response models.

PICA↗

NASA NTRS1994

Multifractal characterizations of nonstationarity and intermittency in geophysical fields: Observed, retrieved, or simulated

Geophysical data rarely show any smoothness at any scale, and this often makes comparison with theoretical model output difficult. However, highly fluctuating signals and fractal structures are typical of open dissipative systems with nonlinear dynamics, the focus of most geophysical research. High levels of variability are excited over a large range of scales by the combined actions of external forcing and internal instability. At very small scales we expect geophysical fields to be smooth, but these are rarely resolved with available instrumentation or simulation tools; nondifferentiable and even discontinuous models are therefore in order. We need methods of statistically analyzing geophysical data, whether measured in situ, remotely sensed or even generated by a computer model, that are adapted to these characteristics. An important preliminary task is to define statistically stationary features in generally nonstationary signals. We first discuss a simple criterion for stationarity in finite data streams that exhibit power law energy spectra and then, guided by developments in turbulence studies, we advocate the use of two ways of analyzing the scale dependence of statistical information: singular measures and qth order structure functions. In nonstationary situations, the approach based on singular measures seeks power law behavior in integrals over all possible scales of a nonnegative stationary field derived from the data, leading to a characterization of the intermittency in this (gradient-related) field. In contrast, the approach based on structure functions uses the signal itself, seeking power laws for the statistical moments of absolute increments over arbitrarily large scales, leading to a characterization of the prevailing nonstationarity in both quantitative and qualitative terms. We explain graphically, step by step, both multifractal statistics which are largely complementary to each other. The geometrical manifestations of nonstationarity and intermittency, 'roughness' and 'sparseness', respectively, are illustrated and the associated analytical (differentiability and continuity) properties are discussed. As an example, the two techniques are applied to a series of recent measurements of liquid water distributions inside marine stratocumulus decks; these are found to be multifractal over scales ranging from approximately 60 m to approximately 60 km. Finally, we define the 'mean multifractal plane' and show it to be a simple yet comprehensive tool with many applications including data intercomparison, (dynamical or stochastic) model and retrieval validations.

Davis, Anthony↗

NASA NTRS1985

The effect of the arbitrary level assignment of satellite cloud motion wind vectors on wind analyses in the pre-thunderstorm environment

The impact of satellite-derived cloud motion vectors on SESAME rawinsonde wind fields was studied in two separate cases. The effect of wind and moisture gradients on the arbitrary assignment of the satellite data is assessed to coordinate surfaces in a severe storm environment marked by strong vertical wind shear. Objective analyses of SESAME rawinsonde winds and combined winds are produced and differences between these two analyzed fields are used to make an assessment of coordinate level choice. It is shown that the standard method of arbitrarily assigning wind vectors to a low level coordinate surface yields systematic differences between the rawinsonde and combined wind analyses. Arbitrary assignment of cloud motions to the 0.9 sigma surface produces smaller differences than assignment to the 825 mb pressure surface. Systematic differences occur near moisture discontinuities and in regions of horizontal and vertical wind shears. The differences between the combined and SESAME wind fields are made smallest by vertically interpolating cloud motions to either a pressure or sigma surface.

Peslen, C. A.↗

NASA NTRS1986

The effect of the arbitrary level assignment of satellite cloud motion wind vectors on wind analyses in the pre-thunderstorm environment

Peslen, C. A.↗

NASA NTRSDate not supplied

Development of a Consistent GEOsat Cloud Property Dataset for the CERES Climate Data Record

Cloud properties are critical for understanding the Earth’s radiation budget and cloud feedbacks. At NASA Langley Research Center, the Satellite ClOud and Radiative Property retrieval System (SatCORPS) provides real-time and historical analyses of clouds derived from Geostationary satellite (GEOsat) data for weather and climate applications. For the Clouds and the Earth’s Radiant Energy System (CERES) program, the global constellation of GEOsats has been analyzed since 2000 to help characterize and account for the diurnal cycle of clouds and their radiative impacts in the CERES climate data record. Obtaining consistent cloud properties over the GEOsat data record during the CERES era is a major objective but a significant challenge considering the diversity of imaging capabilities deployed during that time. The GEOsat data analysis approach for the current CERES Edition-4 (Ed4) data products was focused on accuracy and consistency with MODIS by employing as much spectral information as possible from each satellite. However, the inconsistent use of spectral information across GEOsats led to marked discontinuities in the spatial and temporal record of cloud properties that had to be accounted for post facto in downstream CERES processing. This paper reports progress in developing a new GEOsat analysis system for the next CERES edition (Ed5) that has potential to improve cross-platform consistency and continuity. In this approach, the spectral channel complement is limited to just 3-channels during daytime, ~0.65 µm (VIS), ~3.9 µm (NIR), and ~10.8 µm (IR), common to nearly all of the satellites in the record. At night, a 2-channel approach is taken with the NIR and IR, and ~6.7 µm bands that includes a machine learning approach for optically thick cloud properties. A tradeoff is the potential for reduced accuracy particularly using data from the more advanced satellites that have more spectral channels (e.g. SEVIRI, AHI and ABI) that are known to help improve thin cirrus detection, cloud-aerosol discrimination and estimates in other difficult conditions that challenge cloud remote sensing. The new continuity approach is applied to one month of global GEOSat data for each year of the CERES record since 2000 and compared with the Ed4 GEO and MODIS cloud property time series in order to evaluate the level of improved consistency in the GEOsat record and to assess the accuracy impacts. Cloud fraction will also be assessed with CALIPSO data. Outstanding issues and challenges will be discussed. The results are expected to guide future work needed to develop a more robust GEOsat cloud data record for CERES.

CERES CDR↗

NASA NTRS2002

The Delta(dot) B = O Constraint vs. Minimization of Numerical Errors in MHD Simulations

The MHD equations are a system of non-strictly hyperbolic conservation laws. The non-convexity of the inviscid flux vector resulted in corresponding Jacobian matrices with undesirable properties. On the other hand, the MHD equations can be derived from basic principles in either conservative or non-conservative form. The non-conservative system has a better conditioned eigensystem. The Delta(dot)B = 0 constraint of the A4HD equations is only an initial condition constraint. One does not need the Delta(dot)B condition to close the MHD system. We formulate our new low dissipative high order scheme together with the Cargo & Gallice (1997) form of the MHD approximate Riemann solver in curvilinear grids for both versions of the MHD equations. A novel feature of our new method is that the well-conditioned eigen-decomposition of the non-conservative MHD equations is used to solve the conservative equations. This new feature of the method provides well-conditioned eigenvectors for the conservative formulation, so that correct wave speeds for discontinuities are assured. The justification for using the non-conservative eigen-decomposition to solve the conservative equations is that our scheme has a better control of the numerical error associated with the Delta(dot)B condition. Consequently, computing both forms of the equations with the same eigen-decomposition is almost equivalent. It will be shown that this approach, using the non-conservative eigensystem when solving the conservative equations, also works well in the context of standard shock-capturing schemes.

Yee, H. C.↗

NASA NTRS1988

Rocks, resolution, and the record at the terrestrial K/T boundary, eastern Montana and western North Dakota

Reconstructions of mass extinction events are based upon faunal patterns, reconstructed from numerical and diversity data ultimately derived from rocks. It follows that geological complexity must not be subsumed in the desire to establish patterns. This is exemplified at the Terrestrial Cretaceous-Tertiary (K/T) boundary in eastern Montana and western North Dakota, where there are represented all of the major indicators of the terrestrial K/T transition: dinosaurian and non-dinosaurian vertebrate faunas, pollen, a megaflora, iridium, and shocked quartz. It is the patterns of these indicators that shape ideas about the terrestrial K/T transition. In eastern Montana and western North Dakota, the K/T transition is represented lithostratigraphically by the Cretaceous Hell Creek Formation, and the Tertiary Tullock Formation. Both of these are the result of aggrading, meandering, fluvial systems, a fact that has important consequences for interpretations of fossils they contain. Direct consequences of the fluvial depositional environments are: facies are lenticular, interfingering, and laterally discontinuous; the occurrence of fossils in the Hell Creek and Tullock formations is facies-dependent; and the K/T sequence in eastern Montana and western North Dakota is incomplete, as indicated by repetitive erosional contacts and soil successions. The significance for faunal patterns of lenticular facies, facies-dependent preservation, and incompleteness is discussed. A project attempting to reconstruct vertebrate evolution in a reproducible manner in Hell Creek-type sediments must be based upon a reliable scale of correlations, given the lenticular nature of the deposits, and a recognition of the fact that disparate facies are not comparable in terms of either numbers of preserved vertebrates or depositional rates.

Fastovsky, D. E.↗

NASA NTRS2006

Simulations of Sea-Ice Dynamics Using the Material-Point Method

In recent years, the availability of large volumes of recorded ice motion derived from high-resolution SAR data has provided an amazingly detailed look at the deformation of the ice cover. The deformation is dominated by the appearance of linear kinematic features that have been associated with the presence of leads. These remarkable data put us in a position to begin detailed evaluation of current coupled mechanical and thermodynamic models of sea ice. This presentation will describe the material point method (MPM) for solving these model equations. MPM is a numerical method for continuum mechanics that combines the best aspects of Lagrangian and Eulerian discretizations. The material points provide a Lagrangian description of the ice that models convection naturally. Thus, properties such as ice thickness and compactness are computed in a Lagrangian frame and do not suffer from errors associated with Eulerian advection schemes, such as artificial diffusion, dispersion, or oscillations near discontinuities. This desirable property is illustrated by solving transport of ice in uniform, rotational and convergent velocity fields. Moreover, the ice geometry is represented by unconnected material points rather than a grid. This representation facilitates modeling the large deformations observed in the Arctic, as well as localized deformation along leads, and admits a sharp representation of the ice edge. MPM also easily allows the use of any ice constitutive model. The versatility of MPM is demonstrated by using two constitutive models for simulations of wind-driven ice. The first model is a standard viscous-plastic model with two thickness categories. The MPM solution to the viscous-plastic model agrees with previously published results using finite elements. The second model is a new elastic-decohesive model that explicitly represents leads. The model includes a mechanism to initiate leads, and to predict their orientation and width. The elastic-decohesion model can provide similar overall deformation as the viscous-plastic model; however, explicit regions of opening and shear are predicted. Furthermore, the efficiency of MPM with the elastic-decohesive model is competitive with the current best methods for sea ice dynamics. Simulations will also be presented for an area of the Beaufort Sea, where predictions can be validated against satellite observations of the Arctic.

Sulsky, D.↗