Process Zone Research Articles

Influence of surface integrity on short crack growth behavior in HFMI-treated welded joints

AbstractThis study investigates the influence of surface integrity and the localized fatigue phenomena on the initiation and propagation of short fatigue cracks in high-frequency mechanical impact (HFMI)-treated welded joints. The treated surface region, characterized by a compressive residual stress field, smooth notch geometry, and work hardening layer, improves welded joints’ fatigue strength. However, how these surface conditions influence the fatigue damage process zone during short crack initiation and growth is not yet well known. Therefore, this study systematically investigates the influence of different surface characteristics on fatigue life modeling of HFMI-treated welded joints made of high-strength steel. This is achieved using a non-local continuum damage mechanics-based approach of crack growth and elastic–plastic finite element simulation, explicitly modeling treated surface conditions. The simulated fatigue life is first verified with experiments and then applied to various surface conditions. The simulation results show that most of the fatigue life is spent until a crack size of 0.2 mm. The compressive residual stress field greatly extends both short crack initiation and propagation life, with its degree of contribution highly dependent on loading history and residual stress change. The role of the work hardening layer is mainly concentrated on improving fatigue life during short crack initiation and the very beginning of short crack growth.

Fracture process zone and fracture energy of heterogeneous soft materials

Bio-inspired heterogeneous soft materials are under rapid development due to their superior fracture and fatigue resistance. In the last few years, several kinds of fibrous soft composites in different length scales have been fabricated. However, the fracture behavior and toughening mechanism of this class of materials are still elusive. Here we develop a theoretical model for the crack tip field of fiber reinforced soft composites. The distribution of deformation around the crack tip and released elastic energy during crack propagation are obtained. The fracture process zone and fracture energy are quantified. There is a critical sample height, below which the fracture process zone size and fracture energy are size-dependent, above which they approach material-specific constants: steady-state fracture process zone size and steady-state fracture energy. A formula is derived to relate the steady-state fracture process zone size and parameters of the composite. It is found that both the steady-state fracture process zone size and the critical sample height scale with the fractocohesive length of the composite. The steady-state fracture energy of the composite can be enhanced either by enlarging the fracture process zone size through tuning fiber geometry or by increasing the work to rupture of the fiber through chemical treatment. This work reveals the toughening mechanism of heterogeneous soft materials and paves the way to design soft materials of high fracture energy, high fatigue threshold, and low hysteresis. It also provides a practical guideline for determining the sample size to measure the steady-state fracture energy of heterogeneous soft materials.

Investigation of fracture performance of asphalt mixtures considering size effect law and equivalent crack propagation law

The integration of the Size Effect Law (SEL) and the equivalent Linear Elastic Fracture Mechanics (LEFM) in characterizing the fracture behavior of asphalt mixtures has received limited attention. This study presents a preliminary exploration into the application of equivalent LEFM on the fracture characteristics of two asphalt mixtures, AC10 and AC16, accounting for variations in specimen size. Semi-circular Bending (SCB) tests were conducted on specimens with four different diameters (D = 76 mm, 100 mm, 125 mm, and 150 mm) and three notch-to-radius ratios (α = 0.2, 0.4, and 0.6) at a temperature of −10 °C. The equivalent crack propagation (ECP) length, crack growth resistance, and the size of the Fracture Process Zone (FPZ) were determined and subjected to analysis. The findings revealed a size effect on the ECP length, crack growth resistance, and FPZ size, with an increase in specimen size leading to a corresponding increase in these parameters. Conversely, an increase in α resulted in a decrease in all three types of data. A power law relationship was established between the ECP length and crack growth resistance. It was observed that AC16 demonstrated superior fracture resistance compared to AC10, particularly in terms of crack growth resistance. The initial notch length, ECP length, and FPZ size were found to be of the same order of magnitude, suggesting that the equivalent LEFM can produce more accurate fracture results under low temperatures than the traditional LEFM.

Interlaminar fracture toughness and impact resistance of carbon fiber reinforced composite with magnetic aligned CNTs

This study draws inspiration from Z-pins employed in fiber composites. Carbon fiber reinforced polymer (CFRP) laminates modified with Z-oriented carbon nanotubes (Z-CNTs) were fabricated utilizing a magnetic field orientation technique. The interlaminar fracture toughness, impact resistance and damage tolerance of CFRP modified with Z-CNTs were evaluated through end-notched flexure (ENF), low-velocity impact (LVI) tests, and compression after impact (CAI) testing. The enhancement mechanisms were investigated using multiscale techniques, including digital image correlation (DIC), scanning electron microscopy (SEM), molecular dynamics (MD), and finite element (FE) simulations. CFRP modified with 0.3 wt% CNTs showed optimum interlaminar fracture toughness, which improved by 62.9 % in CNT-random samples compared to non-modified samples and further enhanced by 95.2 % in CNT-oriented samples. The DIC and SEM analysis revealed the enhancement mechanism of Z-CNTs in fracture toughness, including (1) enlarging the toughness deformation of resin; (2) broadening the fracture process zone to dissipate the energy; and (3) shifting the failure mode of CNTs from pull-out dominance to fracture dominance. Furthermore, MD simulations showed that the crack propagation location significantly influences the different enhancement mechanisms exhibited by CNTs, which was accord with SEM observations. LVI and CAI tests indicated that specimens modified with Z-CNTs exhibit greater impact resistance and damage tolerance. FE simulation revealed that the mode II interlaminar fracture toughness plays a dominant role in the impact resistance of the laminate. These findings provide an effective strategy and theoretical supports for the design of impact resistance of composite laminates.

Tensile Fracture Initiation and Propagation of Granite and Gneiss at Wedge Splitting Tests: Part 2—Fracturing Studies Based on Digital Image Correlation and Thin Section Analysis

AbstractThe crack development in quasi-brittle granite and brittle gneiss under mode I loading condition was monitored using digital image correlation (DIC) technique during wedge splitting tests. The cracking behavior was studied on granite specimens that were split in one material direction, perpendicular to the rift plane, and gneiss specimens that were split in three different material directions, parallel and perpendicular to the foliation (along and across a lineation). The granite specimens had a saw cut 5 mm-wide notch and a blunt round notch in form of a 32 mm borehole. The gneiss specimens had a saw cut notch. The results from the DIC measurements revealed a meandering and branching crack path for the granite, whereas a smoother crack path for the gneiss with almost no branching. This behavior was confirmed by microscopy images from thin sections taken from specimens after testing. The thin sections showed that the fractures were prone to propagate across grains in the more coarse-grained granite than in the fine-grained gneiss, where the fractures propagated almost entirely along grain boundaries. The crack initiation occurred mainly in the corner of the saw cut notches and centrally in the round notch. However, the initiation locations in the granite were affected by the medium- to coarse-grained microstructure with grains preventing initiation and propagation which yielded displaced positions out from the ideal ones with respect to the highest stress in some cases. The crack opening displacement was determined along the crack path from the DIC measurements at 12 stages on each specimen of the advancing crack during the splitting progress. The critical crack opening displacement and length of the fracture process zone were assessed and the crack front position yielding the crack length along the tests was determined. The results showed a critical crack length when the deformations in the ligament (also called plastic hinge) affected the cracking process. The average crack velocity in gneiss during the test was more than twice as high as in the granite. This is attributed to a combined effect of the higher brittleness in the gneiss and the effect of a too large elastic energy in the specimen and test setup in relation to the dissipated fracture energy which made the initial crack propagation in the gneiss specimens nearly unstable. The strain energy release rate was calculated along the crack propagation and showed a lower value when the crack lengths were less than 40–60 mm. The calculation of the strain energy release rate was made on crack length measurements from DIC results. The results from the investigation were discussed in relation to the few other similar results found in the literature. The findings give an insight and understanding of the cracking process via both qualitative and quantitative results. Several used methods were novel or not used together in a single study as in this one.

Diffusive-length-scale adjustable phase field fracture model for large/small structures

In phase field fracture (PFF) method, the sharp crack is approximated by a phase field crack zone whose size is characterized by a diffusive length scale. Recently, the diffusive length scale is usually regarded as a constant material parameter determined by the fracture toughness, material strength, and young's modulus. As a result, the application of the PFF method poses challenges when dealing with structures whose sizes are much too large or small compared to the constant diffusive length scale. In details, for a large-scale structure, a significant computational burden is inevitable due to the limitation imposed by the constant diffusive crack length scale on the element size (i.e. the element size is generally at least smaller than half of the diffusive crack length scale to achieve the sufficient precision for the phase field process zone). For a small-scale structure, the crack patterns tend to be unclear and unrealistic due to the excessively large diffusive crack zone. To address these limitations, we propose a novel PFF method to make the relation between the diffusive length scale and the material parameters adjustable via modifying the energetic degradation functions. It is found that with the increase of the diffusive crack length scale, a transition from quasi-brittleness to brittleness is observed in the constitutive relationship curves, which coincides with the classical size effect on structural strength. Further, the Bažant's size effect in classical fracture mechanics can be reproduced by the present PFF method through scaling the size of a geometrically similar structure, i.e. a large-scale structure exhibits toughness-dominated fracture while a small-scale structure behavior strength-dominated fracture. The present PFF method efficiently addresses the mismatch of diffusive length scales for various material constituents by examining phase field crack growth in particle-reinforced composite plates. By adjusting the diffusive crack length scale based on structure size, it avoids high computational costs from large structures and unrealistic crack patterns from small ones. Moreover, it outperforms traditional PFF methods using a constant diffusive length scale in handling composite materials.

Multi-layer continuum model for adhesive FM 300-2 in end-notched flexure carbon laminate

The current investigation analyses the un-stable crack initiation of an adhesively bonded End Notched Flexure (ENF) specimen by means of a continuum damage formulation, based on the Bazant’s crack model, to reproduce accurately the load response and the fracture process zone of the adhesive. The specimen consisted of carbon fibre reinforced polymer adherends and two sheets of FM 300-2 adhesive films. The need of this implementation arose from experimental observations of normal transverse damage mechanisms at the bondline during the crack initiation. The damage model was based on a multi-layer configuration of the adhesive meso-structure that accounted for two damage activation functions in the normal and shear transverse directions. The results demonstrated the capability of the model to predict closely the force drop in the load–displacement curve of the ENF test. Different multi-layer configurations were studied following proper meso-structure simplification of the adhesive to use elastic and damage-capable layers in a correct order and thicknesses. The rise of the modelled shear strength allowed for keeping enough stored strain energy contributed to the low triaxiality stress in the bondline. Finally, the results demonstrated a significant influence of the layers modelling the knitted carrier to the local mechanisms of tensile fracture because of the normal transverse deformations. Unclassified. Based on Foreground information under EDA Contract No B.PRJ.RT.670 covering the Ad Hoc Project entitled “PATCHBOND II”.

Fracture process zone modelling of a magnesia spinel refractory using phase field fracture model

Fracture in quasi-brittle materials, such as refractories and reinforced concrete, involves complex mechanisms due to a progressive micro-cracking process within a fracture process zone (FPZ). This study employs Wu's phase field model (PFM) to simulate fracture behaviour in a magnesia spinel refractory. The PFM integrates fracture mechanics and damage mechanics, predicting tortuous crack patterns when heterogeneous strength distribution is considered. Numerical simulations, including wedge splitting tests typically applied for fracture testing of refractories, demonstrate PFM's effectiveness in capturing fracture behaviour, offering a robust tool for simulation of fracture of refractories with reduced brittleness. Comparative analysis with experimental data confirms the model's accuracy and applicability.

Strain‐Controlled Low‐Cycle Fatigue Behavior of a Superalloy Considering Extensive Temperature Range and Stress Concentration

ABSTRACTThis study aims to investigate the low cycle fatigue (LCF) behavior and fatigue life prediction method of an advanced solid solution–strengthened nickel‐based superalloy at various temperatures and stress concentrations. Cylindrical specimens as well as plates with three types of holes (including straight holes and inclined holes with different orientations) were designed and machined for fatigue tests. The test results revealed fatigue life decreases with the temperature increasing. Meanwhile, the fatigue lives of the plate specimens have decreased differently and are strongly dependent on the geometry and orientation of the holes. The key damage parameter related to the fatigue life was found on the strain gradient path. A criterion based on strain gradient is defined as the boundary of the damage process zone. The results show that the prediction was in good agreement with the experimental data, for the majority of the data are within a scatter band of ±2.

Investigating the effect of dynamic laser beam oscillations in remote fusion cutting process

The latest research on laser beam fusion cutting has revealed significant improvements in process productivity and cut quality through the use of dynamic beam shaping techniques. While many studies have investigated dynamic beam shaping for proximity cutting, the influence of laser beam oscillations on the remote fusion cutting process remains unexplored. The present work aims to study the effect of dynamic beam shaping on the remote fusion cutting process through analytical modeling, experimental investigations, and in situ high-speed monitoring. Initially, an analytical model based on thermodynamic analysis was developed to assess the influence of circular oscillations on the process zone. This model facilitates the evaluation of process performance from an energetic perspective, providing an estimate of the maximum achievable cutting speed for the remote fusion cutting process across various operating conditions. A significant increment in process productivity could be achieved through beam oscillations. Furthermore, based on theoretical findings, the effect of circular laser beam oscillations superimposed on the processing feed direction was experimentally investigated using a 1 mm thick AISI304 stainless steel material. A 6 kW fiber laser was utilized, alongside a high-speed camera-based system for in situ process monitoring. The experimental results demonstrate a significant increase in the process productivity under dynamic beam shaping conditions, consistent with theoretical findings. Specifically, the maximum achievable cutting speed could be increased from 0.13 to 0.20 m/s. Furthermore, the cut quality of produced samples was evaluated in terms of kerf morphology and profile.

Anisotropic Fracture Evolution and Size Effect of the Fracture Process Zone in Laminated Shales Derived from Digital Image Correlation

Anisotropic Fracture Evolution and Size Effect of the Fracture Process Zone in Laminated Shales Derived from Digital Image Correlation

Spatiotemporal Mapping of Auditory Onsets during Speech Production.

The human auditory cortex is organized according to the timing and spectral characteristics of speech sounds during speech perception. During listening, the posterior superior temporal gyrus is organized according to onset responses, which segment acoustic boundaries in speech, and sustained responses, which further process phonological content. When we speak, the auditory system is actively processing the sound of our own voice to detect and correct speech errors in real time. This manifests in neural recordings as suppression of auditory responses during speech production compared with perception, but whether this differentially affects the onset and sustained temporal profiles is not known. Here, we investigated this question using intracranial EEG recorded from seventeen pediatric, adolescent, and adult patients with medication-resistant epilepsy while they performed a reading/listening task. We identified onset and sustained responses to speech in the bilateral auditory cortex and observed a selective suppression of onset responses during speech production. We conclude that onset responses provide a temporal landmark during speech perception that is redundant with forward prediction during speech production and are therefore suppressed. Phonological feature tuning in these "onset suppression" electrodes remained stable between perception and production. Notably, auditory onset responses and phonological feature tuning were present in the posterior insula during both speech perception and production, suggesting an anatomically and functionally separate auditory processing zone that we believe to be involved in multisensory integration during speech perception and feedback control.

Finite element modeling for cohesive/adhesive failure of adhesive structures with a thermosetting resin

In this study, we established a novel method based on the finite element method (FEM) for predicting the strength of adhesive structures. Viscoplasticity, void growth, and cohesive zone model were introduced into the FEM to create a nonlinear damage growth (NDG) model. This model was used to comprehensively analyze the process zones within the adhesive layer. Furthermore, the embedded process zone approach was used to develop an interface constitutive law that averages the mechanical response of the adhesive layer. This modified Ma–Kishimoto (MMK) model can accurately represent the adhesive layer as an interface element and is computationally efficient. Furthermore, the study obtained the necessary interface properties for the MMK model from the NDG model, creating a numerical material test that can approximate the effect of the process zone. To validate the proposed method, single-lap shear tests were performed, and the accuracy of the predicted strength and deformation field was evaluated. The damage evolution in the NDG model and the MMK model were compared, and the scope of application of the MMK model was discussed. The results of this study can be used as a reference for the failure mechanism of thermosetting adhesives and establishment of design indices for adhesive structural strength.

Experimental study of type-I crack propagation in rock monitored by fiber Bragg grating

The occurrence of rock mass engineering disasters can be reduced significantly by obtaining the key information in the process of rock internal fracture accumulation and evolution through advanced monitoring means and taking certain preventive measures. Through fiber Bragg grating monitoring tests of the type-I crack propagation in rock, the spectral characteristics of the fiber Bragg grating in the process of rock failure were analyzed, the bandwidth expansion was extracted to characterize the local tensile strain during the crack propagation process, and the damage variable based on the bandwidth expansion was calculated to describe the type-I crack propagation process. The research results showed the following: (1) The fiber Bragg grating generated a non-uniform local tensile strain during type-I crack propagation. The bandwidth expansion of the wave spectrum was positively correlated with the local tensile strain. The change in the characteristics of the wave spectrum could provide the local real strain in the process of crack propagation. (2) During the steady propagation of the type-I crack, the bandwidth increased obviously, while the number of acoustic emission events was relatively small, and most of them were in the process zone. The gentle increase in the bandwidth corresponded to the unstable growth stage of the type-I crack, which could be used as a precursor of rock sample instability. During loading, the maximum tensile strain of the type-I crack was 6 %, while the maximum shear strain reached 1.8 %, primarily attributed to tensile failure. The average value of fracture toughness KIC, based on the equivalent linear elastic fracture mechanics model, was determined to be 1.21 MPa·m1/2, and the average value of fracture energy Gf was calculated as 57.7 N/m. The variation of peak reflectivity in fiber Bragg grating was associated with the closure of pre-existing defects and the propagation of microcracks in rocks, resulting in unstable shear strain at crack extension sites. (3) The damage variable D defined by the bandwidth broadening began to become concentrated and accumulated in the crack steady growth stage, and it increased slowly in the crack intensification and evolution stage. The grating was sensitive to changes when it was pasted to the surface at 90°, but it could easily fracture or become unable to capture spectral signals in time. Based on a comprehensive comparison, the monitoring effect was the best when the grating was pasted at 60° relative to the crack.

Dynamic tensile intralaminar fracture and continuum damage evolution of 2D woven composite laminates at high loading rate

The intralaminar tensile failure of 2D woven composites under dynamic tensile load was investigated in this paper. Compact tension samples were tested at high loading rate with an electromagnetic Hopkinson bar system. The strain field was obtained with high-speed imaging and digital image correlation, and the J-integral method was employed to obtain the fracture toughness and corresponding R-curve. The continuum damage evolution of intralaminar failure was then analyzed by tracking the opening near the initial crack tip. It is found that the dynamic intralaminar fracture toughness is decreased by 51% compared to the quasi-static condition, the continuum damage evolution and its dependence on loading rate have been reported as well. The failure mechanisms were studied with thermal imaging and scanned electron microscopy, shorter fibre pull-out length and thinner failure process zone may be responsible for the reduced toughness at high loading rate.

Strategic Procurement Practices and Performance of Export Processing Zone Firms in Machakos County, Kenya.

Export Processing Zones (EPZs) were established in Kenya to boost industrialization, attract foreign investment, and drive economic growth. However, despite significant government efforts, the performance of EPZ firms has been below expectations. This study investigates whether strategic procurement practices could improve the performance of EPZ firms in Machakos County, Kenya. Strategic procurement, which aligns procurement with organizational goals and emphasizes long-term value creation, was the focus. The study specifically explored the impact of strategic sourcing, strategic partnerships, strategic inventory management, and strategic contract management on the performance of EPZ firms in Machakos County. Using a descriptive research design, all 99 EPZ firms in Machakos County were surveyed through a census approach. Data was collected via structured questionnaires, with a pilot study conducted to ensure content validity and clarity. Both qualitative and quantitative data were analyzed, with results presented in tables and charts. The findings showed that the four strategic procurement practices adopted improved the performance of export processing zone (EPZ) firms in Machakos County. Strategic sourcing for instance significantly improved the operational efficiency and cost-effectiveness of EPZ firms. Firms that developed strategic partnerships with suppliers experienced enhanced product quality and supply chain reliability. Effective inventory management practices led to reduced waste and better resource allocation, positively impacting overall performance. Proper contract management was associated with improved compliance and reduced legal disputes, contributing to better performance metrics. The study concluded that adopting strategic procurement practices substantially enhances the performance of EPZ firms in Machakos County. By aligning procurement strategies with organizational goals and focusing on long-term value, firms can achieve greater efficiency, reliability, and cost savings. Policymakers should promote the adoption of strategic procurement practices within EPZs to boost economic performance.

Photodiode-based process monitoring for the ultrashort-pulsed laser structuring of the diffusion media for fuel cells

In recent years, the incorporation of process monitoring systems has become an integral part across several laser material processing applications. This is due to the feasibility of inline evaluations, which enable the elimination of downstream quality checks that are typically time-consuming and costly. Simultaneous to recent sensor developments, ultrashort-pulsed (USP) laser material processing has become a promising technology for creating complex microstructures in a wide variety of materials due to its high precision and negligible thermal input. So far, process monitoring has yet to be established in USP applications as it is unknown which sensors are suitable, which parameters can be observed, and which process-relevant information can be detected within the data. Within this work, a state-of-the-art photodiode-based process monitoring system was used to collect spectral data from the process zone and extract relevant information from the process. The laser structuring of the two-layered diffusion media used in polymer electrolyte membrane fuel cells was selected as the use case. This structuring process is challenging due to production-related fluctuations in the material thickness and the surface quality. Within this study, information was acquired from a sensor system employing three photodiodes. Through the application of explorative data analysis techniques and machine learning, AI models were created, focusing on determining the focus position of the laser beam and identifying the transition between the two layers. The developed models were capable of determining the initial focus position with varying accuracies. The best performance was achieved by a convolutional neural network using spectrograms as input data, while the worst accuracy was achieved by a ridge regression model using manually selected features. Furthermore, a clustering model was used, demonstrating the capability to detect the layer transitions within a specified tolerance.

Growth of thermoplastic shear band

Shear localization is a fundamental and widely observed non-equilibrium phenomenon in metallic materials subjected to impulsive loadings. Despite extensive research over several decades, the evolution of shear bands still remains poorly understood. In this study, we propose a two-stage model to describe the evolution of thermoplastic shear band, where the effects of strain hardening, strain-rate hardening, and thermal softening are involved. Our analysis reveals that the growth of thermoplastic shear band is derived by the competition between momentum diffusion and thermal diffusion processes, which is controlled by the solid-state equivalent of the Prandtl number. By incorporating these factors into our model, we find that the presence of hardening effect retards the evolution of shear bands, resulting in a wider shear band and increased dissipation of energy. Theoretical calculation results of shear bandwidth, shear band strain, propagation velocity, process zone geometry and dimensions, critical dissipation energy, shear band toughness, and shear band spacing are well consistent with available experimental data. Additionally, we explore the potential transition to turbulence controlled by Reynolds number within shear bands and investigate the influence of the Taylor-Quinney coefficient on the evolution of shear bands. These findings are expected to offer valuable insights into the mechanism of shear band evolution.

Mode II fracture properties of parallel neosinocalamus affinis bamboo strand lumber

This study investigates the Mode II fracture properties in the RL and TL planes of parallel neosinocalamus affinis bamboo strand lumber. Eighty End-Notched Flexure (ENF) specimens were designed, varying pre-crack lengths and specimen widths. Digital Image Correlation (DIC) was employed to analyze strain variations during the development of the fracture process zone (FPZ) and crack propagation stages. The Compliance Based Beam Method (CBBM) and equivalent crack length calculations were used to determine the strain energy release rate (GIIC) of PNABSL. R curves and crack growth rate curves revealed that crack propagation in the RL and TL planes. The results indicated that Mode II fracture propagation in both the RL and TL planes of PNABSL predominantly exhibited self-similar cracking, though the fracture surfaces were rough and uneven, with noticeable fiber bridging. The study found that Mode II fracture propagation was less stable in the RL plane compared to the TL plane, which is attributed to the easier fracturing of thick-walled fiber cells in the RL plane. The critical strain energy release rate for PNABSL was determined, offering valuable insights into its fracture behavior.

Experimental investigations on the formation of boundary layers in glass-based additive manufacturing of fused silica fibers by Laser Glass Deposition

Laser glass deposition is an additive manufacturing method to produce individualized glass components. This process uses a CO2 laser with a defocused beam as a heat source to additively melt-fused silica filaments. The fiber filament is fed laterally in the process zone and is deposited layer-by-layer to form 3D structures. The formation of boundary layers in conventional 3D-printing methods is a usual byproduct of the process. In this paper, the boundary layer formation of deposited fused silica filaments is investigated in detail by means of varying different process parameters such as, laser power, feed rate, laser spot diameter, and printing strategy. This involves examining both thin-walled and thick-walled test specimens. Quality characteristics like the surface roughness and the optical transmission are analyzed for the printed specimen. Finally, fully transparent structures with surface roughness below 100 nm and a transparency of 90% could be printed boundary layer-free and without post-processing.