Location Of Interface Research Articles

A High-Order Shifted Interface Method for Lagrangian Shock Hydrodynamics

We present a new method for two-material Lagrangian hydrodynamics, which combines the Shifted Interface Method (SIM) with a high-order Finite Element Method. Our approach relies on an exact (or sharp) material interface representation, that is, it uses the precise location of the material interface. The interface is represented by the zero level-set of a continuous high-order finite element function that moves with the material velocity. This strategy allows to evolve curved material interfaces inside curved elements. By reformulating the original interface problem over a surrogate (approximate) interface, located in proximity of the true interface, the SIM avoids cut cells and the associated problematic issues regarding implementation, numerical stability, and matrix conditioning. Accuracy is maintained by modifying the original interface conditions using Taylor expansions. We demonstrate the performance of the proposed algorithms on established numerical benchmarks in one, two and three dimensions.

The influence of the shock-to-reshock time on the Richtmyer–Meshkov instability in reshock

Experiments on the Richtmyer–Meshkov instability (RMI) in a dual driver vertical shock tube (DDVST) are described. An initially planar, stably stratified membraneless interface is formed by flowing air from above and sulfur hexafluoride from below the interface location using the method of Jones & Jacobs (Phys. Fluids, vol. 9, issue 1997, 1997, pp. 3078–3085). A random three-dimensional, multi-modal initial perturbation is imposed by vertically oscillating the gas column to produce Faraday waves. The DDVST design generates two shock waves, one originating above and one below the interface, with these shocks having independently controllable strengths and interface arrival times. The shock waves have nominal strengths of $M_L=1.17$ and $M_H=1.18$ for the shock wave originating in the light and heavy gas, respectively, with these strengths chosen to result in arrested bulk interface motion following reshock. The influence of the length of the shock-to-reshock time, as well as the order of shock arrival, on the post-reshock RMI is examined. The mixing layer width grows according to $h\propto t^\theta$ , where $\theta _H=0.36\pm 0.018$ (95 %) and $\theta _L=0.38\pm 0.02$ (95 %) for heavy and light shock first experiments, respectively, indicating no strong dependence on the order of shock wave arrival. Volume integrated specific turbulent kinetic energy (TKE) in the mixing layer versus time is found to decay according to $E_{tot}/\bar {\rho }\propto t^p$ with $p_H=-0.823\pm 0.06$ (95 %) and $p_L=-1.061\pm 0.032$ (95 %) for heavy and light shock first experiments, respectively. Notably, the 95 % confidence intervals do not overlap. Analysis on the influence of the shock-to-reshock time on turbulent length scales, transition criteria, spectra and mixing layer anisotropy are also presented.

Comprehensive Theoretical Formulation and Numerical Simulation of the Internal Flow in Pressure-Swirl Atomizers Type Screw-Conveyer

The present work shows the development of a comprehensive theoretical formulation for its application in the study of the internal flow of pressure-swirl atomizers with helical channels: “screw-conveyer”, which are characterized by presenting in their inlet channels, an angle of incidence or helix angle ψ. This angle originates a trigonometric factor (cosψ) that must be considered in the geometrical characteristics parameter of pressure-swirl atomizer (Ah), which consequently involves other geometric parameters, such as the annular section coefficient (φ), discharge coefficient (Cd), spray angle (2α), etc., being relevant in the internal flow study and design of the pressure-swirl atomizers type screw-conveyer. This theoretical formulation integrates an internal ideal flow model (Abramovich theory) with a model that considers the influence of the liquid viscosity (Kliachko theory) and hydraulic resistance of Idelchik. For the validation of this theoretical formulation, numerical simulation was used, considering the commercial software Ansys Fluent 2023 R2 furthermore, hexahedral meshes were generated with the ICEM CFD software 2023, for four cases of helix angle ψ (15°, 30°, 45° and 60°), with application of the RNG k-ε turbulence model and VOF multiphase model (volume of fluid) for the location of the liquid-gas interface and spray angle visualization.

A PNP ion channel deep learning solver with local neural network and finite element input data

Abstract This paper presents a deep learning method for solving an improved one-dimensional Poisson–Nernst–Planck ion channel (PNPic) model, called the PNPic deep learning solver. The solver combines a novel local neural network, adapted from the neural network with local converging inputs, with an efficient PNPic finite element solver, developed in this work. In particular, the local neural network is extended to handle the complexities of the PNPic model—a system of nonlinear convection–diffusion and elliptic equations with multiple subdomains connected by interface conditions. The PNPic finite element solver efficiently generates input and reference datasets for fast training the local neural network, as well as input datasets for quickly predicting PNPic solutions with high accuracy for a family of PNPic models. Initial numerical tests, involving perturbations of model parameters and interface locations, demonstrate that the PNPic deep learning solver can generate highly accurate numerical solutions.

Erosion of surfaces by trapped vortices

Two two-dimensional free boundary problems describing the erosion of solid surfaces by the flow of inviscid fluid in the presence of trapped vortices are considered. The first problem tackles an initially flat, infinite fluid-solid interface with uniform flow at infinity and a vortex in equilibrium above the surface. The second involves flow around a finite body with a trailing Föppl-type vortex pair. The conformal invariance of the complex potential permits both problems to be formulated as a Polubarinova–Galin (PG) type equation in which the time-dependent eroding surface in the physical z-plane is mapped to the fixed boundary of the ζ\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\zeta $$\\end{document}-disk. The Hamiltonian governing the equilibrium position of the vortex (or vortex pair in the second problem) is also found from the same map. In each problem, the PG equation giving the conformal map is found numerically and the time-dependent evolution of the interface and vortex location is determined. Different models governing the erosion of the interface are investigated in which the normal velocity of the boundary depends on some given function of the fluid flow velocity at the boundary. Typically, in the infinite surface case, erosion leads to the formation of a symmetric valley beneath the vortex which, in turn, moves downward toward the interface. A finite body undergoes erosion which is asymmetric in the flow direction leading to a flattening of the lee surface of the body so displaying some similarity to the experiments and associated viscous theory of Ristroph et al, Moore et al (Proc Natl Acad Sci 109(48):19606–19609, 2012, Phys Fluids 25(11):116602, 2013).

Characteristics and geological significance of carbon and oxygen isotopes of the Permian Lucaogou Formation dolomite in the southern Junggar Basin, northwestern China

Anomalously positive δ13C values in ancient dolomites are very rare. Dark gray argillaceous rocks of the lacustrine sediments of the Permian Lucaogou Formation are important source rocks in the Junggar Basin, and dolomites of varying thicknesses from 10 cm to 150 cm are often interspersed in argillaceous rocks. Based on the study of petrographic sections, this paper systematically analyzes the carbon and oxygen isotopes of dolomite and discusses the causes of abnormally high carbon isotope values and their significance in reconstructing paleoenvironment and paleoclimate. The results show that carbon isotope values are abnormally high in the dolomite of Lucaogou Formation, and the δ13C value is between +3.2 ‰ PDB and +19.6‰ PDB, with an average of +9.7‰ PDB. The δ18O values range from -17.4‰ PDB to -1.7‰ PDB, with an average of -8.1‰ PDB. From the lower part to the upper part of the Lucaogou Formation, the carbon isotope value gradually increases and becomes increasingly positive, and the carbon isotope of the dolomite deposited near the shore is more positive than that of the dolomite deposited far from the shore. The anomalously positive δ13C of the dolomite is mainly caused by microbial methanogenesis, with some contribution from evaporation. Microorganisms are mainly distributed at the redox interface. Evaporation controls the salinity and fluctuation of the redox interface in sedimentary water. The positive deviation difference in carbon isotopes between nearshore and offshore sedimentary dolomites may be related to the location of the redox interface during deposition. Together, the petrographic features and carbon and oxygen isotope signatures of the sections reflect the gradual evolution of the paleolake from a hydrologically open environment to a hydrologically closed one and the possible transition of the paleoclimate from a relatively warm to an arid condition, which is possibly a geochemical response to global climate change in the Permian period.

Dislocation proliferation at the growth crystal/seed interface of physical vapor transport-grown 4H-SiC crystals

4H-SiC single crystal with [000 1¯ ] direction 4° off axis was prepared by physical vapor transfer method. To enlarge the observed interface area between the seed crystal and the nascent crystal, the wafers were processed at an Angle of approximately 2° in the (000 1¯ ) plane. The longitudinal optical phonon-plasmon coupled (LOPC) mode was measured by a laser Raman spectrometer, and the location of the growth interface was determined by evaluating the free carrier concentration at the seed-crystal interface. The defect structure at the seed-crystal/newly grown crystal interface and in nearby regions was examined by high-resolution x-ray diffractometry. The change in stress at the interface was characterized by a stress tester; the results showed that a greater stress at the interface correlated to more proliferated dislocations. The dislocation morphologies in the interfacial regions of different seed grains etched by molten KOH were observed via microscopy.

Research on the effect of control size and location of human-machine interface display on drivers’ multi-resource demands

Interaction by touching a human–machine interaction (HMI) display has gradually evolved into the mainstream mode of human-vehicle interaction. Limiting the level of driver distraction caused by HMIs is an important concern. To explore the different positions and sizes of controls based on the multi-resource demands of the driver, a driving simulator experiment was conducted. The goal was to collect behavioral and subjective workload data from drivers performing secondary tasks under various conditions. Overall, the resources required by the driver to perform secondary tasks were relatively low on the left side of the HMI display and relatively high on the right side. Specifically, the lower-left position of the HMI display required the least resources, whereas the upper-right position required the most. Consequently, the results indicate that the resource demands of the secondary task on the driver decrease as the size of the controls increases. Most drivers preferred a “compromise” control size (18 mm). To reduce the resource demand, the control size of the HMI display should be set to at least 14 mm. This study has high industrial application value, and the results provide guidance for designing in-vehicle HMI displays with improved safety and usability.

Characterising Horizontal Two-Phase Flows Using Structured-Planar Laser-Induced Fluorescence (S-PLIF) Coupled With Simultaneous Two-Phase PIV (S2P-PIV)

Understanding the physics that underpins the behaviors of two-phase flows is immensely useful for a diverse number of industries, including oil and gas, nuclear, and aerospace. Stratified air-liquid two-phase flows are amenable to the application of two-phase PIV, but this is often only possible if PIV processing parameters can be chosen to match the separate gas and liquid phase image sets. Accurate interface location detection is, therefore, extremely important to facilitate appropriate masking during analysis. This paper examines the effectiveness of using the Structured-Planar Laser Induced Fluorescence (S-PLIF) method to aid in the accurate capture and processing of simultaneous two-phase PIV (S2P-PIV) for a gas-sheared liquid flow. Prior S2P-PIV work made use of direct interface detection, using contrast on images from the air-side camera. However, this method proved limiting in higher gas shear cases when the highly wavy nature of the flow prevented the camera from seeing all of the interface. An alternative prior approach using the liquid-side camera and the planar laser-induced fluorescence (PLIF) technique often led to large errors due to reflections from the free surface. Interface detection using the S-PLIF method facilitates a clearer delineation between the two phases through analysis of the change in the structure of the incident light, significantly reducing errors generated by the traditional PLIF. This then facilitates higher quality PIV analysis, both because interface identification is far more accurate and also because frame specific masking of the individual phases can be more automated. This paper presents air/water data from a rectangular channel with a liquid flowrate of 8 litres per minute and superficial gas velocities in the range 1 - 14 m/s. Data capture of the two phases is concurrent with the developed S-PLIF approach used to locate the interface, enabling separate and more accurate PIV analysis of the two phases over a wider parametric range than was previously possible. The data analysis shows that the transition to the disturbance/roll wave regimes is associated with a change in the liquid’s axial velocity profile and a significant increase in normalized axial velocity root-mean-square velocity fluctuation. Data sets from this work can be used to accurately describe high-speed gas-sheared two-phase flows that will be extremely useful for CFD model development and validation.

Experimental investigation on lateral resistance capacity enhancement effect of a self-centering shape factor in self-centering pier systems

In self-centering (SC) pier systems, a rocking interface is typically employed at the bottom of a pier to achieve a gap-opening mechanism. However, compared with the full-section bending state of the bottom section of a ductile pier (DP), the local compressive state of the rocking interface in the SC pier inevitably leads to a significant reduction in the lateral resistance capacity. Although a prestressing tendon (PT) force can improve the lateral resistance capacity, it also induces local damage. Moreover, its enhancement effect is limited, particularly in the case of a column with a high axial compression ratio. Therefore, this study proposes an innovative measure to significantly improve the hysteretic performance of an SC pier without inducing additional damage, which is to raise the location of the rocking interface. First, the SC shape factor ξSC is innovatively defined to quantify the above-mentioned measure, and its working mechanism is investigated in detail in this study. Second, a high-lateral-resistance SC pier system (HLRSCP), which applies the concept of ξSC, is proposed and introduced for its advantages. Third, hysteretic tests of the DP and HLRSCP with high axial compression ratios are performed and analyzed. Finally, the enhancement effect of ξSC on the lateral resistance capacity is parametrically examined using a finite element model, and the feasible region of ξSC is inferred and presented. The results validate the effectiveness and feasibility of the concept of the SC shape factor in improving the hysteretic performance of the SC pier system.

Spout geometry of fine particle conical spouted beds equipped with internal devices

The size and geometry of the spout are crucial parameters for the design and scaling-up of spouted beds because they condition heat and mass transfer rates. It is well-established in the literature that spouted beds without draft tube show axisymmetric spout shape, but this does not occur with open-sided draft tubes. Particle velocity has been measured at several bed heights and angular coordinates, under different geometric parameters and operating conditions, to determine the location of the spout-annulus interface in beds provided with open-sided draft tubes. Narrow draft tube wall slits and low air velocities lead to spouts of small size and triangular shape. The greatest spout expansion (high air velocities and wide draft tube wall slits) leads to a Reuleaux triangle shape, in which the spout swallows the tube. The results allow a better understanding of fine particle hydrodynamics in conical spouted beds, which is essential information for their scaling-up.

Visual coating inspection framework via self-labeling and multi-stage deep learning strategies

AbstractAn instantaneous and precise coating inspection method is imperative to mitigate the risk of flaws, defects, and discrepancies on coated surfaces. While many studies have demonstrated the effectiveness of automated visual inspection (AVI) approaches enhanced by computer vision and deep learning, critical challenges exist for practical applications in the manufacturing domain. Computer vision has proven to be inflexible, demanding sophisticated algorithms for diverse feature extraction. In deep learning, supervised approaches are constrained by the need for annotated datasets, whereas unsupervised methods often result in lower performance. Addressing these challenges, this paper proposes a novel deep learning-based automated visual inspection (AVI) framework designed to minimize the necessity for extensive feature engineering, programming, and manual data annotation in classifying fuel injection nozzles and discerning their coating interfaces from scratch. This proposed framework comprises six integral components: It begins by distinguishing between coated and uncoated nozzles through gray level co-occurrence matrix (GLCM)-based texture analysis and autoencoder (AE)-based classification. This is followed by cropping surface images from uncoated nozzles, and then building an AE model to estimate the coating interface locations on coated nozzles. The next step involves generating autonomously annotated datasets derived from these estimated coating interface locations. Subsequently, a convolutional neural network (CNN)-based detection model is trained to accurately localize the coating interface locations. The final component focuses on enhancing model performance and trustworthiness. This framework demonstrated over 95% accuracy in pinpointing the coating interfaces within the error range of ± 6 pixels and processed at a rate of 7.18 images per second. Additionally, explainable artificial intelligence (XAI) techniques such as t-distributed stochastic neighbor embedding (t-SNE) and the integrated gradient substantiated the reliability of the models.

Thermo-seismic performance of the novel thermal anchor pipe frame (TAPF) structure for supporting slope in earthquake-prone permafrost regions

The thermal anchor pipe frame (TAPF) structure is a new type of permafrost slope support technology, but its thermo-seismic characteristics are still unclear. In this paper, the coupled hydro-thermal-dynamic model of the thermal anchor pipe frame permafrost slope system is established. The cooling performance of the TAPF structure is investigated and the seismic responses are comparatively analyzed for different freeze-thaw periods. Results show that the presence of “elliptical” freezing cores with lower temperatures (which ranged from −3.94 to −5.36 °C) near the active layer at the end of the thawing period, which raised the permafrost table. The strip-shaped unfrozen moisture accumulation layer appears near the upper part of the thawing front at the highest temperature in the warm season. The peak acceleration response is the largest on October 15, 52.3% greater than the peak acceleration on April 15. The peak axial force of the thermal anchor pipe during the freezing period was 119 kN, which was 32.2% greater than that during the thawing period, and the peak total axial force after the earthquake appeared at the location of the freeze-thaw interface, and the axial force of the heat-transfer free section is greater than that of the anchored section. The frame structure is sensitive to earthquakes during the freezing period, and plastic bending damage may occur, especially at connection locations with thermal anchor pipes. Seasonal effects also have a strong influence on the shape of the acceleration response spectrum. The results can provide useful references for the maintenance and safe operation of permafrost slopes along the Qinghai-Tibet Project as well as for the seismic design of the new structure.

Experimental investigation on the effect of direct contact condensation regime on thermal stratification

Direct contact condensation affects the performance of a small modular reactor with pressure suppression system. Since the heat and mass transfer varies with the condensation regimes as well as steam-water interface shapes, the identification and classification of the condensation regimes are important. The experiment of steam direct contact condensation was carried out and the corresponding condensation regimes were studied. Specifically, the experiment was studied by injecting steam vertically into a wet well with subcooled water. The physical process of steam condensation phenomenon with varying undercooling degree of wet well was recorded by high-speed camera. According to the location of the steam-water interface and the shape of the steam bubble, the four condensation regimes were studied, namely chugging, external necking, spheroidal upward bubble and T-shaped upward bubble. Besides, it is found that thermal stratification occurs in the wet well at low steam mass flux, the influence of each condensation regime on the thermal stratification of wet well was studied. The experimental results show that the internal chugging and external necking cause a large downward circulation of subcooled water, which enhance the flow mixing around the discharge tube and the bottom of the wet well. Consequently, the thermal stratification phenomenon does not happen in those condensation regimes. However, the regimes like spheroidal upward bubble and T-shaped upward bubble cause the upward circulation of subcooled water in wet well, which cause the thermal stratification. The critical Richardson number (Ri) determines the condensation regimes. Furthermore, when the Ri is smaller than 1, the condensation regime is external necking. Besides, when the Ri is greater than 1, the condensation regimes are spheroidal upward bubble and T-shaped upward bubble.

Does backflow occur in forced imbibition into a dual-permeability pore network?

Backflow is a typical flow phenomena in two-phase displacement, but it is rarely studied in the forced imbibition and it remains unknown if the backflow occurs in forced imbibition into a porous medium. In this paper, we focus on the forced imbibition and investigate the backflow behavior in pore doublet (PD) models and in dual-permeability pore networks. Inspired by Gu et al. (2021) and Shan et al. (2023), a general mathematical model is first developed to predict the evolution of the penetration length of wetting fluid in a PD with arbitrary initial interface locations. Based on the analytical predictions and numerical simulations, we then study the influence of initial interface location, capillary number, viscosity ratio, and channel width ratio on the backflow. Results show that the backflow in a PD is mainly attributed to high flow rate in narrow channel, driven by the dominant capillary force. To supplement the insufficiency of inlet flow rate, the fluid in wide channel has to flow backwards to fill the narrow channel. The backflow strengthens with the decrease of capillary number or viscosity ratio, but it exhibits non-monotonic variation with initial interface location and channel width ratio. Moreover, the phase diagrams showing the critical conditions for the occurrence of backflow are established in terms of capillary number, viscosity ratio and channel width ratio. Finally, numerical simulations are conducted to explore the occurrence of backflow in three different dual-permeability pore networks. It is found that the backflow does not occur in dual-permeability pore networks due to the low flow rate caused by overlap event and pronounced tortuosity in the low permeability zone, or occurs only in the early stage when the capillary valve effect does not work.

Investigation of the Two-Channel Feedblock Zone in Co-Extrusion of Polymers

A process when different materials are combined to produce a product with multiple layers is called co-extrusion. During this process, polymers are melted in separate machines and then extrudate from different die channels. Once these channels converge, the polymers meet and flow through a single channel. The surface where the two fluids face is called “interface”. It is crucial to maintain the interface's uniformity and stability in order to achieve the desired multi-layered structure. Most of the issues in co-extrusion are related to issues that can be classified into two categories such as polymer encapsulation/interfacial distortion and die swell. To solve these problems, designers focus on improving the interface's stability. This paper examines effects of cross-section modification of the two-channel feedblock on the interface location and velocity and pressure distributions of the flow. The ANSYS software was used to simulate the co-extrusion of polymers, LLDPE and HDPE, in two-channel feedblock with rectangular, circular, and straight slot cross-sections. The results show that sharp corners increase the thickness of dead zones, while rounding them decreases the thickness. Additionally, stadium-shaped (or straight-slot) cross-section channels can move the flow with a higher maximum velocity and thinner boundary layer combining the results of rectangular and circular feedblocks.

Field Observations of Hydrostatic Pressure Deviations in a Nearshore Sediment Bed

AbstractMomentary liquefaction (ML), an instantaneous loss of effective stress between individual sand grains, is caused by the upward vertical pressure gradient of the excess pore‐water pressure in the sediment during the passage of a wave. Few field experiments have provided detailed observations of vertical pressure gradients within a sediment bed. Here, a novel autonomous pressure‐profiling instrument, the Pressure Stick, was deployed in the surf zone of an ocean beach in Rye, New Hampshire. The Pressure Stick has eight time‐synced absolute pressure sensors distributed over a 70 cm distance to measure absolute pressure in the water column and in the sediment bed. The absence of direct measurements of the sediment‐water interface location and sediment porosity limit direct observations of ML but allow for an examination of the hydrostatic pressure deviations. Observed upward vertical pressure head gradients within the sediment exceed the assumed effective weight of the sediment. Exceedances of two ML thresholds are observed within 16 cm of the sediment‐water interface and consistently occur at the onset of steeper and larger wave crests. Individual waves are characterized with an approximation of the Sleath parameter and the depth normalized wave height. The ML thresholds are exceeded 1%–40% of the time when the Sleath parameter is greater than 0.1 and 1%–33% of the time when the depth normalized wave height is greater than 0.4. These results suggest that wave steepness and nonlinearity may be necessary to induce large vertical pressure gradients that substantially deviate from hydrostatic pressure.

Soft Sensor Development for Real-Time Interface Tracking in Multiple Product Pipelines Based on Knowledge and Data

Summary Petroleum products are usually consecutively transported in the same multiple-product pipeline, and the occurrence of mixed oil happens during the process. Accurately tracking the mixed oil interface is essential for the optimal scheduling of oil transportation. However, complicated operating conditions and unavoidable measurement noise pose challenges to methods for locating the interface. The data-driven modeling method is a potential solution, but it may face limitations due to issues of overfitting, especially when the data set is contaminated with noise. To tackle such problems, in this paper, we propose a knowledge-informed Bayesian-Gaussian mixture regression (KIBGMR) model to enable the real-time tracking of the interface. The KIBGMR employs finite Gaussian distribution to learn the multimode characteristics of input data, including the hydrothermal data, measured density of tail oil, and the velocity of interface transportation and output data including the measured density of tail oil and interface transportation velocity. Subsequently, it utilizes the prior knowledge related to the regression coefficient through the Bayesian treatment. Evaluations demonstrate that the R2 index achieved by the proposed model in predicting the interface arrival time is greater than 0.98, even with the contaminated data set. This research can help operators accurately grasp the location of the mixed oil interface, formulate reasonable valve switch operations for better management of the mixed oil section, and provide a reference for the method of knowledge-data hybrid modeling.

Constrained Vapor Bubble Experiment (CVB) in the Light Microscopy Module (LMM)

Abstract This short article describes the major findings from the CVB experiment performed in the LMM on the International Space Station from 2010–2012. CVB was the first experiment to run in the new facility and focused on understanding the heat transfer and fluid mechanics occurring inside a wickless miniature heat pipe. The LMM was used to map the location of the vapor-liquid interface inside the device and to measure the film thickness profile on the walls of the device. Several interesting and unexpected phenomena were observed in microgravity including flooding of the heater end with liquid as the heat input increased, explosive nucleation of vapor bubbles at the heater end in the shortest version of the heat pipe tested, condensation on highly superheated surfaces, and the spontaneous formation of rip currents as the device tried to enhance the contact line area available for evaporation of the liquid.

Analytical verification of FFT-based micromechanical simulations near an elliptical crack tip by using composite voxels

In this paper, we make use of composite voxels in elasto-viscoplastic formulations based on fast Fourier transforms (EVPFFT) to calculate the stress field around an elliptical crack. The predicted stress field at the vicinity of the crack is verified against recently developed analytical solutions. The representation of an elliptical crack in the image-based method is stair-stepped, resulting in oscillations that are mitigated by the use of partially-filled, i.e. composite voxels. Dependencies of the oscillations in the predicted stress fields based on minor changes in the crack size crack size and aspect ratio are discussed and quantified. New understanding is gained of the factors that affect the oscillations, such as the aspect ratio of the ellipse, and the effect of minor changes in the crack interface location. The results show that the use of composite voxels significantly reduces the spurious oscillations in predicted stress fields, and consequently improves the agreement between the EVPFFT simulations and analytical solutions.