https://community.wolfram.com RSS Feed for Wolfram Community showing any discussions tagged with Chemical Engineering sorted by active. https://community.wolfram.com/groups/-/m/t/3647332 &[Wolfram Notebook][1] [1]: https://www.wolframcloud.com/obj/99571eb2-076b-4c46-b02a-35f7bf96c8b1 Housam Binous 2026-02-28T21:02:03Z https://community.wolfram.com/groups/-/m/t/3524866 ![Thermal runaway in validation of industrial manufacturing with System Modeler and Wolfram Language][1] &[Wolfram Notebook][2] [1]: https://community.wolfram.com//c/portal/getImageAttachment?filename=thermal.png&userId=1618777 [2]: https://www.wolframcloud.com/obj/75063c23-f701-4697-9062-6e529766d792 Sergio Vargas 2025-08-05T07:15:14Z https://community.wolfram.com/groups/-/m/t/3445315 &[Wolfram Notebook][1] [1]: https://www.wolframcloud.com/obj/c4e6f46b-4ff2-43eb-af32-f17e8d2133d2 Housam Binous 2025-04-16T14:33:22Z https://community.wolfram.com/groups/-/m/t/3445135 &[Wolfram Notebook][1] [1]: https://www.wolframcloud.com/obj/ef543934-0d38-429b-904c-e0db22fe8971 Housam Binous 2025-04-16T09:19:20Z https://community.wolfram.com/groups/-/m/t/3345034 &[Wolfram Notebook][1] [1]: https://www.wolframcloud.com/obj/1e682589-cfa7-41db-bcba-3d4e6e3cbd2d Housam Binous 2024-12-26T07:27:44Z https://community.wolfram.com/groups/-/m/t/3103297 &[Wolfram Notebook][1] [1]: https://www.wolframcloud.com/obj/1b94967d-26e5-4e19-b5a3-a4e816ac6d8b Housam Binous 2024-01-16T13:48:53Z https://community.wolfram.com/groups/-/m/t/2904959 &[Wolfram Notebook][1] [1]: https://www.wolframcloud.com/obj/382d3008-5f34-4cd4-b3ae-01c2e45d9419 Housam Binous 2023-04-23T22:31:31Z https://community.wolfram.com/groups/-/m/t/2904950 &[Wolfram Notebook][1] [1]: https://www.wolframcloud.com/obj/f61d6e54-b547-479b-ac48-de9482cc4f22 Housam Binous 2023-04-23T22:25:28Z https://community.wolfram.com/groups/-/m/t/2904594 &[Wolfram Notebook][1] [1]: https://www.wolframcloud.com/obj/1b4a53d1-a869-48ff-b32a-a28c1079317f Housam Binous 2023-04-23T19:42:53Z https://community.wolfram.com/groups/-/m/t/2904572 &[Wolfram Notebook][1] [1]: https://www.wolframcloud.com/obj/d325269f-6f59-4f86-b094-df966afda380 Housam Binous 2023-04-23T17:12:21Z https://community.wolfram.com/groups/-/m/t/2904561 &[Wolfram Notebook][1] [1]: https://www.wolframcloud.com/obj/80bcb33c-4de6-4847-9713-e04d57fc43dd Housam Binous 2023-04-23T16:40:28Z https://community.wolfram.com/groups/-/m/t/2904551 &[Wolfram Notebook][1] [1]: https://www.wolframcloud.com/obj/23eff803-c560-4ed0-b97a-af2bc9fc4da6 Housam Binous 2023-04-23T16:23:29Z https://community.wolfram.com/groups/-/m/t/2904536 &[Wolfram Notebook][1] [1]: https://www.wolframcloud.com/obj/7b5fa2e1-5547-47e2-be55-138c2421f0a5 Housam Binous 2023-04-23T15:35:10Z https://community.wolfram.com/groups/-/m/t/2904467 &[Wolfram Notebook][1] [1]: https://www.wolframcloud.com/obj/074f1769-fcee-4968-8d66-7dfc58b9e4d6 Housam Binous 2023-04-23T13:08:32Z https://community.wolfram.com/groups/-/m/t/2904413 &[Wolfram Notebook][1] [1]: https://www.wolframcloud.com/obj/d3784905-d9fb-41dc-96c5-772a7c71a41c Housam Binous 2023-04-23T12:04:50Z https://community.wolfram.com/groups/-/m/t/2030260 ![enter image description here][1] &[Wolfram Notebook][2] [1]: https://community.wolfram.com//c/portal/getImageAttachment?filename=CapturadeTela2020-11-20a%CC%80s18.42.07.png&userId=2029333 [2]: https://www.wolframcloud.com/obj/64ef0cca-d762-4a13-81fa-609012f78f17 Gino Andrade 2020-07-14T17:20:46Z https://community.wolfram.com/groups/-/m/t/2405785 Hi ExternalMedia and CoolProp is a Modelica third party library that can be compiled and used within Modelica 3.2.3 Has anyone used this within Wolfram SystemModeler and is there instructions for installing the library. Is it just the case of compiling the library and then loading the package? Anyone with success and advice for this would be greatly appreciated. Thanks. Craig Branch 2021-11-13T11:11:08Z https://community.wolfram.com/groups/-/m/t/2405710 ![Typical Case][1] Of course I had to lead with an image. This is a subset of cellular automata, with extra rules. The rules are procedural, and scalable. I get into more detail in the attached notebooks, but the gist of it is the rules are designed to produce universes which are comprehensible to a human, have (some) conservation of energy, have an energy input, and have a way to corral that energy in one spot. In the image, the gray boxes on the ground are just that - non-interacting ground. The colored boxes are cells with mass energy (I've been calling them hadrons), and the pips are light energy. The star moves across the screen, leaving more hadrons in its wake, which then interact. There's a gravity bias, such that it costs light energy to move mass upwards. There's a cap on the amount of energy any cell can have, and cells with more get trimmed back to the cap, which allows for defining the amount of colors. It's inefficient to generate rule numbers for this, which I've gotten around by using a random seed value for the uncalculated parts of the rule. This guarantees sharing a universe with the same seed will result in the same rule number. The code is mostly written for displaying 2D cellular automata with small caps of energy, with these specific special star and ground cells. However, a lot of it should be easily modifiable to work for 3D or above, or with other concepts for special cells. As far as applications - maybe this is a ways out, but I'd like to see if we can study how life could develop in another universe. I'd be happy with finding stable complicated configurations analogous to biochemicals in the meantime. Hence me flagging this post for chemistry & bio sciences - this is your wheelhouse. How to use the notebooks: There are two notebooks attached: an "Instructional" notebook with copious notes, and an "Operational" notebook with far fewer notes. There is also a "Bad Rule Digits.txt" file which stores data for a workaround regarding errors in the CellularAutomata function that may or may not be specific to my machine. You can place the .badlist file in your Documents folder, under a folder (which the code may create automatically) called "Star and Ground CA". The notebooks will also store screenshots of every step in this folder, along with it being the save location for the universe information, currently known rule number, and instance information. I've done about as much checking as sanity permits that everything works, but if you find errors or have questions, please let me know. PS: I was unable to attach a .badlist file until I changed the extension to .txt. It's not critical to have this file, but change the extension back to .badlist to get the notebook to read it. Thanks! [1]: https://community.wolfram.com//c/portal/getImageAttachment?filename=FirstUniverse211112140824.png&userId=2392475 Cameron Kosina 2021-11-12T23:27:05Z https://community.wolfram.com/groups/-/m/t/2369818 &[Wolfram Notebook][1] [1]: https://www.wolframcloud.com/obj/5e6a950d-67eb-4dde-98b7-db38b97d647d Maher Kuzbari 2021-09-20T09:29:40Z https://community.wolfram.com/groups/-/m/t/1378744 ![enter image description here][8] According to the World Health Organization(WHO), more than a quarter billion people today are living with Hepatitis B virus (HBV) and about 890,000 die every year as a result. As such, viruses are responsible for a variety of illnesses, ranging from common cold to some type of cancer and more serious conditions such as AIDS. Although the symptoms of infections might be different, the methods of infection are quite similar. The virus enters the cell, typically through the interaction with a receptor on the cell surface; it must then disassemble its capsid to release the genetic material. The viral genome is then replicated and the proteins are synthesized; new viruses can then assemble, and exit the cell either through the membrane or when the cell dies. With the capsid, the viral protein shell, viruses conceal their genetic material and can slip inside the host's body (like a Trojan horse!). The goal of my project was to create an applicable structural model of the proteins that compose the capsid to better understand the viral life cycle and infection mechanism. First, import data from protein database(PDB): data =Import[URLDownload["https://files.rcsb.org/download/3J3I.pdb1.gz"],"String"]; Since the data arrangement is irregular, replace negative signs(-) with spaces to better align and separate columns: spaced = StringReplace[data,"-"-> " -"] Read data starting from row 180 where atomic information starts: Do[Read[str,"String"],181]; Read[str,"String"] Split the list when consecutive first letter appears to separate between subunits: splitted = Split[fullrecord,First[#1]===First[#2]&]; Remove unnecessary rows(only take ATOM data): removed = DeleteCases[splitted,x_/;First[First[x]]== "TER"]; removed = DeleteCases[removed,x_/;First[First[x]]== "HETATM"]; removed = DeleteCases[removed,x_/;First[First[x]]== "MODEL"]; removed = DeleteCases[removed,x_/;First[First[x]]== "ENDMDL"]; removed = DeleteCases[removed,x_/;First[First[x]]== "MASTER"]; removed = DeleteCases[removed,x_/;First[First[x]]== "END"]; Now check dimensions of the data, the dimension should equal 60 times the T number of the virus structure that is used. For example, 3J3I virus that is used in this code as an example is T = 3, therefore the dimensions should equal 180. This is because of the 60 fold symmetry of icosahedral viruses. Dimensions[removed] Now with the data polished, we need to extract specific coordinates and chain types: xcoor = ToExpression[removed[[All,All, 7]]] ycoor = ToExpression[removed[[All,All, 8]]] zcoor = ToExpression[removed[[All,All, 9]]] atom = removed[[All,All, -1]]; To calculate center of mass, we need to replace chain types with corresponding atomic mass: atomicmasses = Association@ Thread[{"N", "C", "O", "S"} -> Through@{Entity["Element", "Nitrogen"], Entity["Element", "Carbon"], Entity["Element", "Oxygen"], Entity["Element", "Sulfur"]}[ EntityProperty["Element", "AtomicMass"]]] Now we can calculate center of mass: CoM[mass_, pos_] := mass.pos/Total[mass] Put all x, y, z center of mass coordinates together for convenience: totcentercoor[xset_,yset_,zset_] :=List[xset, yset, zset] centercoor = MapThread[totcentercoor, {xset, yset, zset}]; totcoor[x_,y_,z_] :=List[x, y, z] allcoor = MapThread[totcoor, {x, y, z}]; With all coordinate and chain type data, create a 3D point cloud with all protein units. For better analysis, chain types are matched with different colors. colorlist = ToExpression[flatchain/.{"A"-> "Blue", "B" -> "Pink", "C" -> "Green", "D" -> "Red"}]; Graphics3D[{Point[allcoor,VertexColors->colorlist]}, Lighting->Automatic, Boxed->False] And we get a cool image of the entire capsid structure! (This proves my point that all good science is art - look how beautiful the capsid structure is!) ![point cloud][1] Let's also generate a point cloud with center of mass and see how it looks: Graphics3D[{PointSize[0.02],Point[centercoor,VertexColors->colorlist2]}] ![center of mass][2] We now see how the center of mass is distributed and look at each subunits. Just to make sure we're doing everything right, create a surface mesh to see the general structure: ListSurfacePlot3D[centercoor, BoxRatios->{1, 1, 1}] ![mesh][3] Now comes the fun part, we can start tiling on the structure! To create an edge-to-edge tiling on a spherical surface, we are going to use Penrose tiling, an example of non-periodic tiling generated by an aperiodic set of prototiles. Because Penrose tiling is non-periodic, there is no translational symmetry. However it is self similar, which means the same patterns occur at larger and larger scales. Thus, the tiling can be obtained through "inflation" (or "deflation") and every finite patch from the tiling occurs infinitely many times. Such tiling method perfectly matches our purpose of modeling viral capsids because Penrose tiling can be constructed to exhibit both reflection symmetry and fivefold rotational symmetry. We start with creating a Delaunay mesh on the structure with all center of mass points as vertices of triangles. Delaunay triangulations maximize the minimum angle of all the angles of the triangles in the triangulation, thus avoiding silver triangles, triangles with one or two extremely acute angles. Thankfully, Mathematica has a built-in function for creating Delaunay mesh making our job a lot easier. cm =RegionBoundary@DelaunayMesh@centercoor ![enter image description here][4] However, the mesh only creates sets of triangles. To get a uniform tiling, we need to figure out which triangles are most likely to form a subunit with another. Another problem arises with this method because not all triangles are on the same plane, making it harder to classify the triangles. To solve this problem, we are going to use some linear algebra and calculate the normal vectors of each triangles. All triangles with relatively similar normal vectors will be classified together and grouped on the same plane. normal[Polygon[{p1_,p2_,p3_}]]:=Normalize@Cross[p1-p2,p1-p3] delaunaymesh = DelaunayMesh@centercoor; trianglecomponents = MeshPrimitives[RegionBoundary@delaunaymesh,2]; normals=normal/@trianglecomponents; rounded = Round[normals,0.1]; tally = Tally[rounded]; pent = GatherBy[tally,Last][[2]]; quad = GatherBy[tally,Last][[3]]; tri = GatherBy[tally,Last][[1]]; pentcoor = First /@ pent; quadcoor = First /@ quad; tricoor = First /@ tri; value1 = Flatten /@ Map[Position[ rounded, #]&, pentcoor]; value2 = Flatten /@ Map[Position[ rounded, #]&, quadcoor]; value3 = Flatten /@ Map[Position[ rounded, #]&, tricoor]; coors1 = Map[Part[trianglecomponents,#] &, value1]; coors2 = Map[Part[trianglecomponents,#] &, value2]; coors3 = Map[Part[trianglecomponents,#] &, value3]; flat = Flatten[Map[ReplaceAll[Polygon->List][#] &, coors3],2]; subset = Map[Subsets[#, {2}]&,flat]; Since we want to look at the three-fold symmetry of the capsid structure, we are going to extract equilateral triangles(each side represents a rotational axis- therefore an equilateral triangle will be the center of three-fold rotational axis). distance = Round[Apply[EuclideanDistance, subset, {2}],0.1]; equilateral[{x_,y_,z_}]:=(If[x==y==z,Return[1]]) positions = Flatten[ Position[equilateral /@ distance,1],1]; center = Map[Part[coors3,#] &, positions]; Although we grouped the triangles together, graphing them again requires a bit more work since the order of coordinates matter a lot when graphing the polygons. To do this, we will calculate the shortest tour around all vertices and re-orient them in order. extracted[{Polygon[{r1_, r2_, r3_}], Polygon[{t1_, t2_, t3_}], Polygon[{y1_, y2_, y3_}]}] := List[{r1, r2, r3, t1, t2, t3, y1, y2, y3}] pentcoorext = DeleteDuplicates /@ Flatten[extracted /@ coors1, 1]; path1 = Last /@ FindShortestTour /@ pentcoorext; pentfinal = Map[Take[#, 5] &, MapThread[Part, {pentcoorext, path1}]]; extracted2[{Polygon[{r1_, r2_, r3_}], Polygon[{t1_, t2_, t3_}]}] := List[{r1, r2, r3, t1, t2, t3}] quadextract = DeleteDuplicates /@ Flatten[extracted2 /@ coors2, 1]; path2 = Last /@ FindShortestTour /@ quadextract; quadfinal = Map[Take[#, 4] &, MapThread[Part, {quadextract, path2}]]; extracted3[{Polygon[{r1_, r2_, r3_}]}] := List[{r1, r2, r3}] triextract = DeleteDuplicates /@ Flatten[extracted3 /@ coors3, 1]; path3 = Last /@ FindShortestTour /@ triextract; trifinal = Map[Take[#, 3] &, MapThread[Part, {triextract, path3}]]; Now we can create a 3D Penrose tiling on a spherical surface(icosahedral to be specific, but it can be inflated to a sphere)! This is the complete model of simplified subunit structure of the capsid. graphed = Graphics3D[{White, Point[Mean /@ pentcoorext], Blue, EdgeForm[], coors1, Pink, EdgeForm[], coors2, Yellow, EdgeForm[{Thick, Black}], coors3, Green, EdgeForm[], center}, Lighting -> Automatic, Boxed -> False] ![enter image description here][5] But we are not done here yet. To better analyze the structural symmetry, we want to look at the repeated pattern (there would be 20 occurrences of the unit pattern since icosahedral viruses exhibit the identities of icosahedral group) on a flat surface(2D). To do this, we need to flatten all tiles relative to a specific center. We orthogonalize the tiles and find orthogonal vectors relative to the central equilateral triangle we previously extracted. projector[projectingPoint_] := Module[{a, b, c, d, nv, k}, {a, b, c} = Flatten[center[[9]] /. Polygon -> List, 2]; nv = Cross[b - a, c - a]; d = -nv.a; k = (-d - nv.projectingPoint)/Norm[nv]^2; projectingPoint + k*nv ] &; With the vectors calculated, we get a flattened structure of the unit pattern, all aligned on one plane: object = Graphics3D[{ Polygon /@ Map[projector, Part[pentfinal, {1, 2, 3}], {2}], center[[9]], Polygon /@ Map[projector, Part[quadfinal, {1, 2, 8}], {2}], Polygon /@ Map[projector, Part[trifinal, {12, 4, 3, 45, 64, 69}], {2}], Boxed -> False }] Planarize[x_, {u_, v_, w_}] := x.Inverse[Orthogonalize[{v - u, w - u, Cross[v - u, w - u]}]] pentproj = Map[projector, Part[pentfinal, {1, 2, 3}], {2}]; triproj = center[[9, 1, 1]]; quadproj = Map[projector, Part[quadfinal, {1, 2, 8}], {2}]; backproj = Map[projector, Part[trifinal, {12, 4, 3, 45, 64, 69}], {2}]; pentextcoors = Map[Planarize[#, center[[9, 1, 1]]] &, pentproj ]; triextcoors = Map[Planarize[#, center[[9, 1, 1]]] &, triproj ]; quadextcoors = Map[Planarize[#, center[[9, 1, 1]]] &, quadproj ]; backextcoors = Map[Planarize[#, center[[9, 1, 1]]] &, backproj ]; We now see a flattened structure of the unit pattern, with each subunit classified: ![enter image description here][6] One last step to go, we just need to drop the z axis since all coordinates are aligned to put the structure onto a 2D surface. map[x_] := Map[Drop[#, -1] &, x] pentflat = Polygon /@ map /@ pentextcoors; triflat = Polygon@Map[Drop[#, -1] &, triextcoors]; quadflat = Polygon /@ map /@ quadextcoors; backflat = Polygon /@ map /@ backextcoors; Graphics[{Tooltip[ {EdgeForm[{Thick, Black}], Blue, pentflat}, ImageResize[Import[pentimg], {150}]], Tooltip[{EdgeForm[{Thick, Black}], Green, triflat}, "ChainB"], Tooltip[ {EdgeForm[{Thick, Black}], Pink, quadflat}, "ChainB"], Tooltip[ {EdgeForm[{Thick, Black}], Yellow, backflat}, "ChainA"]}] We can now clearly see the unit pattern on a 2D space! Since the pentagons always are on the vertices of the bigger triangle that contain the unit pattern, our structure exhibits all 2,3,5-fold symmetry. ![enter image description here][7] As such, we have analyzed the structural pattern of subunits of viral capsids. The closeness of the obtained idealized tilings to the structures of considered capsids shows that the proposed approach is general for many structures of small and medium-sized viruses. Such tilings provides a general idea about quasi-equivalence of proteins, thus simplify the internal structure of protein molecules and viral genome. Our model provides a better insight into the protein structures of viral capsid, enabling a further research on virus models that do not comply to the conventional Caspar-and-Klug model. Below is the model applied to five different types of icosahedral viruses - 3J3I(T=1), 2VF1(T=2), 1JS9(T=3), 2XFC(T=4), 3J4U(T=7): (Although I did not find the unit patterns for T=4 and T=7 structures, the structural pattern is obvious in the 3D models. ![enter image description here][8] Based on the developed models, I would like to continue researching on viral capsid models, especially on how the stability and structures of different subunits impact the assembly and disassembly pathway. I plan to carry out my analysis on a broader dataset and collect a more generalized modeling algorithm. Author Information: Jihyeon Je(Western Reserve Academy) jej19@wra.net [1]: http://community.wolfram.com//c/portal/getImageAttachment?filename=11.png&userId=1352003 [2]: http://community.wolfram.com//c/portal/getImageAttachment?filename=bb.png&userId=1352003 [3]: http://community.wolfram.com//c/portal/getImageAttachment?filename=cc.png&userId=1352003 [4]: http://community.wolfram.com//c/portal/getImageAttachment?filename=dd.png&userId=1352003 [5]: http://community.wolfram.com//c/portal/getImageAttachment?filename=eee.png&userId=1352003 [6]: http://community.wolfram.com//c/portal/getImageAttachment?filename=ff.png&userId=1352003 [7]: http://community.wolfram.com//c/portal/getImageAttachment?filename=w.png&userId=1352003 [8]: http://community.wolfram.com//c/portal/getImageAttachment?filename=y.png&userId=1352003 Jihyeon Je 2018-07-11T17:06:39Z