Joint-Embedding Predictive Architectures (JEPAs) provide a framework for learning compact representations without directly reconstructing high-dimensional observations. However, in parameterized physical systems, learned representations can entangle geometry with operating conditions and task-specific physical responses, limiting their reuse across prediction tasks. We introduce D-JEPA (Design-recoverable JEPA), a geometry-centric JEPA that computes a compact representation from geometry alone and reuses it across operating conditions and physical response spaces through lightweight physics-specific decoders. An explicit design-recoverability objective encourages the geometry latent to preserve information about the underlying design variables, enabling the representation to support design analysis and optimization. We further identify a case-level collapse failure mode in which target representations become nearly invariant across distinct geometries despite low reconstruction error, and mitigate it using case-level variation constraints and auxiliary target reconstruction. Across four 3D aerodynamic, hydrodynamic, and structural benchmarks, D-JEPA maintains or improves full-field prediction accuracy while achieving near-perfect linear recoverability of design parameters. The frozen geometry representation can be reused at held-out operating conditions and transferred to a structural response task with fewer trainable parameters. Finally, the representation supports differentiable design optimization, with designs validated using high-fidelity CFD, preserving the predicted ranking of candidate designs. These results demonstrate that separating a reusable geometry representation from physics-specific prediction provides a practical representation for scientific surrogate modeling and design.
Figures & tables
Figure 1: Overview of D-JEPA: (a) D-JEPA architecture and training pathways. (b) Design recoverability, condition reuse and cross-physics transfer.
Model
Dataset
Pressure
WSS x
WSS y
WSS z
CD
CL
DoMINO
SHIFT-Wing
0.9678 ± 0.001
0.8972 ± 0.001
0.9456 ± 0.004
0.9654 ± 0.002
0.940 ± 0.002
0.891 ± 0.001
SHIFT-Submarine
0.9321 ± 0.000
0.9332 ± 0.000
0.9421 ± 0.003
0.9243 ± 0.002
0.973 ± 0.001
–
SHIFT-SUV
0.9765 ± 0.000
0.9475 ± 0.000
0.9241 ± 0.002
0.9342 ± 0.002
0.980 ± 0.003
–
GeoTransolver
SHIFT-Wing
0.9950 ± 0.001
0.9350 ± 0.001
0.9645 ± 0.002
0.9899 ± 0.001
0.960 ± 0.002
0.956 ± 0.001
SHIFT-Submarine
0.9921 ± 0.001
0.9432 ± 0.001
0.9621 ± 0.001
0.9743 ± 0.002
0.990 ± 0.002
–
SHIFT-SUV
0.9931 ± 0.001
0.9698 ± 0.001
0.9641 ± 0.001
0.9542 ± 0.001
0.996 ± 0.002
–
Table 1: Surrogate accuracy across parameterized 3D aerodynamic benchmarks. Full-field R2 is reported for surface pressure and wall-shear-stress components, together with integrated force coefficients CD and CL . Values are mean ± std over three seeds; bold indicates the best value for each dataset and quantity.
Dataset
None
Token SIGReg
Case SIGReg
+ Variance
+ Target Recon
SHIFT-Wing
0.0003
0.0003
0.6211
1.0346
1.0317
SHIFT-Submarine
0.0004
0.0005
0.7566
1.0633
1.0620
SHIFT-SUV
0.0002
0.0002
0.0003
0.9892
0.9889
Table 2: Case-level variation in the target representation under progressively stronger representation constraints. Values report stdcase(meanmZt) , measuring variation across geometry cases.
Design R 2 ↑
Pressure R 2 ↑
Dataset
Full
−Llat
Learned-only
−Ldis
Full
−Llat
Learned-only
−Ldis
SHIFT-Wing
0.998
0.993
0.976
0.722
0.9971
0.9912
0.9885
0.9823
SHIFT-Submarine
0.998
0.982
0.757
0.564
0.9967
0.9896
0.9852
0.9872
SHIFT-SUV
0.989
0.984
0.976
0.622
0.9937
0.9851
0.9801
0.9821
Table 3: Component ablation of D-JEPA across the SHIFT datasets, showing the contributions of latent prediction and design supervision to design recoverability and pressure prediction. Learned-only removes the parameter free global and regional descriptors.
Trained operating conditions
Held-out
Mach
0.5000
0.6000
0.7000
0.7500
0.8000
0.8250
0.8375
0.8500
0.8750
0.9000
0.7750
0.8125
Pressure R2
0.9902
0.9935
0.9948
0.9944
0.9941
0.9956
0.9960
0.9967
0.9952
0.9939
0.9922
0.9944
Std
0.001
0.002
0.001
0.003
0.002
0.001
0.001
0.002
0.001
0.002
0.001
0.002
Table 4: SHIFT-Wing pressure prediction accuracy ( R2 ) across trained and held-out Mach numbers using a single frozen geometry representation g . The held-out Mach values are interpolation points within the trained Mach range. Bold values denote held-out conditions.
Training epoch budget
Frozen geometry encoder (0.56 M params)
From scratch (2.30 M params)
10
0.3316 ± 0.020
0.0396 ± 0.007
100
0.8901 ± 0.013
0.8740 ± 0.034
200
0.9327 ± 0.002
0.9297 ± 0.012
300
0.9385 ± 0.001
0.9313 ± 0.002
Table 5: Structural displacement prediction accuracy ( R2 ) on AutoHood 3D comparing frozen geometry encoder transfer versus training from scratch.
Benchmark
Design candidate
Predicted [N]
CFD [N]
Relative error
Predicted → CFD rank
SHIFT-Submarine (drag)
OPT-wide
645.8
637.69
1.3%
1 → 1
OPT-tau
677.0
688.36
1.7%
2 → 2
MID
718.6
723.08
0.6%
3 → 3
ANTI
758.7
773.75
2.0%
4 → 4
SHIFT-Wing (lift)
OPT-Wide
3,426,080
3,193,866
7.27%
1 → 1
OPT-Tau
3,108,536
3,086,237
0.72%
2 → 2
Table 6: High-fidelity CFD validation of surrogate-optimized designs. Integrated drag is reported for SHIFT-Submarine and integrated lift for SHIFT-Wing. Predicted and CFD rankings agree for all evaluated designs.
Appendix figures & tables11 assets
Supplementary material from the paper’s appendix.
Appendix
Dataset
Cases
Design parameters
Operating conditions
Surface fields
SHIFT-Wing
1,966
7
angle of attack, Mach
pressure, WSS x , WSS y , WSS z
SHIFT-Submarine
4,000
7
fixed
pressure, WSS x , WSS y , WSS z
SHIFT-SUV
1,996
7
fixed
pressure, WSS x , WSS y , WSS z
AutoHood 3D
12,000
–
load case
pressure, displacement
Appendix
Table A.1: Benchmark dataset properties, parameter dimensions, operating conditions and surface targets. A dash indicates that the quantity is not applicable.
Loss term
Latent
Recon.
Design
Token
Case
Var.
Tgt. rec.
( Llat )
( λrLrec )
( λdLdis )
( λsLsig )
( λcLsig-case )
( λvLvar )
( λtLrec-t )
Weight
1.0
1.0
1.0
0.01
1.0
1.0
1.0
Appendix
Table A.2: Hyperparameter weights used for the composite D-JEPA objective.
Setting
Value
Optimizer
AdamW ( Loshchilov and Hutter, 2019 )
Base LR
1×10−3
Weight decay
1×10−3
LR schedule
Cosine decay
Gradient clipping
10.0
Epochs
400
Appendix
Table A.3: Optimization configuration and training hyperparameters.
Dataset
Configuration
σcase
Llat
Lrec
Pressure R2
Design R2
SHIFT-Wing
None
0.0003
0.0000
0.0354
0.9906
0.992
Token SIGReg
0.0003
0.0000
0.0352
0.9906
0.992
Case SIGReg
0.6211
0.0403
0.0373
0.9893
0.987
+ Variance
1.0346
0.0498
0.0381
0.9886
0.986
+ Target Recon.
1.0317
0.0447
0.0379
0.9971
0.998
SHIFT-Submarine
None
0.0004
0.0000
0.0348
0.9896
0.982
Appendix
Table A.4: Complete representation-collapse diagnostics. The case-level standard deviation measures between-case variation, while Llat and Lrec report latent prediction and physical reconstruction losses. Pressure R2 and design-recovery R2 quantify downstream performance.
Figure A.1: Qualitative surface pressure prediction for SHIFT-SUV. Predicted pressure and absolute point-wise error are shown for representative median and worst-case held-out geometries.
Figure A.2: Qualitative surface pressure prediction for SHIFT-Submarine. Predicted pressure and absolute point-wise error are shown for representative median and worst-case held-out geometries.
Figure A.3: Qualitative surface pressure prediction for SHIFT-Wing. Predicted pressure and absolute point-wise error are shown for representative median and worst-case held-out geometries.
SHIFT-SUV
SHIFT-Wing
SHIFT-Submarine
Parameter
R2
Parameter
R2
Parameter
R2
hood_height
0.995
AR (aspect ratio)
0.999
cross_section_shape
1.000
hood_base_height
0.987
Lambda_c_4 (quarter-chord sweep)
0.998
fore_body_shape
0.996
windshield_angle
0.986
twist_20
0.997
fore_length
0.999
vehicle_height
0.991
twist_k
0.997
mid_length
0.999
backlight_angle
0.990
twist_60
0.999
aft_tail_half_angle
0.997
Appendix
Table A.5: Linear-probe recovery accuracy ( R2 ) for individual geometric design parameters from the frozen 128-dimensional geometry representation g .
Figure A.4: Design-variable probe recoverability across benchmark datasets. True design parameters are compared with values recovered from the frozen geometry representation using linear probes.
Figure A.5: Operating-condition generalization across Mach numbers on SHIFT-Wing. Surface pressure R2 is evaluated across trained and held-out Mach conditions using the same geometry representation.
Execution pipeline stage
Compute hardware
Wall-clock time
D-JEPA/AeroJEPA training (400 epochs)
1 × NVIDIA T4 (16 GB)
4–8 hours / dataset
Decoder transfer fine-tuning (400 epochs)
1 × NVIDIA T4 (16 GB)
∼ 2 hours
GeoTransolver baseline (500 epochs)
8 × NVIDIA H100 (80 GB)
∼ 12 hours / dataset
Mach-condition fine-tuning (200 epochs)
1 × NVIDIA T4 (16 GB)
∼ 3 hours
Appendix
Table A.6: Primary computational hardware configurations and approximate wall-clock runtimes across pipeline stages.