SEARCH · Engineering Papers

Results for “Surface Deformation”

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 235 records · Page 13

All sources NASA NTRS DOE OSTI

NASA NTRS1978

Local precision nets for monitoring movements of faults and large engineering structures

Along Bocono Fault were installed local high precision geodetic nets to observe the possible horizontal crustal deformations and movements. In the fault area there are few big structures which are also included in the mentioned investigation. In the near future, measurements shall be extended to other sites of Bocono Fault and also to the El Pilar Fault. In the same way and by similar methods high precision geodetic nets are applied in Venezuela to observe the behavior of big structures, as bridges and large dams and of earth surface deformations due to industrial activities.

Henneberg, H. G.↗

NASA NTRS2004

Status of a UAV SAR Designed for Repeat Pass Interferometry for Deformation Measurements

Under the NASA ESTO sponsored Instrument Incubator Program we have designed a lightweight, reconfigurable polarimetric L-band SAR designed for repeat pass deformation measurements of rapidly deforming surfaces of geophysical interest such as volcanoes or earthquakes. This radar will be installed on an unmanned airborne vehicle (UAV) or a lightweight, high-altitude, and long endurance platform such as the Proteus. After a study of suitable available platforms we selected the Proteus for initial development and testing of the system. We want to control the repeat track capability of the aircraft to be within a 10 m tube to support the repeat deformation capability. We conducted tests with the Proteus using real-time GPS with sub-meter accuracy to see if pilots could fly the aircraft within the desired tube. Our results show that pilots are unable to fly the aircraft with the desired accuracy and therefore an augmented autopilot will be required to meet these objectives. Based on the Proteus flying altitude of 13.7 km (45,000 ft), we are designing a fully polarimetric L-band radar with 80 MHz bandwidth and 16 km range swath. This radar will have an active electronic beam steering antenna to achieve Doppler centroid stability that is necessary for repeat-pass interferometry (RPI). This paper will present are design criteria, current design and expected science applications.

radar↗

NASA NTRS1999

Modeling the Volcanic Source at Long Valley, CA, Using a Genetic Algorithm Technique

In this project, we attempted to model the deformation pattern due to the magmatic source at Long Valley caldera using a real-value coded genetic algorithm (GA) inversion similar to that found in Michalewicz, 1992. The project has been both successful and rewarding. The genetic algorithm, coded in the C programming language, performs stable inversions over repeated trials, with varying initial and boundary conditions. The original model used a GA in which the geophysical information was coded into the fitness function through the computation of surface displacements for a Mogi point source in an elastic half-space. The program was designed to invert for a spherical magmatic source - its depth, horizontal location and volume - using the known surface deformations. It also included the capability of inverting for multiple sources.

Tiampo, Kristy F.↗

DOE OSTI2024

Enhanced joining strength in additive-manufactured polylactic-acid structures fused by embedded heated metallic meshes

Additively manufactured thermoplastic polymers, such as polylactic acid (PLA), hold significant promise for sustainable engineering structures, including wind turbine blades. Upscaling these structures beyond the limitations of 3D printer build volumes is a challenge; fusion joining presents a potential solution. This paper introduces a displacement-controlled resistance welding process for PLA, as an alternative to the typical force-controlled methods. Here, we investigated the bonding quality of resistance-welded and adhesive-bonded PLA beams through three-point bending and measured the surface deformations using digital image correlation. Different metal meshes (30%/0.11 mm Ni—Cu, 34%/0.07 mm Ni—Cu, and 36%/0.25 mm Co—Ni) served as heating elements. The process parameters were varied for the 34%/0.07 mm Ni—Cu mesh to identify an optimum set of parameters. Results showed that this optimized displacement-controlled welding achieved 94% of the original strength of monolithic samples. This indicates that the new welding process not only ensures high-quality bonding and fine surface finishing but also promotes sustainability, recyclability, and economic efficiency in various polymer and composite structural applications.

Additive manufacturing↗

NASA NTRS2000

Computations of Boiling in Microgravity

The absence (or reduction) of gravity, can lead to major changes in boiling heat transfer. On Earth, convection has a major effect on the heat distribution ahead of an evaporation front, and buoyancy determines the motion of the growing bubbles. In microgravity, convection and buoyancy are absent or greatly reduced and the dynamics of the growing vapor bubbles can change in a fundamental way. In particular, the lack of redistribution of heat can lead to a large superheat and explosive growth of bubbles once they form. While considerable efforts have been devoted to examining boiling experimentally, including the effect of microgravity, theoretical and computational work have been limited. Here, the growth of boiling bubbles is studied by direct numerical simulations where the flow field is fully resolved and the effects of inertia, viscosity, surface deformation, heat conduction and convection, as well as the phase change, are fully accounted for. Boiling involves both fluid flow and heat transfer and thus requires the solution of the Navier-Stokes and the energy equations. The numerical method is based on writing one set of governing transport equations which is valid in both the liquid and vapor phases. This local, single-field formulation incorporates the effect of the interface in the governing equations as source terms acting only at the interface. These sources account for surface tension and latent heat in the equations for conservation of momentum and energy as well as mass transfer across the interface due to phase change. The single-field formulation naturally incorporates the correct mass, momentum and energy balances across the interface. Integration of the conservation equations across the interface directly yields the jump conditions derived in the local instant formulation for two-phase systems. In the numerical implementation, the conservation equations for the whole computational domain (both vapor and liquid) are solved using a stationary grid and the phase boundary is followed by a moving unstructured two-dimensional grid. While two-dimensional simulations have been used for preliminary studies and to examine the resolution requirement, the focus is on fully three-dimensional simulations. The numerical methodology, including the parallelization and grid refinement strategy is discussed, and preliminary results shown. For buoyancy driven flow, the heat transfer is in good agreement with experimental correlations. The changes when gravity is turned off and/or fluid shear is added are discussed, as well as the difference between simulations of a layer freely releasing bubbles versus simulations using only one wavelength initial perturbation. Figure 1 shows the early stages of the formation of a three-dimensional bubble from a thin vapor layer. The boundary conditions are periodic in the x and y direction, the bottom is a hot and the top allows a free outflow. The jagged edge of the surface close to the bottom of the computational domain is due to some of the surface elements being on the other side of the domain and some elements not plotted by our plotting routine. In the second figure, we show the temperature distribution through two perpendicular planes.

Tryggvason, G.↗

NASA NTRS2022

Surface Biology and Geology (SBG) Thermal Infrared (TIR) Free -Flyer Concept

Scientists and engineers from the United States and Italy are developing a thermal infrared (TIR) free-flyer concept that forms an essential element of the NASA Surface Biology & Geology (SBG) designated observable and the greater NASA Earth System Observatory (ESO), which also includes observations of aerosols, clouds, convection, and precipitation; mass change; and surface deformation & change.In response to the recommendations from the 2017 National Academies decadal survey for Earth science, NASA initiated the SBG designated observable with five key research and applications focus areas: ecosystems and natural resources, hydrology, weather, climate, and solid Earth. SBG includes spaceborne measurements of hyperspectral imagery in the visible to shortwave infrared (0.4-2.5 um) and multispectral imagery in the mid and thermal infrared (3-12 um) that provide the remote sensing data needed to inform each research and applications area. High-level TIR data products include Earth surface temperature & emissivity, evapotranspiration, substrate composition, volcanic plumes, and high-temperature features. While the science, applications, and technology build on prior mission studies conducted for more than a decade, it has only been recently that a joint team of scientists and engineers from the NASA Jet Propulsion Laboratory (JPL), Agenzia Spaziale Italiana (ASI), Istituto Nazionale Geofisica e Volcanologia (INGV), and the Istituto Nazionale Astrofisica (INAF) have developed a compelling concept for a TIR free-flyer. In this concept, the TIR instrument is an eight-band radiometer. Seven of the bands are between 3 and 12 um with an additional band at 1.65 um. The instrument’s 68.8 deg total viewing angle and the free-flyer’s 665 km operational altitude results in a ground sampling distance (GSD) of <60 m at nadir, produces a swath width of 935 km, and allows for a 2-3-day revisit time. A two-band visible and near infrared (VNIR) camera, with <30 m GSD at nadir and 935 km swath width, will complement the TIR instrument. VNIR high-level data products include top & bottom of atmosphere reflectance and normalized difference vegetation indices (NDVI). The TIR instrument and VNIR camera will be mounted on a PLATiNO+ spacecraft, a multi-application, reconfigurable, and scalable platform. The launcher selected for the TIR free-flyer is the VEGA-C. To maximize the science and applications benefits the SBG-TIR free-flyer team are collaborating with personnel from the ESA Land Surface Temperature Monitoring mission (LSTM) as well individuals from the Thermal Infrared Imaging Satellite for High-resolution Natural resource Assessment (TRISHNA), a joint mission by CNES and ISRO. Pre-Phase A concept studies for the SBG-TIR free-flyer commenced in the spring 2021. This allows for the concept to be refined; for technical risks to be identified and mitigated, as necessary; and strategies for reducing the development schedule to be defined and assessed before conducting a NASA-ASI Mission Concept Review (MCR) in the spring of 2022 and a NASA Key Decision Point-A (KDP-A) as early as the summer of 2022. These are key milestones in being able to launch and begin operating the SBG-TIR free-flyer in the second half of the decade, so that many of the questions posed by the National Academies in the 2017 decadal survey can be answered.

Fabrizia-Buongiorno, Maria↗

NASA NTRS2006

Crustal Deformation and the Seismic Cycle across the Kodiak Islands, Alaska

The Kodiak Islands are located approx.120 to 250 km from the Alaska-Aleutian Trench and are within the southern extent of the 1964 Prince William Sound (M(sub W) = 9.2) earthquake rupture and aftershock zone. Here we report new campaign GPS results (1993-2001) from northeastern Kodiak and reprocessed GPS results (1993-1997) from southwestern Kodiak. The rate and orientation of the horizontal velocities, relative to a fixed North America, range from 29.7 +/- 1.7 mm/yr at N30.3degW +/- 3.3deg, located approx.120 km from the deepest point of the trench, to 8.0 +/- 1.3 mm/yr at N62.4degW +/- 9.3deg, located approx.230 km from the trench. We evaluated alternate models of coseismic and interseismic slip to test the importance of the mechanisms that account for surface deformation rates. Near the Gulf of Alaska coastal region of Kodiak the horizontal velocity can be accounted for primarily by the viscoelastic response to plate motion and a locked main thrust zone (MTZ), down-dip creep, and to a lesser extent, slip in the 1964 earthquake. Further inland the dominant mechanisms that account for post-1964 uplift rates are time-dependent, down-dip creep and a locked MTZ; for the horizontal velocity component southwest translation of western Kodiak may be important as well. Based on the pre-1964 and post-1964 earthquake pattern of interseismic earthquakes, we suggest that between the occurrences of great earthquakes like the 1964 event, more moderate to large earthquakes occur in the southwestern Kodiak region than near northeastern Kodiak .

Sauber, Jeanne↗

NASA NTRS1995

Coseismic fault slip associated with the 1992 M(sub w) 6.1 Joshua Tree, California, earthquake: Implications for the Joshua Tree-Landers earthquake sequence

Coseismic surface deformation associated with the M(sub w) 6.1, April 23, 1992, Joshua Tree earthquake is well represented by estimates of geodetic monument displacements at 20 locations independently derived from Global Positioning System and trilateration measurements. The rms signal to noise ratio for these inferred displacements is 1.8 with near-fault displacement estimates exceeding 40 mm. In order to determine the long-wavelength distribution of slip over the plane of rupture, a Tikhonov regularization operator is applied to these estimates which minimizes stress variability subject to purely right-lateral slip and zero surface slip constraints. The resulting slip distribution yields a geodetic moment estimate of 1.7 x 10(exp 18) N m with corresponding maximum slip around 0.8 m and compares well with independent and complementary information including seismic moment and source time function estimates and main shock and aftershock locations. From empirical Green's functions analyses, a rupture duration of 5 s is obtained which implies a rupture radius of 6-8 km. Most of the inferred slip lies to the north of the hypocenter, consistent with northward rupture propagation. Stress drop estimates are in the range of 2-4 MPa. In addition, predicted Coulomb stress increases correlate remarkably well with the distribution of aftershock hypocenters; most of the aftershocks occur in areas for which the mainshock rupture produced stress increases larger than about 0.1 MPa. In contrast, predicted stress changes are near zero at the hypocenter of the M(sub w) 7.3, June 28, 1992, Landers earthquake which nucleated about 20 km beyond the northernmost edge of the Joshua Tree rupture. Based on aftershock migrations and the predicted static stress field, we speculate that redistribution of Joshua Tree-induced stress perturbations played a role in the spatio-temporal development of the earth sequence culminating in the Landers event.

Bennett, Richard A.↗

NASA NTRS1987

The effect of overburden thickness on tension fracture patterns above an uplifting dome

Four experiments demonstrate that tension-fracture patterns above an uplifting dome depend on the thickness of the overburden layer being deformed. Four layers of increasing thicknesses (4.92 cm, 6.92 cm, 9.05 cm, and 11.12 cm) of a very fine sand (85%) and silt-clay (15%) mixture were updomed by slowly inflating a 1.22 m-diameter circular rubber pillow. The upper 2 cm of each layer was wetted and air dried to make it brittle and susceptible to fracture. The fractures that formed during these experiments exhibited a continuum of patterns from dominantly arcuate to dominantly radial as the overburden thickness increased. However, fracture density remained constant in each case for a given amount of surface deformation.

Pranger, H., II↗

NASA NTRS2001

Evaluation of Instability Phenomena in Sands: Plane Strain Versus Triaxial Conditions

Extensive research was carried out in the 1950s on theories of plasticity to extend the concepts developed for metals to materials that failed according to the Mohr-Coulomb criterion. The new ideas made it possible to merge the two distinct concepts (strength and deformation techniques) into one that relies on better understanding of plasticity and resulted in a rapid growth in the field of constitutive modeling of soil behavior. At the same time advanced experimental apparatuses and laboratory procedures were developed to calibrate the models. However, most laboratory experiments on granular materials are performed under Conventional Triaxial Conditions (CTC) for the purposes of evaluating constitutive behavior and stability properties, whereas most geotechnical field problems are closer to the Plane Strain (PS) condition. The triaxial tests performed in most laboratories comprise a simplification over in situ states and allow easier and robust experimentation. Most landslide problems, failure of soils beneath shallow and deep foundations, and failure of retaining structures, are cases that can generally be considered as plane strain. Strength and deformation characteristics of granular materials loaded in plane strain may be considerably different from those observed in CTC. Most studies on sands were limited to evaluating the constitutive behavior and in some cases extended to briefly describing the associated instability phenomena. This paper presents the results of a series of PS and CTC experiments performed on fine uniform silica sand known as F-75 Ottawa sand. Advanced analysis techniques were used to study the instability phenomena, which yielded very accurate measurements of shear bands occurrences and patterns. Destructive thin-sectioning technique along with monitoring the specimen surface deformation was used in the PS experiments and Computed Tomography (CT) was used to investigate the progress of primary and secondary shear bands in specimens subjected to CTC. Comparison between the two cases will be presented and discussed.

Alshibli, Khalid A.↗

NASA NTRS1992

A prototype automatic phase compensation module

The growing demands for high gain and accurate satellite communication systems will necessitate the utilization of large reflector systems. One area of concern of reflector based satellite communication is large scale surface deformations due to thermal effects. These distortions, when present, can degrade the performance of the reflector system appreciable. This performance degradation is manifested by a decrease in peak gain, and increase in sidelobe level, and pointing errors. It is essential to compensate for these distortion effects and to maintain the required system performance in the operating space environment. For this reason the development of a technique to offset the degradation effects is highly desirable. Currently, most research is direct at developing better material for the reflector. These materials have a lower coefficient of linear expansion thereby reducing the surface errors. Alternatively, one can minimize the distortion effects of these large scale errors by adaptive phased array compensation. Adaptive phased array techniques have been studied extensively at NASA and elsewhere. Presented in this paper is a prototype automatic phase compensation module designed and built at NASA Lewis Research Center which is the first stage of development for an adaptive array compensation module.

Terry, John D.↗

DOE OSTI2024

Interaction-enhanced nesting in spin-fermion and Fermi-Hubbard models

The spin-fermion (SF) model postulates that the dominant coupling between low-energy fermions in near critical metals is mediated by collective spin fluctuations (paramagnons) peaked at the Néel wave vector, Q N , connecting hot spots on opposite sides of the Fermi surface. It has been argued that strong correlations at hot spots lead to a Fermi surface deformation (FSD) featuring flat regions and increased nesting. This conjecture was confirmed in the perturbative self-consistent calculations when the paramagnon propagator dependence on momentum deviation from Q N is given by χ − 1 ∝ | Δ q | . Using diagrammatic Monte Carlo (diagMC) technique we show that such a dependence holds only at temperatures orders of magnitude smaller than any other energy scale in the problem, indicating that a different mechanism may be at play. Instead, we find that a χ − 1 ∝ | Δ q | 2 dependence yields a robust finite- T scenario for achieving FSD. To link phenomenological and microscopic descriptions, we applied the connected determinant diagMC method to the ( t − t ′ ) Hubbard model and found that at large U / t > 5.5 before the formation of electron and hole pockets (i) the FSD defined as a maximum of the spectral function is not very pronounced; instead, it is the lines of zeros of the renormalized dispersion relation that deforms toward nesting, and (ii) the static spin susceptibility is well described by χ − 1 ∝ | Δ q | 2 . Flat FS regions yield a nontrivial scenario for realizing a non-Fermi liquid state. Published by the American Physical Society 2024

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

NASA NTRSDate not supplied

NASA’s Earth System Observatory Formulation Overview

NASA’s Earth System Observatory (ESO) is an array of Earth-focused, interconnected satellite missions focused on five core study areas: Surface Biology and Geology, Mass Change, Aerosols, Surface Deformation and Change, and Clouds, Convection, and Precipitation. Observatory development follows recommendations from the National Academies of Sciences, Engineering and Medicine’s 2017 Earth Science Decadal Survey. Data gathered by ESO missions will provide actionable science to inform decisions related to climate change, disaster mitigation, and wildfires, improve real-time agricultural processes, and many other applications. Targeting launch dates in the late 2020s and early 2030s, each ESO satellite will deliver valuable information, but by working together as a single observatory system, their combined data and imagery will provide the global community with a 4D, holistic view of Earth, from bedrock to atmosphere. The ESO is also building on the legacy of international collaboration in Earth science, with initial participation and collaboration on these missions across space agency partners, including the Japanese Aerospace Exploration Agency (JAXA), Centre National D'Etudes Spatiales (CNES), Canadian Space Agency (CSA), Deutsches Zentrum für Luft- und Raumfahrt (DLR), and Agenzia Spaziale Italiana (ASI). This paper provides an overview of ESO, as well as an update on ESO missions currently in formulation – the Atmosphere Observing System, Mass Change, and Surface Biology and Geology elements – including current mission architectures, international partner collaboration, science community engagement, and applications efforts, as well as how NASA and its partners will make ESO data accessible to users all over the world.

NASA Earth Science↗

NASA NTRS1999

Computations of Boiling in Microgravity

Tryggvason, Gretar↗

NASA NTRSDate not supplied

NASA’s Earth System Observatory Formulation Progress

The 2017 Earth Science Decadal Survey by the National Academies of Science, Engineering, and Medicine identified five foundational observations, or Designated Observables (DO), to be implemented as missions by NASA. To address these five DO areas – Surface Biology and Geology; Mass Change; Aerosols; Cloud, Convection, and Precipitation; and Surface Deformation and Change – NASA is developing the Earth System Observatory (ESO) as an array of Earth-focused, interconnected satellite missions. The ESO will build on the strong history of international partnerships in NASA Earth science, with initial participation and collaboration on these missions across space agency partners including the Japanese Aerospace Exploration Agency (JAXA), Centre National D'Etudes Spatiales (CNES), Canadian Space Agency (CSA), Deutsches Zentrum für Luft- und Raumfahrt (DLR), and Agenzia Spaziale Italiana (ASI). Targeting launch dates in the late 2020s and early 2030s, the ESO will work as a single observatory, with each satellite in the ESO delivering its own valuable information, but taken together, the data and imagery will provide the global community with a multidimensional, holistic view of the Earth. In addition to furthering the world’s scientific knowledge, data gathered by ESO missions will provide actionable science to inform decisions related to climate change, disaster mitigation, and wildfires, improve real-time agricultural processes, and many other applications. This paper will provide an overview of the ESO development activities, as well as an update on the individual ESO missions currently in formulation – Atmosphere Observing System, Gravity Recovery and Climate Experiment-Continuity, and Surface Biology and Geology – including current mission architectures and concepts, international partner collaboration, science community engagement, and applications efforts, as well as how ESO data will be freely open and accessible to users all over the world.

NASA Earth Science↗

DOE OSTI2023

Helfrich-Hurault elastic instabilities driven by geometrical frustration

Here, the Helfrich-Hurault (HH) elastic instability is a well-known mechanism behind patterns that form as a result of strain upon liquid crystal systems with periodic ground states. In the HH model, layered structures undulate and buckle in response to local, geometric incompatibilities in order to maintain the preferred layer spacing. Classic HH systems include cholesteric liquid crystals under electromagnetic field distortions and smectic liquid crystals under mechanical strains, where both materials are confined between rigid substrates. However, richer phenomena are observed when undulation instabilities occur in the presence of deformable interfaces and variable boundary conditions. Understanding how the HH instability is affected by deformable surfaces is imperative for applying the instability to a broader range of materials. In this review, the HH mechanism is reexamined and special focus is given to how the boundary conditions influence the response of lamellar systems to geometrical frustration. Lamellar liquid crystals confined within a spherical shell geometry are used as the model system. Made possible by the relatively recent advances in microfluidics within the past 15 years, liquid crystal shells are composed entirely of fluid interfaces and have boundary conditions that can be dynamically controlled at will. Past and recent work that exemplifies how topological constraints, molecular anchoring conditions, and boundary curvature can trigger the HH mechanism in liquid crystals with periodic ground states is examined. The review ends by identifying similar phenomena across a wide variety of materials, both biological and synthetic. The fact that the HH mechanism is a generic and often overlooked response of periodic materials to geometrical frustration is highlighted.

42 ENGINEERING↗

NASA NTRS1994

A spectral formalism for computing three-dimensional deformations due to surface loads. 1: Theory

We outline a complete spectral formalism for computing high spatial resolution three-dimensional deformations arising from the surface mass loading of a spherically symmetric planet. The main advantages of the formalism are that all surface mass loads are always described using a consistent mathematical representation and that calculations of deformation fields for various spatial resolutions can be performed by simpley altering the spherical harmonic degree truncation level of the procedure. The latter may be important when incorporating improved observational constraints on a particular surface mass load, when considering potential errors in the computed field associated with mass loading having a spatial scale unresolved by the observational constraints, or when treating a number of global surface mass loads constrained with different spatial resolutions. The advantages do not extend to traditional 'Green's function' approaches which involve surface element discretizations of the global mass loads. Another advantage of the spectral formalism, over the Green's function approach, is that a posteriori analyses of the computed deformation fields are easily performed. In developing the spectral formalism, we consider specific cases where the Earth's mantle is assumed to respond as an elastic, slightly anelastic, or linear viscoelastic medium. In the case of an elastic or slightly anelastic mantle rheology the spectral response equations incorporate frequency dependent Love numbers. The formalism can therefore be used, for example, to compute the potentially resonant deformational response associated with the free core nutation and Chandler wobble eigenfunctions. For completeness, the spectral response equations include both body forces, as arise from the gravitational attraction of the Sun and the Moon, and surface mass loads. In either case, and for both elastic and anelastic mantle rheologies, we outline a pseudo-spectral technique for computing the ocean adjustment associated with the total gravitational perturbation induced by the external forcing. Three-dimensional deformations computed using the usual Love number approach are generally referenced to an origin at the center of mass of the undeformed planet. We derive a spectral technique for transforming the results to an origin located at the center of mass of the deformed planet.

Mitrovica, J. X.↗

DOE OSTI2025

Polarization Switching on the Open Surfaces of the Wurtzite Ferroelectric Nitrides: Ferroelectric Subsystems and Electrochemical Reactivity

Binary ferroelectric nitrides are promising materials for information technologies and power electronics. However, polarization switching in these materials is highly unusual. From the structural perspective, polarization reversal is associated with the change of the effective polarity at the surfaces and interfaces from N‐to‐M terminated, suggesting strong coupling between ferroelectric and chemical phenomena. Phenomenologically, macroscopic studies demonstrate the presence of complex time dependent phenomena including wake‐up. Here, in this study, the polarization switching using the multidimensional high‐resolution piezoresponse force microscopy (PFM) and spectroscopy is explored, detecting both the evolution of induced ferroelectric domain, electromechanical response, and surface deformation during first‐order reversal curve measurements. The presence of two weakly coupled ferroelectric subsystems are demonstrated and the bias‐induced electrochemical reactivity. The observed behaviors are very similar to the recent studies of other wurtzite system but additionally include electrochemical reactivity, suggesting the universality of these behaviors for the wurtzite binary ferroelectrics. These studies suggest potential of high‐resolution multimodal PFM spectroscopies to resolve complex coupled polarization dynamics in materials. Furthermore, these PFM based studies are fully consistent with the recent electron microscopy observations of the shark‐teeth like ferroelectric domains in nitrides. Hence, it is believed that these studies establish the universal phenomenological picture of polarization switching in binary wurtzite.

36 MATERIALS SCIENCE↗