Cuboid Geometry Research Articles

Effective method for photonic jet shaping in reflection mode

We report an effective method to shape a photonic jet (PJ) generated by a dielectric cuboid scatterer on a hollow reflection screen. The study focuses on the shaping of PJ by hollow and cuboid geometries, including side length, depth, and position. The results show that all the geometric parameters can effectively shape the PJ in characteristic parameters of intensity, focal length (FL), and lateral size. The comparison shows that the present shaping method is much superior to the previously reported one based on a change of cuboid height, characterized by three times larger of the varying ranges of the PJ characteristic parameters. The shaping by the side length and depth of hollow and cuboid is explained from viewpoints of ray optics and electromagnetic field theory in terms of hollow-induced phase singularity in Poynting vector distribution. The hollow position mainly shapes the propagation trajectory of PJ. The deviation of the hollow position results in the curvature of the propagation trajectory of PJ, and the curvature is correlated with the asymmetric extent of the scattering system.

In-situ Selective Laser Heat Treatment for microstructural control of additively manufactured Ti-6Al-4V

As-built Laser Powder Bed Fusion (LPBF) Ti-6Al-4V typically exhibits a fully acicular α′-martensite microstructure, and requires post-process heat treatment in order to decompose the martensite and achieve sufficient ductility. In the present study, we demonstrate a simple concept based on in-situ Selective Laser Heat Treatment (SLHT) that can effectively alter the microstructure and activate the decomposition of the α′-martensite into a lamellar α+β microstructure within a short time scale (~30s). SLHT consists of multiple rescanning of the printed part, with low energy density, triggering solid-state phase transformations. Operando X-ray diffraction has been performed on cuboid and thin wall geometries, and was augmented by thermal finite element simulations. Upon SLHT, a gradual formation of the β phase as well as an α/α′peak narrowing trend have been evidenced through X-ray diffraction, as an indication of the diffusional nature of α′-martensite decomposition. Moreover, through fine tuning of the process parameters at the final stage of SLHT, a controlled temperature evolution during cooling was achieved, leading to preservation of the β phase, a product of the decomposition, down to room temperature. Complementary microstructural characterizations via EBSD, SEM, and TEM confirm the presence of a lamellar α+β microstructure after SLHT. Our results evidence, for the first time, the fast kinetics of α′-martensite decomposition under in-situ SLHT. The approach is meant to be implemented at selected locations during the LPBF process, avoiding time-consuming post processing steps, and leading to composite-like, architected microstructures.

ARCHER - a new Three-Dimensional method of characteristics neutron transport code for Pebble-bed HTR with coarse mesh finite difference acceleration

A new deterministic code ARCHER using Three-Dimensional (3-D) Method of Characteristics (MOC) has been developed by Institute of Nuclear and new Energy Technology (INET), Tsinghua University, aiming to obtain high-fidelity transport solution of pebble-bed High Temperature gas-cooled Reactor (HTR) with explicit pebble-bed geometry. However, 3-D MOC usually suffers from the slow convergence rate without efficient acceleration methods. In this work, the Coarse Mesh Finite Difference (CMFD) based on regular mesh is implemented to accelerate 3-D MOC calculation, which is available for the problems with cylinder and cuboid geometry. The complex geometry of reactor core with randomly stacked spherical fuel elements and irregular helium flow area, as well as the topological relationships between fine Flat Source Region (FSR) mesh for MOC and coarse CMFD mesh, are the key but challenging tasks for high-fidelity transport solution. Constructive Solid Geometry (CSG) and Boolean operation are constructed to describe problem geometry, which has the ability to deal with complex topological relationships between meshes. In addition, tree grid structure is utilized to reduce the time complexity of ray tracing and obtain the grid mapping relationship between coarse mesh and fine mesh. Spatial domain decomposition parallel based on MPI is utilized to alleviate the memory pressure and pursue high computational performance. Ray parallel is further achieved through OpenMP with shared memory. Several problems have been used to verify the geometric modeling capability and computational performance of this new 3-D MOC code for pebble-bed HTR.

Chaotic cavity design of a UV-C disinfection chamber for uniform radiation distribution.

UV-C radiation has been used extensively for disinfection wherein the dosage of radiation varies with the microorganism concerned. The efficacy of disinfection largely relies on the radiation reaching uniformly across the chamber without simultaneously leaving out dark regions within the chamber. The regular geometries of the chambers, such as with circular or rectangular cross sections, invariably lead to hot spots accompanied by regions that remain dark, thereby severely limiting the germicidal effects of the UV-C radiation. We propose and demonstrate the use of a chaotic Bunimovich stadium design chamber that ensures uniform UV-C exposure without any accompanying "dark spots." The design incorporates the shape parameters of a chaotic Bunimovich stadium, wherein ray-tracing simulations confirm the space-filling nature of the photon trajectories, typical of chaotic dynamics. This leads to uniform radiation distribution across the whole volume. Experimental results are presented for a biosafety cabinet that can be used for UV-C disinfection wherein the UV-C light intensity distribution is more evenly distributed than the conventional regular (cuboidal) geometry. Adaptation of this simple design consideration is paramount for air disinfection, wherein dark regions in conventional air ducts can severely compromise the efficacy of the disinfection.

Optimization of detector placements in reduction of multiple counts for μSR measurements at China Spallation Neutron Source

An Experimental Muon Source (EMuS) has been proposed to conduct muon spin rotation/relaxation/resonance (μSR) measurements at China Spallation Neutron Source (CSNS). To make better use of muons in each pulse, a highly segmented μSR spectrometer with more than 2000 detector channels is under design. Due to such high granularity of detectors, multiple counting events generated from particle scattering or spiral motion of positrons in a strong longitudinal field should be carefully considered in the design. According to the simulation, long scintillators have a good capability of angular discrimination. Detectors with cuboid geometries are better than those with frustum shapes. The cuboid detector with a length of 50 mm is longer enough to get the optimal range of discrimination angle. In a real μSR spectrometer, detectors can be placed parallelly along the beam direction or pointing to the sample. A figure of merit (FoM) has been proposed to compare such two arrangements by integrating their impacts on multiple counts and total counting loss in zero and longitudinal fields. The outstanding performance of multiple counting rejection due to the angular discrimination capability makes the pointing arrangement achieve much higher FoM. The simulation results can provide good support for the design of the highly segmented μSR spectrometer.

MagTense: A micromagnetic framework using the analytical demagnetization tensor

We present the open source micromagnetic framework, MagTense, which utilizes a novel discretization approach of rectangular cuboid or tetrahedron geometry “tiles” to analytically calculate the demagnetization field. Each tile is assumed to be uniformly magnetized, and from this assumption only, the demagnetization field can be analytically calculated. Using this novel approach we calculate the solution to the μmag standard micromagnetic problems 2, 3 and 4 and find that the MagTense framework accurately predicts the solution to each of these. Finally, we show that simulation time can be significantly improved by performing the dense demagnetization tensor matrix multiplications using NVIDIA CUDA.

Study of two-layer tapered depth of interaction PET detector

Improving detection efficiency in small animal PET scanners without degrading spatial resolution is one of the main problems of these scanners. Commercial small animal PET scanners use different methods to achieve desirable levels of sensitivity and spatial resolution. GE Healthcare eXplore VISTA PET scanner uses double layer (LYSO-GSO) depth-of-interaction (DOI) capable cuboid detector modules. In this work, the design of GE Healthcare eXplore VISTA PET scanner is improved using tapered detector geometry instead of cuboid geometry. Using tapered detector geometry, the gaps between adjacent modules are filled and the sensitive volume has increased about 11.5%. The new designed PET scanner sensitivity and spatial resolution are studied for different crystal layer configurations (LYSO-GSO and GSO-LYSO with different thicknesses). As expected, average sensitivity over FOV is improved. Spatial resolution is slightly degraded but it is still uniform over FOV.

Numerical and experimental studies of the hypersonic flow around a cube at incidence

In order to improve predictions of the on-ground casualty risk associated with the uncontrolled atmospheric re-entry of satellites from Low Earth Orbit, there is significant research interest in the development of engineering models of hypersonic heating rates to faceted shapes. A key part of developing such models is generating accurate datasets of the heat fluxes experienced by faceted shapes at various orientations in hypersonic flows. In this work, we use wind tunnel experiments and CFD simulations to study the hypersonic flow around a cuboid geometry at 5° incidence in a Mach 5 flow at Reynolds numbers of 79.5×103, 109×103 and 148×103. The wind tunnel data are obtained in the University of Manchester's High SuperSonic Tunnel and consist of schlieren images and temperature histories collected using infrared thermography. These temperature histories are used to calculate experimental heat fluxes by solving a three-dimensional inverse heat conduction problem. CFD simulations around the same geometry at equivalent free-stream conditions are calculated with the DLR-TAU code. The experimental and CFD results show good agreement both in terms of heat fluxes as well as flow structure. Notable flow structures include wedge-shaped regions of high heat flux which emanate from the windward corners of the cube. Analysis of numerical Q-criterion contours show that these high heat flux regions are caused by vortex structures generated by the expansion at the cube corner. Analysis of the numerical skin friction coefficient shows that even at incidence there is no breakaway separation from the expansion edges of the cube and the flow remains attached throughout. We show that although there is little change in the average heat flux experienced by a cube at 5° incidence to the free-stream compared to one at 0° incidence, there are significant changes in the heat flux contours over the cubes at these two incidences. Finally, we calculate a number of heating shape factors which can easily be implemented in satellite re-entry and demise prediction analysis tools.

ANALYSIS OF HISTORICAL BUILDINGS NGAWEN TEMPLE IN MAGELANG

The analysis of historical buildings Ngawen Temple has a purpose of knowing the history of Ngawen Temple and knowing the approach to the Ngawen Temple's geometry. This research type is descriptive research with a qualitative approach. The subject of this research is the historical building of Ngawen Temple. The data collection techniques used in this research are interviews and documentation (observation) directly. The interview was done with the temple officers around the building, and the documentation was taken directly. The data analysis technique is done by analyzing the buildings related to the geometry, such as the building's shape on the Ngawen Temple. The results of this study were obtained; namely, Ngawen Temple set in Buddhism. The Ngawen Temple complex consists of 5 (five) temples that line the parallel from north to south. Temple building facing the east. From the south of Ngawen Temple I, II, III, IV, and V, each temple plans the square. One of Ngawen Temple's uniqueness is the existence of 4 lion statues in every corner of temple II and Candi IV. Another uniqueness is the temple's architecture, found in a lion's statue that supports the four sides of the temple building reconstructed from the five buildings. The lion statue's carved style resembles the lion symbol of Singapore's country and serves the rainwater that comes out through the statue's mouth. The Ngawen Temple's shape resembles the cuboid's geometry, the rectangular pyramidal frustum, and the rectangular pyramid

Experimental characterization of the hypersonic flow around a cuboid

Understanding the hypersonic flow around faceted shapes is important in the context of the fragmentation and demise of satellites undergoing uncontrolled atmospheric entry. To better understand the physics of such flows, as well as the satellite demise process, we perform an experimental study of the Mach 5 flow around a cuboid geometry in the University of Manchester High SuperSonic Tunnel. Heat fluxes are measured using infrared thermography and a 3D inverse heat conduction solution, and flow features are imaged using schlieren photography. Measurements are taken at a range of Reynolds numbers from {40.0 times 10^3} to {549 times 10^3}. The schlieren results suggest the presence of a separation bubble at the windward edge of the cube at high Reynolds numbers. High heat fluxes are observed near corners and edges, which are caused by boundary-layer thinning. Additionally, on the side (off-stagnation) faces of the cube, we observe wedge-shaped regions of high heat flux emanating from the windward corners of the cube. We attribute these to vortical structures being generated by the strong expansion around the cube’s corners. We also observe that the stagnation point of the cube is off-centre of the windward face, which we propose is due to sting flex under aerodynamic loading. Finally, we propose a simple method of calculating the stagnation point heat flux to a cube, as well as relations which can be used to predict hypersonic heat fluxes to cuboid geometries such as satellites during atmospheric re-entry.Graphic abstract

High-Quality MAPbBr3 Cuboid Film with Promising Optoelectronic Properties Prepared by a Hot Methylamine Precursor Approach.

Though CH3NH3PbBr3 single crystals are frequently applied in various optoelectronic devices due to their favorable cuboid geometry, superior optoelectronic properties, and better stability than CH3NH3PbI3, CH3NH3PbBr3 polycrystalline films normally show poorer morphology with scattered crystals than their iodide counterparts, inherently due to their different crystallization habits. In this work, a facile process based on a hot methylamine-based precursor with high viscosity and concentration is demonstrated to counteract rapid ion diffusion. The precursor also has special features including a large colloidal size, a solid form at room temperature, and fast crystallization offered by the easy evacuation of methylamine. CH3NH3PbBr3 films composed of tightly aligned CH3NH3PbBr3 cuboids on micron scale are obtained. Wide channel (100 μm) photodetectors made from the CH3NH3PbBr3 films show promising photoresponse and fast response speeds on par with those based on single crystals, suggesting high film quality and good optoelectronic connections between neighboring cuboids.

2-D DOA Estimation in a Cuboid Array Based on Metaheuristic Algorithms and Maximum Likelihood

This paper proposes to apply the genetic algorithm and the firefly algorithm to enhance the estimation of the direction of arrival (DOA) angle of electromagnetic signals of a smart antenna array. This estimation is essential for beamforming, where the antenna array radiating pattern is steered to provide faster and reliable data transmission with increased coverage. This work proposes using metaheuristics to improve a maximum likelihood DOA estimator for an antenna array arranged in a uniform cuboidal geometry. The DOA estimation performance of the proposed algorithm was compared to that of MUSIC on different two dimensions scenarios. The metaheuristic algorithms present better performance than the well-known MUSIC algorithm.

Boric acid as converter material for semiconductor neutron detector- A GEANT4 simulation study

Monte Carlo based GEANT4 simulations have been carried out to estimate the efficiency (η) of thermal neutron detection with boric acid (H3BO3) as a converter material. Different geometries of detector configurations — planar stack, sphere, cylinder and cuboid were considered and their efficiency were computed for relevant geometric parameters along with different 10B enrichment (BE) levels content in boric acid, with Low Level Discriminator (LLD) value fixed at 300 keV. The planar stack configuration detector was designed by replicating two dimensional detector units (DUs) – each consisting of boric acid as a converter material and Silicon (Si) as a detector material – in the third dimension. The simulation carried out for varied layer thicknesses (0.25μm to 10μm) of boric acid showed that efficiency (η) increased with increasing thickness, reached a maximum at a critical thickness (tsc) and thereafter found to decrease. Typically, at a critical thickness of ∼3.5μm, the η for 100% BE content in boric acid was found to be 15.96% for 30 DUs. For the case of sphere, cylinder and cuboid of boric acid embedded in Si, the simulations were performed for diameter and width both varying from 0.5μm to 9.5μm in steps of 0.5μm. In particular for cylinder and cuboid geometries, the depth was also varied from 25μm to 275μm with a step size of 25μm. In all these configurations η was found to increase initially and was maximum at a critical diameter and width of ∼8–9μm, and thereafter it was found to reduce. Apart from the simulated η estimation, histogram plots depicting the energy deposited in Si detector region by the generated charged particles (7Li = 0.839 MeV and α=1.47MeV) upon thermal neutron interaction in boric acid have also been studied for various geometrical parameters such as thickness, diameter and width of detector.

Development and clinical implementation of eclipse scripting-based automated patient-specific collision avoidance software.

PurposeIncreased use of Linac‐based stereotactic radiosurgery (SRS), which requires highly noncoplanar gantry trajectories, necessitates the development of efficient and accurate methods of collision detection during the treatment planning process. This work outlines the development and clinical implementation of a patient‐specific computed tomography (CT) contour‐based solution that utilizes Eclipse Scripting to ensure maximum integration with clinical workflow.MethodsThe collision detection application uses triangle mesh structures of the gantry and couch, in addition to the body contour of the patient taken during CT simulation, to virtually simulate patient treatments. Collision detection is performed using Binary Tree Hierarchy detection methods. Algorithm accuracy was first validated for simple cuboidal geometry using a calibration phantom and then extended to an anthropomorphic phantom simulation by comparing the measured minimum distance between structures to the predicted minimum distance for all allowable orientations. The collision space was tested at couch angles every 15° from 90 to 270 with the gantry incremented by 5° through the maximum trajectory. Receiver operating characteristic curve analysis was used to assess algorithm sensitivity and accuracy for predicting collision events. Following extensive validation, the application was implemented clinically for all SRS patients.ResultsThe application was able to predict minimum distances between structures to within 3 cm. A safety margin of 1.5 cm was sufficient to achieve 100% sensitivity for all test cases. Accuracy obtained was 94.2% with the 5 cm clinical safety margin with 100% true positive collision detection. A total of 88 noncoplanar SRS patients have been currently tested using the application with one collision detected and no undetected collisions occurring. The average time for collision testing per patient was 2 min 58 s.ConclusionsA collision detection application utilizing patient CT contours was developed and successfully clinically implemented. This application allows collisions to be detected early during the planning process, avoiding patient delays and unnecessary resource utilization if detected during delivery.

Thermal Modeling in Metal Additive Manufacturing Using Graph Theory

Abstract The goal of this work is to predict the effect of part geometry and process parameters on the instantaneous spatiotemporal distribution of temperature, also called the thermal field or temperature history, in metal parts as they are being built layer-by-layer using additive manufacturing (AM) processes. In pursuit of this goal, the objective of this work is to develop and verify a graph theory-based approach for predicting the temperature distribution in metal AM parts. This objective is consequential to overcome the current poor process consistency and part quality in AM. One of the main reasons for poor part quality in metal AM processes is ascribed to the nature of temperature distribution in the part. For instance, steep thermal gradients created in the part during printing leads to defects, such as warping and thermal stress-induced cracking. Existing nonproprietary approaches to predict the temperature distribution in AM parts predominantly use mesh-based finite element analyses that are computationally tortuous—the simulation of a few layers typically requires several hours, if not days. Hence, to alleviate these challenges in metal AM processes, there is a need for efficient computational models to predict the temperature distribution, and thereby guide part design and selection of process parameters instead of expensive empirical testing. Compared with finite element analyses techniques, the proposed mesh-free graph theory-based approach facilitates prediction of the temperature distribution within a few minutes on a desktop computer. To explore these assertions, we conducted the following two studies: (1) comparing the heat diffusion trends predicted using the graph theory approach with finite element analysis, and analytical heat transfer calculations based on Green’s functions for an elementary cuboid geometry which is subjected to an impulse heat input in a certain part of its volume and (2) simulating the laser powder bed fusion metal AM of three-part geometries with (a) Goldak’s moving heat source finite element method, (b) the proposed graph theory approach, and (c) further comparing the thermal trends predicted from the last two approaches with a commercial solution. From the first study, we report that the thermal trends approximated by the graph theory approach are found to be accurate within 5% of the Green’s functions-based analytical solution (in terms of the symmetric mean absolute percentage error). Results from the second study show that the thermal trends predicted for the AM parts using graph theory approach agree with finite element analyses, and the computational time for predicting the temperature distribution was significantly reduced with graph theory. For instance, for one of the AM part geometries studied, the temperature trends were predicted in less than 18 min within 10% error using the graph theory approach compared with over 180 min with finite element analyses. Although this paper is restricted to theoretical development and verification of the graph theory approach, our forthcoming research will focus on experimental validation through in-process thermal measurements.

A computational framework for cortical microtubule dynamics in realistically shaped plant cells.

Plant morphogenesis is strongly dependent on the directional growth and the subsequent oriented division of individual cells. It has been shown that the plant cortical microtubule array plays a key role in controlling both these processes. This ordered structure emerges as the collective result of stochastic interactions between large numbers of dynamic microtubules. To elucidate this complex self-organization process a number of analytical and computational approaches to study the dynamics of cortical microtubules have been proposed. To date, however, these models have been restricted to two dimensional planes or geometrically simple surfaces in three dimensions, which strongly limits their applicability as plant cells display a wide variety of shapes. This limitation is even more acute, as both local as well as global geometrical features of cells are expected to influence the overall organization of the array. Here we describe a framework for efficiently simulating microtubule dynamics on triangulated approximations of arbitrary three dimensional surfaces. This allows the study of microtubule array organization on realistic cell surfaces obtained by segmentation of microscopic images. We validate the framework against expected or known results for the spherical and cubical geometry. We then use it to systematically study the individual contributions of global geometry, cell-edge induced catastrophes and cell-face induced stability to array organization in a cuboidal geometry. Finally, we apply our framework to analyze the highly non-trivial geometry of leaf pavement cells of Arabidopsis thaliana, Nicotiana benthamiana and Hedera helix. We show that our simulations can predict multiple features of the microtubule array structure in these cells, revealing, among others, strong constraints on the orientation of division planes.

Risk Assessment Integrated QbD Approach for Development of Optimized Bicontinuous Mucoadhesive Limicubes for Oral Delivery of Rosuvastatin.

Statins are widely prescribed for hyperlipidemia, cancer, and Alzheimer's disease but are facing some inherent challenges such as low solubility and drug loading, higher hepatic metabolism, as well as instability at gastric pH. So, relatively higher circulating dose, required for exerting the therapeutic benefits, leads to dose-mediated severe toxicity. Furthermore, due to low biocompatibility, high toxicity, and other regulatory caveats such as product conformity, reproducibility, and stability of conventional formulations as well as preferentially higher bioabsorption of lipids in their favorable cuboidal geometry, enhancement in in vivo biopharmaceutical performance of Rosuvastatin could be well manifested in Quality by Design (QbD) integrated cuboidal-shaped mucoadhesive microcrystalline delivery systems (Limicubes). Here, quality-target-product-profile (QTPPs), critical quality attributes (CQAs), Ishikawa fishbone diagram, and integration of risk management through risk assessment matrix for failure mode and effects analysis (FMEA) followed by processing of Plackett-Burman design matrix using different statistical test for the first time established an approach to substantiate the claims that controlling levels of only these three screened out independent process variables, i.e., Monoolein (B = 800-1100μL), Poloxamer (C = 150-200mg), and stirring speed (F = 700-1000rpm) were statistically significant to modulate and improve the biopharmaceutical performance affecting key attributes, viz., average particle size (Y1 = 1.40-2.70μ), entrapment efficiency (Y2 = 62.60-88.80%), and drug loading (Y3 = 0.817-1.15%), in QbD-enabled process. The optimal performance of developed Limicubes exhibited an average particle size of 1.8 ± 0.2μ, entrapment efficiency 80.32 ± 2.88%, and drug loading 0.93 ± 0.08% at the level of 1100μL (+ 1), 200mg (+ 1), and 700rpm (- 1) for Monoolein, Poloxamer, and stirring speed, respectively.

Wall slipping behavior of foam with nanoparticle-armored bubbles and its flow resistance factor in cracks

In this work, wall slipping behavior of foam with nanoparticle-armored bubbles was first studied in a capillary tube and the novel multiphase foam was characterized by a slipping law. A crack model with a cuboid geometry was then used to compare with the foam slipping results from the capillary tube and also to evaluate the flow resistance factor of the foam. The results showed that the slipping friction force FFR in the capillary tube significantly increased by addition of modified SiO2 nanoparticles, and an appropriate power law exponents by fitting FFR vs. Capillary number, Ca, was 1/2. The modified nanoparticles at the surface were bridged together and formed a dense particle “armor” surrounding the bubble, and the interconnected structures of the “armor” with strong steric integrity made the surface solid-like, which was in agreement with the slip regime associated with rigid surface. Moreover, as confirmed by 3D microscopy, the roughness of the bubble surface increased with nanoparticle concentration, which in turn increased the slipping friction force. Compared with pure SDBS foam, SDBS/SiO2 foam shows excellent stability and high flow resistance in visual crack. The resistance factor of SiO2/SDBS foam increased as the wall surface roughness increased in core cracks.

Static and dynamic compression strength of hot-pressed boron carbide using a dumbbell-shaped specimen

Dumbbell-shaped specimens were machined from a tile of a commercially available boron carbide and subjected to static and dynamic compressive loads. Static experiments were performed using a screw-driven load frame, and the fracture process was recorded with a high-speed camera. Dynamic experiments were performed using a split-Hopkinson pressure bar and an ultra-high-speed camera to record the fracture process. The average static (~10−1 s−1) compressive strength of this boron carbide was determined to be 6.1 ± 0.3 GPa, and the dynamic compression strength, up to a strain rate of 300 s−1, was 6.2 ± 0.3 GPa. These strength values are twice the value reported by the manufacturer and approximately 30% higher than the values reported by others for the same boron carbide but when using a cuboidal specimen geometry.

A multiple anisotropic material finite element and its application to stiffness tensor predictions for woven composite lamina

The effective anisotropic stiffness tensor Cijkl of a woven composite lamina can be calculated by performing a finite element analysis (FEA) of the representative volume element (RVE); however, the complex geometry of heterogeneous RVEs makes meshing tedious. This paper develops two unique multiphase voxel elements (MVEs) that account for multiple materials within an elemental domain by utilizing the material properties at each integration point to develop local strain corrections. Studies performed on simple cuboid geometries show exceptional agreement with standard FEA results, and results from the new voxel elements are more accurate than previous MVEs presented in literature. The application for these new MVEs is demonstrated on plain and twill weave composite laminae, and the results generated are in exceptional agreement with experimental results available in the literature and computationally expensive studies using standard finite elements.