Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Field stability”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 1,783 records · Page 99

Environmental Testing of SolarSCADA's SSI Data Acquisition System

In April of 2022, SolarSCADA, Inc., was awarded voucher funding to conduct research and analysis at the National Renewable Energy Laboratory (NREL) per the American-Made Solar Prize. SolarSCADA produces a fully integrated data acquisition system that is designed specifically for PV fields and includes a data monitoring platform. SolarSCADA is focused providing the SSI unit with all sensors/measurements per requirements in IEC-61724-1 [1] and with utility scale bankability. The objective of the SolarSCADA voucher funding was established as follows: Leverage NREL's existing calibration and environmental testing capabilities to evaluate the stability of the SSI channels and overall systeme in response to extreme environmental conditions associated with fielded PV.

14 SOLAR ENERGY↗

HERMeS Thruster Magnetic Field Topology Optimization Study: Plasma Plume Results

NASA’s Hall Effect Rocket with Magnetic Shielding (HERMeS) 12.5-kW Technology Demonstration Unit-1 (TDU-1) has been the subject of extensive technology maturation in preparation for flight system development. The magnetic optimization study aims to explore and investigate alternate magnetic field topologies. Four candidate magnetic field topologies that reduce the effectiveness of the shielding along the discharge channel walls with the intent to also reduce the erosion rates along the front pole covers were designed. Three of the four candidate magnetic field topologies (B1, B2, and B4) have been manufactured and were subjected to an extensive test campaign that included laser induced fluorescence (LIF), performance, stability, wear, and plume characterization. In Phase I, LIF measurements along the discharge chamber centerline found that upstream retraction of the thruster’s peak magnetic field does result in an upstream shift of the acceleration zone, but the magnitude of the shift does not correspond one-to-one to the shift in the location of the peak radial magnetic field magnitude. Phase II test segment results found that at a normalized thruster magnetic field setting of 1, the thruster performance and discharge current oscillations were very similar for configurations B0-B2. Also in Phase II abbreviated wear tests found that average inner front pole cover erosion rates for configuration B1 were approximately 65% of those at B0, and the average erosion rates for configuration B2 were approximately 35% of those at B0. This paper reports on the detailed plasma plume mapping results for configurations B0-B4. Faraday probe results indicate that, at 600 V operation, the plume divergence increased slightly for the B4 configuration when compared to B0-B2. Wien filter spectrometer results found that at 600 V, lower multiply-charged species fractions were attained as the thruster transitioned from configurations B0 to B4. Finally, retarding potential analyzer results found that at 600 V operation and at angles above 60,˚ significantly lower high-energy ion content was found for configurations B2 and B4 when compared to B0 and B1.

Electric propulsion,.Hall thruster, laser diagnost↗

Stability Analysis of the Flow Over a Swept Forward-Facing Step Using PIV Base Flows

Step excrescences are a type of surface imperfection encountered on swept wings of commercial aircraft, commonly because accessibility requirements prevent creating the wing’s surface from a single panel. If the step height is too large, super-critical, the flow can undergo an early transition to turbulence, which can be highly detrimental for the performance of the wing. In the case of a forward-facing step on a swept wing in a low-disturbance environment, stationary crossflow vortices can develop a significant amplitude upstream of the step and hence dominate the structure of the boundary-layer flow over the step. The main goal of the present investigation is to illuminate the path to transition supported by the fascinatingly complex flow field in the direct downstream vicinity of a step with a super critical height. The high-resolution, stereographic Particle Image Velocimetry (PIV) measurement dataset presently available for this flow field provides a complete description of the laminar flow for the execution of BiGlobal stability analysis in a plane parallel to the step. Although the notorious sensitivity of stability results to the description of the base flow demands a very careful uncertainty analysis of those results, it is argued that this very fact can be leveraged to produce new insight into the supported perturbation dynamics. In performing the analysis, several unsteady mode families are discovered that display the explosive perturbation expected for early transition to be induced. In considering domain widths equal to an integer-multiple of the incident crossflow-vortex wavelength and analyzing an extent of 5 crossflow-vortex wavelengths parallel to the step, it is found that the stability results converge while increasing the domain width. It is demonstrated, moreover, that the results for the wider domains can be approximated by appropriately averaging the results on neighboring single-crossflow-vortex-wavelength domains covering the same region. Besides being useful for computational purposes, this observed property suggests interpreting the instability mechanism as a distorted primary mechanism rather than a “proper” secondary mechanism. This follows in the context of the secondary in-stability analysis of three-dimensional boundary layers, because the secondary mechanism is usually characterized by being localized in a pocket of strong shear, while the distorted primary mechanism typically has an infinite support in the direction parallel to the step. Even though the growth rates are found to be sensitive to the interrogation-window size inherent to the PIV post-processing procedure, the spatial structure of the eigenfunctions is found to be relatively insensitive. Lastly, the spatial structure of the eigenfunctions corresponding to all velocity components are matched with the shape functions determined by computing the Spectral Proper Orthogonal Decomposition (SPOD) of a time-resolved measurement of the perturbation content.

forward-facing step↗

Evaluation of Mode Meters Robust to Forced Oscillations using Field-Measured Data

Mode meters are tools used by power system operators to continuously monitor a system's small-signal stability. They do so by estimating the system's electromechanical modes of oscillation. When a system undergoes a forced oscillation, mode meters may become biased because the two types of oscillation cannot be distinguished. Modified mode meter algorithms robust to this bias have been proposed in prior research, but these studies were based primarily on simulated data. In this paper, modified least squares and Yule-Walker mode meter algorithms are evaluated using field-measured data from phasor measurement units (PMUs). Results show that the sensitivities of the least squares algorithm make it impractical for use given the complexities of real-world forced oscillations. However, the modified Yule-Walker algorithm is shown to perform well and has significant potential for practical deployment in mode meter tools.

Follum, James D.↗

Out-of-plane chiral domain wall spin-structures in ultrathin in-plane magnets

Chiral spin textures in ultrathin films, such as skyrmions or chiral domain walls, are believed to offer large performance advantages in the development of novel spintronics technologies. While in-plane magnetized films have been studied extensively as media for current- and field-driven domain wall dynamics with applications in memory or logic devices, the stabilization of chiral spin textures in in-plane magnetized films has remained rare. Here we report a phase of spin structures in an in-plane magnetized ultrathin film system where out-of-plane spin orientations within domain walls are stable. Moreover, while domain walls in in-plane films are generally expected to be non-chiral, we show that right-handed spin rotations are strongly favoured in this system, due to the presence of the interfacial Dzyaloshinskii-Moriya interaction. These results constitute a platform to explore unconventional spin dynamics and topological phenomena that may enable high-performance in-plane spin-orbitronics devices.

36 MATERIALS SCIENCE↗

Surface Engineering by Simultaneous Action of Multiple External Fields (Final Technical Report)

This research project aims at enabling surface engineering strategies based on the fundamental understanding and predictive modeling of the surface morphological response of solid materials subjected to the combined action of multiple external forces. Surface morphological response to mismatch strain in strained-layer heterostructures also is analyzed. Our research focuses on the surface morphological evolution and stability of electrically conducting (metallic, such as Cu) and semiconducting (such as Si) solids under the simultaneous application of mechanical stresses, electric fields, and temperature gradients. We emphasize on identifying the conditions under which the multiply driven surface morphology is stable and on exploring the complexity of the corresponding various morphologically stable surface patterns. We also emphasize on the driven evolution of nanoscale features aiming at stabilizing and controlling nanoscale patterns on surfaces by their manipulation through simultaneously applied multiple external forces. Specifically, we address systematically the morphological response to the combined action of multiple external forces of bulk solid surfaces, surfaces of thin films grown epitaxially on thick or thin substrates, as well as nanoscale surface features such as coherently strained islands grown epitaxially on substrate surfaces. The research is based on a modeling approach that combines theoretical analyses of surface morphological stability with self-consistent dynamical simulations of surface evolution based on properly parameterized continuum and multiscale surface transport models.

36 MATERIALS SCIENCE↗

INSTABILITY OF MAGNETIZED IONIZATION FRONTS SURROUNDING H II REGIONS

An ionization front (IF) surrounding an H II region is a sharp interface where a cold neutral gas makes the transition to a warm ionized phase by absorbing UV photons from central stars. We investigate the instability of a plane-parallel D-type IF threaded by parallel magnetic fields, by neglecting the effects of recombination within the ionized gas. We find that weak D-type IFs always have the post-IF magnetosonic Mach number M{sub M2}≤1. For such fronts, magnetic fields increase the maximum propagation speed of the IFs, while reducing the expansion factor α by a factor of 1 + 1/(2β{sub 1}) compared to the unmagnetized case, with β{sub 1} denoting the plasma beta in the pre-IF region. IFs become unstable to distortional perturbations owing to gas expansion across the fronts, exactly analogous to the Darrieus-Landau instability of ablation fronts in terrestrial flames. The growth rate of the IF instability is proportional linearly to the perturbation wavenumber, as well as the upstream flow speed, and approximately to α{sup 1/2}. The IF instability is stabilized by gas compressibility and becomes completely quenched when the front is D-critical. The instability is also stabilized by magnetic pressure when the perturbations propagate in the direction perpendicular to the fields. When the perturbations propagate in the direction parallel to the fields, on the other hand, it is magnetic tension that reduces the growth rate, completely suppressing the instability when M{sub M2}{sup 2}<2/(2β{sub 1}−1). When the front experiences an acceleration, the IF instability cooperates with the Rayleigh-Taylor instability to make the front more unstable.

79 ASTRONOMY AND ASTROPHYSICS↗

Multi-extreme conditions at the Second Target Station

Three concepts for the application of multi-extreme conditions under in situ neutron scattering are described here. The first concept is a neutron diamond anvil cell made from a non-magnetic alloy. It is shrunk in size to fit existing magnets and future magnet designs and is designed for best pressure stability upon cooling. This will allow for maximum pressures above 10 GPa to be applied simultaneously with (steady-state) high magnetic field and (ultra-)low temperature. Additionally, an implementation of miniature coils for neutron diamond cells is presented for pulsed-field applications. The second concept presents a set-up for laser-heating a neutron diamond cell using a defocused CO 2 laser. Cell, anvil, and gasket stability will be achieved through stroboscopic measurements and maximum temperatures of 1500 K are anticipated at pressures to the megabar. The third concept presents a hybrid levitator to enable measurements of solids and liquids at temperatures in excess of 4000 K. This will be accomplished by a combination of bulk induction and surface laser heating and hyperbaric conditions to reduce evaporation rates. The potential for deployment of these multi-extreme environments within this first instrument suite of the Second Target Station is described with a special focus on VERDI, PIONEER, CENTAUR, and CHESS. Furthermore, considerations for deployment on future instruments, such as the one proposed as TITAN, are discussed. Overall, the development of these multi-extremes at the Second Target Station, but also beyond, will be highly advantageous for future experimentation and will give access to parameter space previously not possible for neutron scattering.

47 OTHER INSTRUMENTATION↗

Confined partial filament eruption and its reformation within a stable magnetic flux rope

We present observations of a confined partial eruption of a filament on 2012 August 4, which restores its initial shape within ≈2 hr after eruption. From the Global Oscillation Network Group Hα observations, we find that the filament plasma turns into dynamic motion at around 11:20 UT from the middle part of the filament toward the northwest direction with an average speed of ≈105 km s{sup –1}. A little brightening underneath the filament possibly shows the signature of low-altitude reconnection below the filament eruptive part. In Solar Dynamics Observatory/Atmospheric Imaging Assembly 171 Å images, we observe an activation of right-handed helically twisted magnetic flux rope that contains the filament material and confines it during its dynamical motion. The motion of cool filament plasma stops after traveling a distance of ≈215 Mm toward the northwest from the point of eruption. The plasma moves partly toward the right foot point of the flux rope, while most of the plasma returns after 12:20 UT toward the left foot point with an average speed of ≈60 km s{sup –1} to reform the filament within the same stable magnetic structure. On the basis of the filament internal fine structure and its position relative to the photospheric magnetic fields, we find filament chirality to be sinistral, while the activated enveloping flux rope shows a clear right-handed twist. Thus, this dynamic event is an apparent example of one-to-one correspondence between the filament chirality (sinistral) and the enveloping flux rope helicity (positive). From the coronal magnetic field decay index, n, calculation near the flux rope axis, it is evident that the whole filament axis lies within the domain of stability (i.e., n < 1), which provides the filament stability despite strong disturbances at its eastern foot point.

79 ASTRONOMY AND ASTROPHYSICS↗

HALL EFFECT CONTROLLED GAS DYNAMICS IN PROTOPLANETARY DISKS. II. FULL 3D SIMULATIONS TOWARD THE OUTER DISK

We perform three-dimensional stratified shearing-box magnetohydrodynamic (MHD) simulations on the gas dynamics of protoplanetary disks with a net vertical magnetic flux of B {sub z0}. All three nonideal MHD effects, Ohmic resistivity, the Hall effect, and ambipolar diffusion, are included in a self-consistent manner based on equilibrium chemistry. We focus on regions toward outer disk radii, from 5 to 60 AU, where Ohmic resistivity tends to become negligible, ambipolar diffusion dominates over an extended region across the disk height, and the Hall effect largely controls the dynamics near the disk midplane. We find that at around R = 5 AU the system launches a laminar or weakly turbulent magnetocentrifugal wind when the net vertical field B {sub z0} is not too weak. Moreover, the wind is able to achieve and maintain a configuration with reflection symmetry at the disk midplane. The case with anti-aligned field polarity (Ω⋅B{sub z0}<0) is more susceptible to the magnetorotational instability (MRI) when B {sub z0} decreases, leading to an outflow oscillating in radial directions and very inefficient angular momentum transport. At the outer disk around and beyond R = 30 AU, the system shows vigorous MRI turbulence in the surface layer due to far-UV ionization, which efficiently drives disk accretion. The Hall effect affects the stability of the midplane region to the MRI, leading to strong/weak Maxwell stress for aligned/anti-aligned field polarities. Nevertheless, the midplane region is only very weakly turbulent in both cases. Overall, the basic picture is analogous to the conventional layered accretion scenario applied to the outer disk. In addition, we find that the vertical magnetic flux is strongly concentrated into thin, azimuthally extended shells in most of our simulations beyond 15 AU, leading to enhanced radial density variations know as zonal flows. Theoretical implications and observational consequences are briefly discussed.

79 ASTRONOMY AND ASTROPHYSICS↗

A survey of the role of thermodynamic stability in viscous flow

The stability of near-equilibrium states has been studied as a branch of the general field of nonequilibrium thermodynamics. By treating steady viscous flow as an open thermodynamic system, nonequilibrium principles such as the condition of minimum entropy-production rate for steady, near-equilibrium processes can be used to generate flow distributions from variational analyses. Examples considered in this paper are steady heat conduction, channel flow, and unconstrained three-dimensional flow. The entropy-production-rate condition has also been used for hydrodynamic stability criteria, and calculations of the stability of a laminar wall jet support this interpretation.

Horne, W. C.↗

A subdivision-stabilized B-spline mixed material point method

Subjected to external loadings, polymeric materials, e.g., biological tissues, hydrogels, and elastomers, may undergo extreme, nearly incompressible, (self-)contact deformations. For numerical modeling employing mesh-based techniques such as the finite element method (FEM), these deformations pose significant challenges due to large distortions in the deformed geometry, accuracy issues stemming from volumetric locking effects, and increased computational cost from complex contact searches. As an alternative to mesh-based methods, the material point method (MPM), a continuum-based particle technique, is gaining attention for its ability to handle extreme distortions and capture no-slip contact without added cost. For nearly incompressible material behaviors, while mixed formulations can address locking effects by treating displacements and pressure as independent fields, they can suffer from numerical instabilities close to the incompressibility limit due to the violation of the inf-sup condition, leading to inaccurate nodal pressure solutions. Here we propose an efficient and stable mixed B-spline material point method with highest achievable regularity for quasi-compressible polymeric materials. Using the two-scale relation of B-splines, we introduce a subdivision-stabilization for the two-field mixed MPM and obtain numerically stable, oscillation-free nodal solutions with equal-order interpolations with optimal regularity. Building on the Eulerian-Lagrangian nature of MPM, a previously-converged solution framework is adopted to mitigate issues related to cell-crossing and numerical fracture artifact present in standard MPM. We assess the stability and accuracy of the developed mixed MPM at large deformations for soft materials through the benchmark Cook’s membrane problem. Additionally, we test the robustness of the proposed MPM by modeling several examples, including the compression and indentation of a circular block into a quasi-compressible substrate and the twisting deformation of a rectangular block. The findings demonstrate the MPM’s capabilities for modeling practical soft material applications.

36 MATERIALS SCIENCE↗

Vorticity, shocks, and magnetic fields in subsonic, ICM-like turbulence

We analyze data from high resolution simulations of the generation of compressible, MHD turbulence with properties chosen to resemble conditions in galaxy clusters. In particular, the flow is driven to have turbulence Mach number M{sub t}∼1/2 in an isothermal medium with an initially very weak, uniform seed magnetic field (β=P{sub g}/P{sub B}=10{sup 6}). Since cluster turbulence is likely to result from a mix of sheared (solenoidal) and compressive forcing processes, we examine the distinct turbulence properties for both cases. In one set of simulations velocity forcing is entirely solenoidal (∇⋅δu=0), while in the other it is entirely compressive (∇×δu=0). Both cases develop a mixture of solenoidal and compressive turbulent motions, since each generates the other. The development of compressive turbulent motions leads to shocks, even when the turbulence is solenoidally forced and subsonic. Shocks, in turn, produce and amplify vorticity, which is especially important in compressively forced turbulence. To clarify those processes we include a pair of appendices that look in detail at vorticity evolution in association with shocks. From our simulation analyses we find that magnetic fields amplified to near saturation levels in predominantly solenoidal turbulence can actually enhance vorticity on small scales by concentrating and stabilizing shear. The properties, evolution rates, and relative contributions of the kinetic and magnetic turbulent elements depend strongly on the character of the forcing. Specifically, shocks are stronger, but vorticity evolution and magnetic field amplification are slower and weaker when the turbulence is compressively forced. We identify a simple relation to estimate characteristic shock strengths in terms of the turbulence Mach number and the character of the forcing. Our results will be helpful in understanding flow motions in galaxy clusters.

79 ASTRONOMY AND ASTROPHYSICS↗

The Mu2e crystal calorimeter

The Mu2e experiment at Fermilab will search for the charged-lepton flavour violating conversion of negative muons into electrons in the coulomb field of an Al nucleus, planning to reach a single event sensitivity of about $3\cdot10^{−17}$, four orders of magnitude beyond the current best limit.The conversion electron has a clear monoenergetic signature at 104.967 MeV, slightly below the muon rest mass, and will be identified by a complementary measurement carried out by a high-resolution tracker and an electromagnetic calorimeter (EMC). The calorimeter is composed of 1348 pure CsI crystals, each read by two custom UV-extended SiPMs, that are arranged in two annular disks. The EMC should achieve better than 10% energy resolution and 500 ps timing resolution for 100 MeV electrons and maintain extremely high levels of reliability and stability and in a harsh operating environment with high vacuum, 1 T B-field and radiation exposures up to 100 krad and $10^{12} \rm{n}_\rm{1MeVeq}/\rm{cm}^{2}$.The calorimeter technological choice, along with the design of the custom front-end electronics, cooling and mechanical systems were validated through an electron beam test on a large-scale 51-crystals prototype (Module-0) and extensive test campaigns that characterised and verified the performance of crystals, photodetectors, analogue and digital electronics. This included hardware stress tests and irradiation campaigns with neutrons, protons, and photons.The production and QC phases of all calorimeter components is completed apart for the digital electronics that has been revised for improving its resilience to SEU and SEL. A series of vertical slice tests with the final electronics is being carried out on the Module-0 at LNF along with the implementation and validation of the calibration procedures. The first disk has been assembled in 2022 while the second disk is under assembly at the moment of writing. Status of construction will be summarised, along with plans for commissioning and first calibration of the fully assembled disks.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Magnetotail instability

The stability of the geomagnetic tail is investigated on the basis of three dimensional resistive magnetohydrodynamic simulations, using different dynamic constraints and different initial equilibria. Different forms of the energy equation for isotropic pressure are found to have no significant effect on the dynamic growth of a resistive tearing instability, which is responsible for near Earth reconnection, plasmoid formation and ejection, and the generation of fast plasma flows. The constraints of a modified double adiabatic approach, however, can quench the tearing instability through the development of large, mirror type, anisotropies in the boundary regions of the plasma sheet, unless isotropization occurs on fast, nearly Alfvenic, time scales. The presence of a net cross tail magnetic field component B(sub yN) can reduce the growth of the instability without complete stabilization. An increase of B(sub z) from midnight toward the tail flanks, however, by more than a factor of about 3, apparently completely stabilizes the tearing mode. Stabilization and destabilization thus may depend on properties and constraints (and their release) in regions other than the neutral sheet where reconnection is initiated.

Birn, Joachim↗

Controlling the stability of Pickering emulsions by pH-responsive nanoparticles

Electrostatic dissipative particle dynamics simulations were conducted to model the interactions between emulsion droplets stabilized by pH-sensitive polyelectrolyte-grafted nanoparticles. Using a steered molecular dynamics approach, a mechanistic study of forced coalescence was performed to probe the resistance between two particle-covered droplets. The degree of ionization of the grafted polyelectrolytes was adjusted to capture the pH responsiveness. The maximal resistance forces were measured to quantitatively discriminate the efficacy of particles in stabilizing emulsions at different degrees of ionization. Through analyzing droplet dynamics, resistance force variation, and electric field, we discovered that the resistance is attributed to direct electrostatic repulsion, the image charge effect near the water–oil interface, and steric hindrance among extended polymers. When the particle density on the droplet surface is relatively low, the increasing resistance forces at higher degrees of ionization can effectively prevent droplet coalescence. Oppositely, the ionization compromises emulsion stability when the particle surface coverage is high. Substantial desorption of particles from the interface was triggered as the degree of ionization increases. Finally, this in turn reduces resistance force and facilitates coalescence. Moreover, the nanoparticles prevent coalescence at high surface coverages by forming dense layers at individual interfaces, while the particle bridges straddling two interfaces were found at low surface coverages, which can also keep the droplets apart.

42 ENGINEERING↗

Affinity propagation clustering of full-field, high-spatial-dimensional measurements for robust output-only modal identification: A proof-of-concept study

Determination of the model order is a challenging problem in system identification, especially in output-only or operational modal identification where some modes are weakly excited. While existing methods such as the stabilization diagram method (spectral information) are effective, they do not scale to high-dimensional data, which is usually needed for high-fidelity characterization of structural dynamics and has been made available in the emerging full-field measurement techniques using optical methods such as photogrammetry and laser vibrometers. In this proof-of-concept study we present a new non-parametric, data-driven approach for robust output-only identification of high-spatial-dimensional modal parameters of basic structures by efficiently processing and interactively exploiting the full-field measurement (i.e., very dense spatial measurement points). Specifically, we first over-estimate the system model once, producing a pool of candidate modes associated with their modal frequencies and full-field, high-spatial-dimensional mode shapes. This is accomplished by a data-driven method termed affinity propagation clustering (APC), where the active clusters, which are the active modes in our formulations, emerge from the "message-passing" procedure and does not require a pre-determination of the cluster number (mode or model order). Next, rather than using the spectral information to distinguish the physical and spurious modes in the stabilization diagram method, we exploit and visualize the spatial, full-field mode shape associated with each candidate mode to do so. We conduct extensive experiments on basic structural models with comparisons to a few existing methods. The results indicate that the new method is computationally efficient for identifying high-spatial-dimensional modal parameters, and robust to identify weak modes by exploiting the full-field measurement. We also discuss its applicability and limitations for structures with complex geometry (shapes).

42 ENGINEERING↗

Potential Flow Interactions With Directional Solidification

The effect of convective melt motion on the growth of morphological instabilities in crystal growth has been the focus of many studies in the past decade. While most of the efforts have been directed towards investigating the linear stability aspects, relatively little attention has been devoted to experimental and numerical studies. In a pure morphological case, when there is no flow, morphological changes in the solid-liquid interface are governed by heat conduction and solute distribution. Under the influence of a convective motion, both heat and solute are redistributed, thereby affecting the intrinsic morphological phenomenon. The overall effect of the convective motion could be either stabilizing or destabilizing. Recent investigations have predicted stabilization by a flow parallel to the interface. In the case of non-parallel flows, e.g., stagnation point flow, Brattkus and Davis have found a new flow-induced morphological instability that occurs at long wavelengths and also consists of waves propagating against the flow. Other studies have addressed the nonlinear aspects (Konstantinos and Brown, Wollkind and Segel)). In contrast to the earlier studies, our present investigation focuses on the effects of the potential flow fields typically encountered in Hele-Shaw cells. Such a Hele-Shaw cell can simulate a gravity-free environment in the sense that buoyancy-driven convection is largely suppressed, and hence negligible. Our interest lies both in analyzing the linear stability of the solidification process in the presence of potential flow fields, as well as in performing high-accuracy nonlinear simulations. Linear stability analysis can be performed for the flow configuration mentioned above. It is observed that a parallel potential flow is stabilizing and gives rise to waves traveling downstream. We have built a highly accurate numerical scheme which is validated at small amplitudes by comparing with the analytically predicted results for the pure morphological case. We have been able to observe nonlinear effects at larger times. Preliminary results for the case when flow is imposed also provide good validation at small amplitudes.

Buddhavarapu, Sudhir S.↗